A computationally efficient approach that computes the optimal regularization parameter for the Tikhonov-minimization scheme is developed for photoacoustic imaging. This approach is based on the least squares-QR decomposition which is a well-known dimensionality reduction technique for a large system of equations. It is shown that the proposed framework is effective in terms of quantitative and qualitative reconstructions of initial pressure distribution enabled via finding an optimal regularization parameter. The computational efficiency and performance of the proposed method are shown using a test case of numerical blood vessel phantom, where the initial pressure is exactly known for quantitative comparison.

A combined two-photon microscopy (TPM) and angiographic optical coherence tomography (OCT) is developed, which can provide molecular, cellular, structural, and vascular information of tissue specimens in vivo. This combined system is implemented by adding an OCT vasculature visualization method to the previous combined TPM and OCT, and then is applied to in vivo tissue imaging. Two animal models, a mouse brain cranial window model and a mouse ear cancer model, are used. Both molecular, cellular information at local regions of tissues, and structural, vascular information at relatively larger regions are visualized in the same sections. In vivo tissue microenvironments are better elucidated by the combined TPM and angiographic OCT.

A new approach to retrieve the attenuation-corrected fluorescence using spatial frequency-domain imaging is demonstrated. Both in vitro and ex vivo experiments showed the technique can compensate for the fluorescence attenuation from tissue absorption and scattering. This approach has potential in molecular image-guided surgery.

