Optical-resolution photoacoustic microscopy (OR-PAM) becomes a premier microscopic imaging tool in biomedicine because it provides agent-free optical absorption information in tissues. By tightly focusing light to a spot, a significantly improved lateral resolution can be achieved in OR-PAM. The focal spot size is typically determined by the numerical aperture of the used objective lens. Here, we demonstrate objective-free OR-PAM using a fiber optic Bessel beam generator. In this approach, no objective lens is required and, beneficially, the complexities of conventional OR-PAM systems can be greatly relieved. We have obtained photoacoustic images of a carbon fiber with a diameter of ~ 6 μm, whose lateral resolution was measured to be better than 6 to 7 μm.

A photoacoustic (PA) and ultrasound (US) dual modality system, for imaging human peripheral joints, is introduced. The system utilizes a commercial US unit for both US control imaging and PA signal acquisition. Preliminary in vivo evaluation of the system, on normal volunteers, revealed that this system can recover both the structural and functional information of intra- and extra-articular tissues. Confirmed by the control US images, the system, on the PA mode, can differentiate tendon from surrounding soft tissue based on the endogenous optical contrast. Presenting both morphological and pathological information in joint, this system holds promise for diagnosis and characterization of inflammatory joint diseases such as rheumatoid arthritis.

The morphogenetic relationship between early patterning and polarity formation is of fundamental interest and remains a controversial issue in preimplantation embryonic development. We use a label-free three-dimensional (3-D) imaging technique of full-field optical coherence tomography (FF-OCT) successfully for the first time to study the dynamics of developmental processes in mouse preimplantation lives. Label-free 3-D subcellular time-lapse images are demonstrated to investigate 3-D spatial relationship between the second polar body (2PB) and the first cleavage plane. By using FF-OCT together with quantitative study, we show that only 25% of the predicted first cleavage planes, defined by the apposing plane of two pronuclei, pass through the 2PB. Also only 27% of the real cleavage planes pass through the 2PB. These results suggest that the 2PB is not a convincing spatial cue for the event of the first cleavage. Our studies demonstrate the feasibility of FF-OCT in providing new insights and potential breakthroughs to the controversial issues of early patterning and polarity in mammalian developmental biology.

Light-tissue interactions that influence vascular contrast enhancement in narrow band imaging (NBI) have not been the subject of extensive theoretical study. In order to elucidate relevant mechanisms in a systematic and quantitative manner we have developed and validated a Monte Carlo model of NBI and used it to study the effect of device and tissue parameters, specifically, imaging wavelength (415 versus 540 nm) and vessel diameter and depth. Simulations provided quantitative predictions of contrast-including up to 125% improvement in small, superficial vessel contrast for 415 over 540 nm. Our findings indicated that absorption rather than scattering-the mechanism often cited in prior studies-was the dominant factor behind spectral variations in vessel depth-selectivity. Narrow-band images of a tissue-simulating phantom showed good agreement in terms of trends and quantitative values. Numerical modeling represents a powerful tool for elucidating the factors that affect the performance of spectral imaging approaches such as NBI.

A semi-infinite space Green’s function of the time-dependent radiative transfer equation is derived based upon an exact solution for anisotropic scattering in infinite space and the approximated extrapolated boundary condition. The obtained solution is compared to Monte Carlo simulations and used for retrieving the optical parameters of simulated data in a nonlinear fit. It is shown that the solution performs well for boundaries with mismatched refractive indices. The relative errors in the fitted optical parameters are considerably smaller than by using the diffusion theory.

Surfactant stabilized microbubbles are widely used clinical contrast agents for ultrasound imaging. In this work, the light propagation through a turbid medium in the presence of microbubbles has been investigated. Through a series of experiments, it has been found that the optical attenuation is increased when the microbubbles in a turbid medium are insonified by ultrasound. Such microbubble enhanced optical attenuation is a function of both applied ultrasound pressure and microbubble concentration. To understand the mechanisms involved, a Monte Carlo (MC) model has been developed. Under ultrasound exposure, the sizes of microbubbles vary in space and time, and their dynamics are modeled by the Rayleigh-Plesset equation. By using Mie theory, the spatially and temporally varying optical scattering and scattering efficiency of microbubbles are determined based on the bubble sizes and internal refractive indices. The MC model is shown to effectively describe a medium with rapidly changing optical scattering, and the results are validated against both computational results using an N-layered diffusion equation model and experimental results using a clinical microbubble contrast agent (SonoVue®).