Recently, we tested the vision of spectacle corrected subjects using a special hologram (http://dx.doi.org/10.1364/BOE.3.001173). A subject viewing through this hologram sees the images of various numbers at different distances from the eye. Each of these images subtends an angle of 50 ′ at the eye. Each image corresponds to a different amount of divergence or convergence at the eye. The limit of convergence (positive blur) with which a subject can recognize the numbers was measured. Hyperopes could recognize numbers with 0.9 D more of positive blur in comparison with myopes. This difference was not obtained with a standard test chart under white light illumination using positive lenses to provide the blur. The observed difference with the hologram is then attributed to the multivergence nature of the target in the hologram and/or the monochromaticity of the illumination used. Using a hologram of a logMAR chart at a single distance of infinity we have now found that the observed difference is induced by the multivergence target and is due to an effect that is similar to the Mandelbaum effect. We present the details of this experiment and compare the results obtained with our earlier experiment.

We present simulations and measurements with an optimized goniometer for determination of the scattering phase function of suspended particles. We applied the Monte Carlo method, using a radially layered cylindrical geometry and mismatched boundary conditions, in order to investigate the influence of reflections caused by the interfaces of the glass cuvette and the scatterer concentration on the accurate determination of the scattering phase function. Based on these simulations we built an apparatus which allows direct measurement of the phase function from ϑ=7  deg to ϑ=172  deg without any need for correction algorithms. Goniometric measurements on polystyrene and SiO₂ spheres proved this concept. Using the validated goniometer, we measured the phase function of yeast cells, demonstrating the improvement of the new system compared to standard goniometers. Furthermore, the scattering phase function of different fat emulsions, like Intralipid, was determined precisely.

Photoacoustic tomography (PAT) offers structural and functional imaging of living biological tissue with highly sensitive optical absorption contrast and excellent spatial resolution comparable to medical ultrasound (US) imaging. We report the development of a fully integrated PAT and US dual-modality imaging system, which performs signal scanning, image reconstruction, and display for both photoacoustic (PA) and US imaging all in a truly real-time manner. The back-projection (BP) algorithm for PA image reconstruction is optimized to reduce the computational cost and facilitate parallel computation on a state of the art graphics processing unit (GPU) card. For the first time, PAT and US imaging of the same object can be conducted simultaneously and continuously, at a real-time frame rate, presently limited by the laser repetition rate of 10 Hz. Noninvasive PAT and US imaging of human peripheral joints in vivo were achieved, demonstrating the satisfactory image quality realized with this system. Another experiment, simultaneous PAT and US imaging of contrast agent flowing through an artificial vessel, was conducted to verify the performance of this system for imaging fast biological events. The GPU-based image reconstruction software code for this dual-modality system is open source and available for download from http://sourceforge.net/projects/patrealtime.

A high-speed (100 kHz A-scans/s ) complex conjugate resolved 1 μm swept source optical coherence tomography (SS-OCT) system using coherence revival of the light source is suitable for dense three-dimensional (3-D) imaging of the anterior segment. The short acquisition time helps to minimize the influence of motion artifacts. The extended depth range of the SS-OCT system allows topographic analysis of clinically relevant images of the entire depth of the anterior segment of the eye. Patients with the type 1 Boston Keratoprosthesis (KPro) require evaluation of the full anterior segment depth. Current commercially available OCT systems are not suitable for this application due to limited acquisition speed, resolution, and axial imaging range. Moreover, most commonly used research grade and some clinical OCT systems implement a commercially available SS (Axsun) that offers only 3.7 mm imaging range (in air) in its standard configuration. We describe implementation of a common swept laser with built-in k-clock to allow phase stable imaging in both low range and high range, 3.7 and 11.5 mm in air, respectively, without the need to build an external MZI k-clock. As a result, 3-D morphology of the KPro position with respect to the surrounding tissue could be investigated in vivo both at high resolution and with large depth range to achieve noninvasive and precise evaluation of success of the surgical procedure.

Collagen fibers can be imaged with second harmonic generation (SHG) and are associated with efficient tumor cell locomotion. Preferential locomotion along these fibers correlates with a more aggressively metastatic phenotype, and changes in SHG emission properties accompany changes in metastatic outcome. We therefore attempted to elucidate the cellular and molecular machinery that influences SHG in order to understand how the microstructure of tumor collagen fibers is regulated. By quantifying SHG and immunofluorescence (IF) from tumors grown in mice with and without stromal tumor necrosis factor (TNF)-α and in the presence or absence of tumor-associated macrophages (TAMs), we determined that depletion of TAMs alters tumor collagen fibrillar microstructure as quantified by SHG and IF. Furthermore, we determined that abrogation of TNF-α expression by tumor stromal cells also alters fibrillar microstructure and that subsequent depletion of TAMs has no further effect. In each case, metastatic burden correlated with optical readouts of collagen microstructure. Our results implicate TAMs and stromal TNF-α as regulators of breast tumor collagen microstructure and suggest that this regulation plays a role in tumor metastasis. Furthermore, these results indicate that quantification of SHG represents a useful strategy for evaluating the cells and molecular pathways responsible for manipulating fibrillar collagen in breast tumor models.

Lymphatic vessels are a part of the circulatory system that collect plasma and other substances that have leaked from the capillaries into interstitial fluid (lymph) and transport lymph back to the circulatory system. Since lymph is transparent, lymphatic vessels appear as dark hallow vessel-like regions in optical coherence tomography (OCT) cross sectional images. We propose an automatic method to segment lymphatic vessel lumen from OCT structural cross sections using eigenvalues of Hessian filters. Compared to the existing method based on intensity threshold, Hessian filters are more selective on vessel shape and less sensitive to intensity variations and noise. Using this segmentation technique along with optical micro-angiography allows label-free noninvasive simultaneous visualization of blood and lymphatic vessels in vivo. Lymphatic vessels play an important role in cancer, immune system response, inflammatory disease, wound healing and tissue regeneration. Development of imaging techniques and visualization tools for lymphatic vessels is valuable in understanding the mechanisms and studying therapeutic methods in related disease and tissue response.

A microelectromechanical system (MEMS) mirror based endoscopic swept-source optical coherence tomography (SS-OCT) system that can perform three-dimensional (3-D) imaging at high speed is reported. The key component enabling 3-D endoscopic imaging is a two-axis MEMS scanning mirror which has a 0.8×0.8  mm ² mirror plate and a 1.6×1.4  mm ² device footprint. The diameter of the endoscopic probe is only 3.5 mm. The imaging rate of the SS-OCT system is 50  frames/s . OCT images of both human suspicious oral leukoplakia tissue and normal buccal mucosa were taken in vivo and compared. The OCT imaging result agrees well with the histopathological analysis.

A method, is presented for blind unmixing spectrally resolved fluorescence lifetime images. The method is based on the combined analysis of spectral and lifetime phasors and allows unmixing of up to three components without any prior knowledge. Fractional intensities, spectra and decay curves of the individual components can be extracted with this new technique. The reliability and sensitivity are investigated and the possibility of extending the method to unmix more components is discussed. The method is evaluated on mixtures of fluorescent dyes and labeled cells.

Multimodal approaches that combine near-infrared (NIR) and conventional imaging modalities have been shown to improve optical parameter estimation dramatically and thus represent a prevailing trend in NIR imaging. These approaches typically involve applying anatomical templates from magnetic resonance imaging/computed tomography/ultrasound images to guide the recovery of optical parameters. However, merging these data sets using current technology requires multiple software packages, substantial expertise, significant time-commitment, and often results in unacceptably poor mesh quality for optical image reconstruction, a reality that represents a significant roadblock for translational research of multimodal NIR imaging. This work addresses these challenges directly by introducing automated digital imaging and communications in medicine image stack segmentation and a new one-click three-dimensional mesh generator optimized for multimodal NIR imaging, and combining these capabilities into a single software package (available for free download) with a streamlined workflow. Image processing time and mesh quality benchmarks were examined for four common multimodal NIR use-cases (breast, brain, pancreas, and small animal) and were compared to a commercial image processing package. Applying these tools resulted in a fivefold decrease in image processing time and 62% improvement in minimum mesh quality, in the absence of extra mesh postprocessing. These capabilities represent a significant step toward enabling translational multimodal NIR research for both expert and nonexpert users in an open-source platform.

Current fluorescence diffuse optical tomography (fDOT) systems can provide large data sets and, in addition, the unknown parameters to be estimated are so numerous that the sensitivity matrix is too large to store. Alternatively, iterative methods can be used, but they can be extremely slow at converging when dealing with large matrices. A few approaches suitable for the reconstruction of images from very large data sets have been developed. However, they either require explicit construction of the sensitivity matrix, suffer from slow computation times, or can only be applied to restricted geometries. We introduce a method for fast reconstruction in fDOT with large data and solution spaces, which preserves the resolution of the forward operator whilst compressing its representation. The method does not require construction of the full matrix, and thus allows storage and direct inversion of the explicitly constructed compressed system matrix. The method is tested using simulated and experimental data. Results show that the fDOT image reconstruction problem can be effectively compressed without significant loss of information and with the added advantage of reducing image noise.

We have developed an imaging system to extract high contrast images from different layers of biological organisms. Utilizing a digital holographic approach, the system works without scanning through layers of the specimen. In dark-field illumination, scattered light has the main contribution in image formation, but in the case of coherent illumination, this creates a strong speckle noise that reduces the image quality. To remove this restriction, the specimen has been illuminated with various speckle-fields and a hologram has been recorded for each speckle-field. Each hologram has been analyzed separately and the corresponding intensity image has been reconstructed. The final image has been derived by averaging over the reconstructed images. A correlation approach has been utilized to determine the number of speckle-fields required to achieve a desired contrast and image quality. The reconstructed intensity images in different object layers are shown for different sea urchin larvae. Two multimedia files are attached to illustrate the process of digital focusing.

The hemodynamic response to neuronal activation is a well-studied phenomenon in the brain, due to the prevalence of functional magnetic resonance imaging. The retina represents an optically accessible platform for studying lamina-specific neurovascular coupling in the central nervous system; however, due to methodological limitations, this has been challenging to date. We demonstrate techniques for the imaging of visual stimulus-evoked hyperemia in the rat inner retina using Doppler optical coherence tomography (OCT) and OCT angiography. Volumetric imaging with three-dimensional motion correction, en face flow calculation, and normalization of dynamic signal to static signal are techniques that reduce spurious changes caused by motion. We anticipate that OCT imaging of retinal functional hyperemia may yield viable biomarkers in diseases, such as diabetic retinopathy, where the neurovascular unit may be impaired.

We present an endoscopic probe for optical coherence tomography (OCT) equipped with a miniaturized hollow ultrasonic motor that rotates the objective lens and provides an internal channel for the fiber to pass through, enabling 360 deg unobstructed circumferential scanning. This probe has an outer diameter of 1.5 mm, which is ultra-small for motorized probes with an unobstructed view in distal scanning endoscopic OCT. Instead of a mirror or prism, a customized aspheric right-angle lens is utilized, leading to an enlargement of the numerical aperture and thus high transverse resolution. Spectral-domain OCT imaging of bio-tissue and a phantom are demonstrated with resolution of 7.5  μm(axial)×6.6  μm(lateral) and sensitivity of 96 dB.

The number of molecular species suitable for multispectral fluorescence imaging is limited due to the overlap of the emission spectra of indicator fluorophores, e.g., dyes and nanoparticles. To remove fluorophore emission cross-talk in wide-field multispectral fluorescence molecular imaging, we evaluate three different solutions: (1) image stitching, (2) concurrent imaging with cross-talk ratio subtraction algorithm, and (3) frame-sequential imaging. A phantom with fluorophore emission cross-talk is fabricated, and a 1.2-mm ultrathin scanning fiber endoscope (SFE) is used to test and compare these approaches. Results show that fluorophore emission cross-talk could be successfully avoided or significantly reduced. Near term, the concurrent imaging method of wide-field multispectral fluorescence SFE is viable for early stage cancer detection and localization in vivo. Furthermore, a means to enhance exogenous fluorescence target-to-background ratio by the reduction of tissue autofluorescence background is demonstrated.

Purpose: We evaluated the feasibility of Raman spectroscopy (RS) in diagnosis and prognosis of castration-resistant prostate cancer (CRPC) in patients with prostate cancer (PC). Materials and methods: Raman spectra are detected from PC cell lines (LNCaP and C4-2) and tissues using a Labram HR 800 RS. Then, principal component analysis (PCA) and support vector machine (SVM) are applied for prediction. A leave-one-out cross-validation is used to train and test the SVM. Results: There are 50 qualified patients, including 33 with androgen-dependent prostate cancer (ADPC) and 17 with CRPC. The spectral changes at 1126, 1170, 1315 to 1338, and 1447  cm ⁻¹ between CRPC and ADPC are detected in both cells and tissues models, which are assigned to specific amino acids and DNA. PCA/SVM algorithm provided a sensitivity of 88.2% and a specificity of 87.9% for diagnosing CRPC tissues. Furthermore, 14 patients with ADPC progressed to CRPC within 12 months. These patients are separated into two groups depending on whether their cancers progressed to CRPC within 12 months. PCA/SVM could differentiate these two groups with a sensitivity of 85.7% and a specificity of 88.9%. Conclusions: RS has the potential in diagnosis and prognosis of CRPC in clinical practice.

The spectral changes of native fluorophores among normal fibroblasts and cancer cell lines of different metastatic ability are investigated by fluorescence spectroscopy. The normal (fibroblast), moderately metastatic (DU-145), and advanced metastatic (PC-3) cell lines were each selectively excited at 300 nm, and their fluorescence emission spectra are analyzed using principal component analysis to explore the differences of the relative contents of tryptophan and reduced nicotinamide adenine dinucleotide in these cell lines. The results show that the tryptophan emission featured predominantly in the fluorescence spectra of the advanced metastatic cancer cells in comparison with the moderately metastatic cancer and normal cells.

A first accurate measurement of the complex refractive index in an intralipid emulsion is demonstrated, and thereby the average scatterer particle size using standard Mie scattering calculations is extracted. Our method is based on measurement and modeling of the reflectance of a divergent laser beam from the sample surface. In the absence of any definitive reference data for the complex refractive index or particle size in highly turbid intralipid emulsions, we base our claim of accuracy on the fact that our work offers several critically important advantages over previously reported attempts. First, our measurements are in situ in the sense that they do not require any sample dilution, thus eliminating dilution errors. Second, our theoretical model does not employ any fitting parameters other than the two quantities we seek to determine, i.e., the real and imaginary parts of the refractive index, thus eliminating ambiguities arising from multiple extraneous fitting parameters. Third, we fit the entire reflectance-versus-incident-angle data curve instead of focusing on only the critical angle region, which is just a small subset of the data. Finally, despite our use of highly scattering opaque samples, our experiment uniquely satisfies a key assumption behind the Mie scattering formalism, namely, no multiple scattering occurs. Further proof of our method’s validity is given by the fact that our measured particle size finds good agreement with the value obtained by dynamic light scattering.

Patients with diabetes mellitus and hypertension (HT) diseases are predisposed to kidney diseases. The objective of this study was to identify potential biomarkers in the urine of diabetic and hypertensive patients through Raman spectroscopy in order to predict the evolution to complications and kidney failure. Urine samples were collected from control subjects (CTR) and patients with diabetes and HT with no complications (lower risk, LR), high degree of complications (higher risk, HR), and doing blood dialysis (DI). Urine samples were stored frozen (−20°C) before spectral analysis. Raman spectra were obtained using a dispersive spectrometer (830-nm, 300-mW power, and 20-s accumulation). Spectra were then submitted to principal component analysis (PCA) followed by discriminant analysis. The first PCA loading vectors revealed spectral features of urea, creatinine, and glucose. It has been found that the amounts of urea and creatinine decreased as disease evoluted from CTR to LR/HR and DI (PC1, p<0.05), and the amount of glucose increased in the urine of LR/HR compared to CTR (PC3, p<0.05). The discriminating model showed better overall classification rate of 70%. These results could lead to diagnostic information of possible complications and a better disease prognosis.

In abdominal trauma patients, monitoring intestinal perfusion and oxygen consumption is essential during the resuscitation period. Photoplethysmography is an optical technique potentially capable of monitoring these changes in real time to provide the medical staff with a timely and quantitative measure of the adequacy of resuscitation. The challenges for using optical techniques in monitoring hemodynamics in intestinal tissue are discussed, and the solutions to these challenges are presented using a combination of Monte Carlo modeling and theoretical analysis of light propagation in tissue. In particular, it is shown that by using visible wavelengths (i.e., 470 and 525 nm), the perfusion signal is enhanced and the background contribution is decreased compared with using traditional near-infrared wavelengths leading to an order of magnitude enhancement in the signal-to-background ratio. It was further shown that, using the visible wavelengths, similar sensitivity to oxygenation changes could be obtained (over 50% compared with that of near-infrared wavelengths). This is mainly due to the increased contrast between tissue and blood in that spectral region and the confinement of the photons to the thickness of the small intestine. Moreover, the modeling results show that the source to detector separation should be limited to roughly 6 mm while using traditional near-infrared light, with a few centimeters source to detector separation leads to poor signal-to-background ratio. Finally, a visible wavelength system is tested in an in vivo porcine study, and the possibility of monitoring intestinal perfusion changes is showed.

The aim of this study is to compare a two-wavelength light emitting diode–based tissue oximeter (INVOS), which is designed to show trends in tissue oxygenation, with a four-wavelength laser–based oximeter (FORE-SIGHT), designed to deliver absolute values of tissue oxygenation. Simultaneous values of cerebral tissue oxygenation (StO₂) are measured using both devices in 15 term and 15 preterm clinically stable newborns on the first and third day of life. Values are recorded simultaneously in two periods between which oximeter sensor positions are switched to the contralateral side. Agreement between StO₂ values before and after the change of sensor position is analyzed. We find that mean cerebral StO₂ values are similar between devices for term and preterm babies, but INVOS shows StO₂ values spread over a wider range, with wider standard deviations than shown by the FORE-SIGHT. There is relatively good agreement with a bias up to 3.5% and limits of agreement up to 11.8%. Measurements from each side of the forehead show better repeatability for the FORE-SIGHT monitor. We conclude that performance of the two devices is probably acceptable for clinical purposes. Both performed sufficiently well, but the use of FORE-SIGHT may be associated with tighter range and better repeatability of data.

Diffuse reflectance spectra are used to determine the optical properties of biological samples. In medicine and forensic science, the turbid objects under study often possess large absorption and/or scattering properties. However, data analysis is frequently based on the diffusion approximation to the radiative transfer equation, implying that it is limited to tissues where the reduced scattering coefficient dominates over the absorption coefficient. Nevertheless, up to absorption coefficients of 20  mm ⁻¹ at reduced scattering coefficients of 1 and 11.5  mm ⁻¹ , we observed excellent agreement (r ² =0.994 ) between reflectance measurements of phantoms and the diffuse reflectance equation proposed by Zonios et al. [Appl. Opt.38, 6628–6637 (1999)], derived as an approximation to one of the diffusion dipole equations of Farrell et al. [Med. Phys.19, 879–888 (1992)]. However, two parameters were fitted to all phantom experiments, including strongly absorbing samples, implying that the reflectance equation differs from diffusion theory. Yet, the exact diffusion dipole approximation at high reduced scattering and absorption also showed agreement with the phantom measurements. The mathematical structure of the diffuse reflectance relation used, derived by Zonios et al. [Appl. Opt.38, 6628–6637 (1999)], explains this observation. In conclusion, diffuse reflectance relations derived as an approximation to the diffusion dipole theory of Farrell et al. can analyze reflectance ratios accurately, even for much larger absorption than reduced scattering coefficients. This allows calibration of fiber-probe set-ups so that the object’s diffuse reflectance can be related to its absorption even when large. These findings will greatly expand the application of diffuse reflection spectroscopy. In medicine, it may allow the use of blue/green wavelengths and measurements on whole blood, and in forensic science, it may allow inclusion of objects such as blood stains and cloth at crime scenes.

This article [J. Biomed. Opt.. 17, , 076028 (2012)] was originally published online on 27 July 2012 with Figures 1 and 3 reversed. The corrected figures are reprinted below.

This article [J. Biomed. Opt.. 18, , 086002 (2013)] was originally published online on 2 August 2013 with an error in the author list. The first author’s name was reversed, and should be “Raju Poddar” as it appears above. This article was corrected online on 27 August 2013. It appears correctly in print.

Research on single-molecule cellular biophysics is an emerging field that is rapidly advancing based on the development of new instrumentation and experimental designs. To understand complex biological systems, researchers require instrumentation and techniques that span many orders of scale. The measured units include length, force, energy, time, and concentration and the validation of measurements for both in vivo and in vitro studies, as well as the limitations of the techniques and an understanding of the perturbations that the techniques induce in the experimental systems.