Brain tissue is highly vulnerable to ischemia/hypoxia, and real-time monitoring of its viability is important. By fiber-based measurements for rat brain, we previously observed a unique triphasic reflectance change (TRC) after a certain period of time after hypoxia. After TRC, rats could not be rescued, suggesting that TRC can be used as an indicator of loss of brain tissue viability. In this study, we investigated this diffuse-reflectance change due to hypoxia in three parts. First, we developed and validated a theoretical method to quantify changes in the absorption and reduced scattering coefficients involved in TRC. Second, we performed charge-coupled-device-based reflectance imaging of the rat brain during hypoxia followed by reoxygenation to examine spatiotemporal characteristics of the reflectance and its correlation with reversibility of brain tissue damage. Third, we made simultaneous imaging and fiber-based measurement of the reflectance for the rat to compare signals obtained by these two modalities. We observed a nontriphasic reflectance change by the imaging, and it was associated with brain tissue viability. We found that TRC measured by the fibers preceded the reflectance-signal change captured by the imaging. This time difference is attributable to the different observation depths in the brain with these two methods.

Photoacoustic imaging, using targeted plasmonic metallic nanoparticles, is a promising noninvasive molecular imaging method. Analysis of the photoacoustic signal generated by plasmonic metallic nanoparticles is complex because of the dependence upon physical properties of both the nanoparticle and the surrounding environment. We studied the effect of the aggregation of gold nanoparticles on the photoacoustic signal amplitude. We found that the photoacoustic signal from aggregated silica-coated gold nanoparticles is greatly enhanced in comparison to disperse silica-coated gold nanoparticles. Because cellular uptake and endocytosis of nanoparticles results in their aggregation, these results have important implications for the application of plasmonic metallic nanoparticles towards quantitative molecular imaging.

We achieved fast and automatic retinal vessel segmentation by employing the active contours method in Doppler spectral-domain optical coherence tomography (SD-OCT). In a typical OCT B-scan image, we first extracted the phase variations between adjacent A-lines and removed bulk motion. Then we set the initial contour as the boundary of the whole image and iterated until all of the segmented vessel contours became stabilized. Using a typical office computer, the whole segmentation took no more than 50 s, making real-time retinal vessel segmentation possible. We tested the active contours method segmentation in both controlled phantom and in vivo rodent eye images.

Diffuse fluorescence tomography requires high contrast-to-background ratios to accurately reconstruct inclusions of interest. This is a problem when imaging the uptake of fluorescently labeled molecularly targeted tracers in tissue, which can result in high levels of heterogeneously distributed background uptake. We present a dual-tracer background subtraction approach, wherein signal from the uptake of an untargeted tracer is subtracted from targeted tracer signal prior to image reconstruction, resulting in maps of targeted tracer binding. The approach is demonstrated in simulations, a phantom study, and in a mouse glioma imaging study, demonstrating substantial improvement over conventional and homogenous background subtraction image reconstruction approaches.

Oxidative stress (OS) and mitochondrial dysfunction contribute to photoreceptor cell loss in retinal degenerative disorders. The metabolic state of the retina in a rodent model of retinitis pigmentosa (RP) was investigated using a cryo-fluorescence imaging technique. The mitochondrial metabolic coenzymes nicotinamide adenine dinucleotide (NADH) and flavin adenine dinucleotide (FAD) are autofluorescent and can be monitored without exogenous labels using optical techniques. The cryo-fluorescence redox imaging technique provides a quantitative assessment of the metabolism. More specifically, the ratio of the fluorescence intensity of these fluorophores (NADH/FAD), the NADH redox ratio (RR), is a marker of the metabolic state of the tissue. The NADH RR and retinal function were examined in an established rodent model of RP, the P23H rat compared to that of nondystrophic Sprague-Dawley (SD) rats. The NADH RR mean values were 1.11±0.03 in the SD normal and 0.841±0.01 in the P23H retina, indicating increased OS in the P23H retina. Electroretinographic data revealed a significant reduction in photoreceptor function in P23H animals compared to SD nozrmal rats. Thus, cryo-fluorescence redox imaging was used as a quantitative marker of OS in eyes from transgenic rats and demonstrated that alterations in the oxidative state of eyes occur during the early stages of RP.

The accepted screening technique for colon cancer is white light endoscopy. While most abnormal growths (lesions) are detected by this method, a significant number are missed during colonoscopy, potentially resulting in advanced disease. Missed lesions are often flat and inconspicuous in color. A prototype ultraviolet spectral imager measuring autofluorescence (AF) and reflectance has been developed and applied in a study of 21 fresh human colon surgical specimens. Six excitation wavelengths from 280 to 440 nm and formulaic ratio imaging were utilized to increase lesion contrast and cause neoplasms to appear bright compared to normal tissue. It was found that in the subset of lesions which were most difficult to visualize in standard color photographs [low contrast lesions, (LCLs)] a ratio image (F340/F440) of AF images excited at 340 and 440 nm produced extraordinary images and was effective in about 70% of these difficult cases. Contrast may be due to increased levels of reduced nicotinamide adenine dinucleotide, increased hemoglobin absorption, and reduced signal from submucosal collagen. A second successful ratio image (R480/R555) combined two reflectance images to produce exceptional images especially in particular LCLs where F340/F440 was ineffective. The newly discovered ratio images can potentially improve detection rate in screening with a novel AF colonoscope.

Ultrasound-guided diffuse optical tomography (DOT) is a promising method for characterizing malignant and benign lesions in the female breast. We introduce a new two-step algorithm for DOT inversion in which the optical parameters are estimated with the global optimization method, genetic algorithm. The estimation result is applied as an initial guess to the conjugate gradient (CG) optimization method to obtain the absorption and scattering distributions simultaneously. Simulations and phantom experiments have shown that the maximum absorption and reduced scattering coefficients are reconstructed with less than 10% and 25% errors, respectively. This is in contrast with the CG method alone, which generates about 20% error for the absorption coefficient and does not accurately recover the scattering distribution. A new measure of scattering contrast has been introduced to characterize benign and malignant breast lesions. The results of 16 clinical cases reconstructed with the two-step method demonstrates that, on average, the absorption coefficient and scattering contrast of malignant lesions are about 1.8 and 3.32 times higher than the benign cases, respectively.

Polarization gating is a popular technique in biomedical optics. It is widely used to inspect the surface of the tissues (under colinear or cocircular detection) or instead to probe the volume (cross-linear detection), without information on the probed depth. Elliptical polarization is introduced to explore the possibility of probing diffuse tissues at selective depths. A thorough Monte Carlo simulation study shows complete correlation between the probed depths and the ellipticity of the polarized light, for a medium with known optical properties. Within a wide range of optical parameters, a linear relation between the backscattered intensity and the depth extension of the probed volume was found whatever the polarization used, but with a controlled extension depending on the ellipticity.

Fluorescence microendoscopy can potentially be a powerful modality in minimally invasive percutaneous intervention for cancer diagnosis because it has an exceptional ability to provide micron-scale resolution images in tissues inaccessible to traditional microscopy. After targeting the tumor with guidance by macroscopic images such as computed tomorgraphy or magnetic resonance imaging, fluorescence microendoscopy can help select the biopsy spots or perform an on-site molecular imaging diagnosis. However, one challenge of this technique for percutaneous lung intervention is that the respiratory and hemokinesis motion often renders instability of the sequential image visualization and results in inaccurate quantitative measurement. Motion correction on such serial microscopy image sequences is, therefore, an important post-processing step. We propose a nonlinear motion compensation algorithm using a cubature Kalman filter (NMC-CKF) to correct these periodic spatial and intensity changes, and validate the algorithm using preclinical imaging experiments. The algorithm integrates a longitudinal nonlinear system model using the CKF in the serial image registration algorithm for robust estimation of the longitudinal movements. Experiments were carried out using simulated and real microendoscopy videos captured from the CellVizio 660 system in rabbit VX2 cancer intervention. The results show that the NMC-CKF algorithm yields more robust and accurate alignment results.

Even though surface morphology is always taken into account when assessing clinically pigmented skin lesions, it is not captured by most modern imaging systems using digital imaging. Our aim is to develop a novel three-dimensional (3D) imaging technique to record detailed information of the surface anatomy of melanocytic lesions that will enable improved classification through digital imaging. The apparatus consists of three high-resolution cameras, a light source, and accompanying software. Volume measurements of specific phantoms using volumetric tubes render slightly lower values than those obtained by our 3D imaging system (mean%±SD, 3.8%±0.98, P<0.05). To examine the reproducibility of the method, sequential imaging of melanocytic lesions is carried out. The mean%±SD differences of area, major axis length, volume, and maximum height are 2.1%±1.1, 0.9%±0.8, 3.8%±2.9, and 2.5%±3.5, respectively. Thirty melanocytic lesions are assessed, including common and dysplastic nevi and melanomas. There is a significant difference between nevi and melanomas in terms of variance in height and boundary asymmetry (P<0.001). Moreover, dysplastic nevi have significantly higher variances in pigment density values than common nevi (P<0.001). Preliminary data suggest that our instrument has great potential in the evaluation of the melanocytic lesions. However, these findings should be confirmed in larger-scale studies.

We describe a novel, noninvasive measurement approach for recording a small involuntary tremor of the eye known as ocular microtremor. The method is based on measuring out-of-plane angular displacements of a target by using laser-speckle correlation of images recorded in the Fourier plane of a lens. The system has a dynamic range of 4 to 5000 μrad, resolution of 4 μrad, and a bandwidth of 250 Hz. The design and optimization of the system is presented with an in vitro validation of the system against its specification.

Angiogenesis is a hallmark of many conditions, including cancer, stroke, vascular disease, diabetes, and high-altitude exposure. We have previously shown that one can study angiogenesis in animal models by using total hemoglobin (tHb) as a marker of cerebral blood volume (CBV), measured using broadband near-infrared spectroscopy (bNIRS). However, the method was not suitable for patients as global anoxia was used for the calibration. Here we determine if angiogenesis could be detected using a calibration method that could be applied to patients. CBV, as a marker of angiogenesis, is quantified in a rat cortex before and after hypoxia acclimation. Rats are acclimated at 370-mmHg pressure for three weeks, while rats in the control group are housed under the same conditions, but under normal pressure. CBV increased in each animal in the acclimation group. The mean CBV (%volume/volume) is 3.49%±0.43% (mean±SD) before acclimation for the experimental group, and 4.76%±0.29% after acclimation. The CBV for the control group is 3.28%±0.75%, and 3.09%±0.48% for the two measurements. This demonstrates that angiogenesis can be monitored noninvasively over time using a bNIRS system with a calibration method that is compatible with human use and less stressful for studies using animals.

We demonstrate a 1050-nm spectral domain optical coherence tomography (OCT) system with a 12 mm imaging depth in air, a 120 kHz A-scan rate and a 10 μm axial resolution for anterior-segment imaging of human eye, in which a new prototype InGaAs linescan camera with 2048 active-pixel photodiodes is employed to record OCT spectral interferograms in parallel. Combined with the full-range complex technique, we show that the system delivers comparable imaging performance to that of a swept-source OCT with similar system specifications.

Ultrathin flexible fiberscopes typically have separate illumination and imaging channels and are available in diameters ranging from 0.5 to 2.5 mm. Diameters can potentially be reduced by combining the illumination and imaging paths into a single fiberoptic channel. Single-channel fiberscopes must incorporate a system to minimize Fresnel reflections from air–glass interfaces within the common illumination and detection path. The Fresnel reflection at the proximal surface of the fiber bundle is particularly problematic. This paper describes and compares methods to reduce the background signal from the proximal surface of the fiber bundle. Three techniques are evaluated: (1) antireflective (AR)-coating the proximal face of the fiber, (2) incorporating crossed polarizers into the light path, and (3) a novel technique called numerical aperture sharing, whereby a portion of the image numerical aperture is devoted to illumination and a portion to detection.

Flow chamber assays, in which blood is perfused over surfaces of immobilized extracellular matrix proteins, are used to investigate the formation of platelet thrombi and aggregates under shear flow conditions. Elucidating the dynamic response of thrombi/aggregate formation to different coagulation pathway perturbations in vitro has been used to develop an understanding of normal and pathological cardiovascular states. Current microscopy techniques, such as differential interference contrast (DIC) or fluorescent confocal imaging, respectively, do not provide a simple, quantitative understanding of the basic physical features (volume, mass, and density) of platelet thrombi/aggregate structures. The use of two label-free imaging techniques applied, for the first time, to platelet aggregate and thrombus formation are introduced: noninterferometric quantitative phase microscopy, to determine mass, and Hilbert transform DIC microscopy, to perform volume measurements. Together these techniques enable a quantitative biophysical characterization of platelet aggregates and thrombi formed on three surfaces: fibrillar collagen, fibrillar collagen +0.1  nM tissue factor (TF), and fibrillar collagen +1  nM TF. It is demonstrated that label-free imaging techniques provide quantitative insight into the mechanisms by which thrombi and aggregates are formed in response to exposure to different combinations of procoagulant agonists under shear flow.

We introduce a novel approach for determination of volume and shape of individual blood platelets modeled as an oblate spheroid from angle-resolved light scattering with flow-cytometric technique. The light-scattering profiles (LSPs) of individual platelets were measured with the scanning flow cytometer and the platelet characteristics were determined from the solution of the inverse light-scattering problem using the precomputed database of theoretical LSPs. We revealed a phenomenon of parameter compensation, which is partly explained in the framework of anomalous diffraction approximation. To overcome this problem, additional a priori information on the platelet refractive index was used. It allowed us to determine the size of each platelet with subdiffraction precision and independent of the particular value of the platelet aspect ratio. The shape (spheroidal aspect ratio) distributions of platelets showed substantial differences between native and activated by 10 μM adenosine diphosphate samples. We expect that the new approach may find use in hematological analyzers for accurate measurement of platelet volume distribution and for determination of the platelet activation efficiency.

We have assessed the value of single fiber reflectance (SFR) spectroscopy in prediction of cervical squamous intraepithelial lesions (SIL). SFR was used to measure reflected light from 32 patients undergoing standard colposcopy. Seven parameters extracted from the spectra in addition to two biographic parameters were compared in biopsy-confirmed SIL versus nonSIL. The significant parameters in the model were determined using stepwise logistic regression. The classification performance was evaluated by a leave-one-out cross-validation method and reported by receiver operating characteristic (ROC) curves. Light absorption properties and biographic characteristics of the patient contributed significantly to the accuracy of the model. Combining important parameters, the best retrospective sensitivity, specificity and area under the ROC curve for SIL sites versus nonSIL were 89%, 80% and 0.89%, respectively. SFR spectroscopy shows promise as a noninvasive, real-time method to guide the clinician in reducing the number of unnecessary biopsies. Discrimination of SIL from other abnormalities compares favorably with that obtained by fluorescence alone and by fluorescence combined with reflectance spectroscopy while the simplicity and low cost of the presented system are dominant.

The reduction of trial-to-trial variability (TTV) in task-evoked functional near-infrared spectroscopy signals by considering the correlated low-frequency spontaneous fluctuations that account for the resting-state functional connectivity in the brain is investigated. A resting-state session followed by a task-state session of a right hand finger-tapping task has been performed on five subjects. Significant ipsilateral and bilateral resting-state functional connectivity has been detected at the subjects’ motor cortex using the seed correlation method. The correlation coefficients obtained during the resting-state are used to reduce the TTV in the signals measured during the task sessions. The results suggest that correlated spontaneous low-frequency fluctuations contribute significantly to the TTV in the task evoked fNIRS signals.

Our current research in the framework of an interdisciplinary project focuses on modelling the dynamics of the hemoglobin reoxygenation process in post-mortem human skin by reflectance spectrometry. The observations of reoxygenation of hemoglobin in livores after postmortem exposure to a cold environment relate the reoxygenation to the commonly known phenomenon that the color impression of livores changes from livid to pink under low ambient temperatures. We analyze the spectra with respect to a physical model describing the optical properties of human skin, discuss the dynamics of the reoxygenation, and propose a phenomenological model for reoxygenation. For additional characterization of the reflectance spectra, the curvature of the local minimum and maximum in the investigated spectral range is considered. There is a strong correlation between the curvature of specra at a wavelength of 560 nm and the concentration of O2-Hb. The analysis is carried out via C programs, as well as MySQL database queries in Java EE, JDBC, Matlab, and Python.

Stokes shift spectroscopy (S3) offers a novel and simpler way to rapidly recognize spectral fingerprints of multiple fluorophores in complex media such as in tissue. This spectroscopic technique can be used as an effective approach to detect cancer in tissue. The alterations of the measured S3 spectra between cancerous and normal tissues were observed in human breast and prostate samples. In order to obtain the optimal Stokes shift interval, Δλ i , for the purpose of breast/prostate cancer detection using S3, the S3 spectra of a mixed aqueous solution of tryptophan, nicotinamide adenine dinucleotide, and flavin were measured with different Δλ i values. The experimental results analyzed using nonnegative least square method show that there is a reduced contribution from collagen and an increased contribution from tryptophan to the S3 signal of the cancerous tissue as compared with those of the normal tissue. This study indicates that the changes of relative contents of tryptophan and collagen in tissue shown by the S3 spectra may present potential native biomarkers for breast and prostate cancer detection. S3 has the potential to be a new armamentarium.

The arterial stiffness evaluation is largely reported as an independent predictor of cardiovascular diseases. The central pulse waveform can provide important data about arterial health and has been studied in patients with several pathologies, such as diabetes mellitus, coronary artery disease and hypertension. The implementation and feasibility studies of a fiber Bragg grating probe for noninvasive monitoring of the carotid pulse are described based on fiber Bragg grating technology. Assessment tests were carried out in carotids of different volunteers and it was possible to detect the carotid pulse waveform in all subjects. In one of the subjects, the sensor was also tested in terms of repeatability. Although further tests will be required for clinical investigation, the first studies suggest that the developed sensor can be a valid alternative to electromechanical tonometers.

Noninvasive glucose monitoring is being investigated as a tool for effectively managing diabetes mellitus. Optical polarimetry has emerged as one such method, which can potentially be used to ascertain blood glucose levels by measuring the aqueous humor glucose levels in the anterior chamber of the eye. The key limitation for realizing this technique is the presence of sample noise due to corneal birefringence, which in the presence of motion artifact can confound the glucose signature in the aqueous humor of the eye. We present the development and characterization of a real-time, closed-loop, dual-wavelength polarimetric system for glucose monitoring using both a custom-built plastic eye phantom (in vitro) and isolated rabbit corneas (ex vivo) mounted in an artificial anterior chamber. The results show that the system can account for these noise sources and can monitor physiologic glucose levels accurately for a limited range of motion-induced birefringence. Using the dual-wavelength system in vitro and ex vivo, standard errors were 14.5  mg/dL and 22.4  mg/dL, respectively, in the presence of birefringence with motion. The results indicate that although dual-wavelength polarimetry has a limited range of compensation for motion-induced birefringence, when aligned correctly, it can minimize the effect of time-varying corneal birefringence for a range of motion larger than what has been reported in vivo.

Modeling behavior of broadband (30 to 1000 MHz) frequency modulated near-infrared (NIR) photons through a phantom is the basis for accurate extraction of optical absorption and scattering parameters of biological turbid media. Photon dynamics in a phantom are predicted using both analytical and numerical simulation and are related to the measured insertion loss (IL) and insertion phase (IP) for a given geometry based on phantom optical parameters. Accuracy of the extracted optical parameters using finite element method (FEM) simulation is compared to baseline analytical calculations from the diffusion equation (DE) for homogenous brain phantoms. NIR spectroscopy is performed using custom-designed, broadband, free-space optical transmitter (Tx) and receiver (Rx) modules that are developed for photon migration at wavelengths of 680, 780, and 820 nm. Differential detection between two optical Rx locations separated by 0.3 cm is employed to eliminate systemic artifacts associated with interfaces of the optical Tx and Rx with the phantoms. Optical parameter extraction is achieved for four solid phantom samples using the least-square-error method in MATLAB (for DE) and COMSOL (for FEM) simulation by fitting data to measured results over broadband and narrowband frequency modulation. Confidence in numerical modeling of the photonic behavior using FEM has been established here by comparing the transmission mode’s experimental results with the predictions made by DE and FEM for known commercial solid brain phantoms.

This article [J. Biomed. Opt.. 17, (8 ), 081417 (2012)] was originally published online on 2 August 2012 with an error in the reference citation of Fig. 1. The correct figure caption and references for Fig. 1 are shown below. A new reference was added to the reference list, shown here.