Image correction scheme applied to functional diffuse optical tomography scattering images

We have extended our investigation on the use of a linear algorithm for enhancing the accuracy of diffuse optical tomography (DOT) images, to include spatial maps of the diffusion coefficient. The results show that the corrected images are markedly improved in terms of estimated size, spatial resolution, two-object resolving power, and quantitative accuracy. These image-enhancing effects are significant at expected levels of diffusion-coefficient contrast in tissue and noise levels typical of experimental DOT data. Overall, the types and magnitudes of image-enhancing effects obtained here are qualitatively similar to those seen in previous studies on mu(a) perturbations. The implications for practical implementations of DOT time-series imaging are discussed.     

Effects of incident illumination on in-line phase-contrast imaging

Effects of incident illumination on phase-contrast images obtained by means of free-space propagation are investigated under the "transport-of-intensity" approximation. Analytical expressions for image intensity distribution are derived in the cases of coherent quasi-plane and quasi-spherical incident waves, as well as for spatially incoherent and quasi-homogeneous sources and some other types of sources. Practical methods for measuring the relevant parameters of the incident radiation are discussed together with formulas allowing one to calculate the effect of these parameters on the image intensity distribution. The results are expected to be useful in quantitative in-line imaging, phase retrieval, and tomography with polychromatic and spatially partially coherent radiation. As an application we present a method for simultaneous "automatic" phase retrieval and spatial deconvolution in in-line imaging of homogeneous objects using extended polychromatic x-ray sources.     

Improvement of image quality in diffuse optical tomography by use of full time-resolved data

In the field of diffuse optical tomography (DOT), it is widely accepted that time-resolved (TR) measurement can provide the richest information on photon migration in a turbid medium, such as biological tissue. However, the currently available image reconstruction algorithms for TR DOT are based mostly on the cw component or some featured data types of original temporal profiles, which are related to the solution of a time-independent diffusion equation. Although this methodology can greatly simplify the reconstruction process, it suffers from low spatial resolution and poor quantitativeness owing to the limitation of effectively applicable data types. To improve image quality, it has been argued that exploiting the full TR data is essential. We propose implementation of a DOT algorithm by using full TR data and furthermore a variant algorithm with time slices of TR data to alleviate the computational complexity and enhance noise robustness. Compared with those algorithms where the featured data types are used, our evaluations on the spatial resolution and quantitativeness show that a significant improvement in imaging quality can be achieved when full TR data are used, which convinces the DOT community of the potential advantage of the TR domain over cw and frequency domains.     

Optical tomography by the temporally extrapolated absorbance method

The concept of the temporally extrapolated absorbance method (TEAM) for optical tomography of turbid media has been verified by fundamental experiments and image reconstruction. The TEAM uses the time-resolved spectroscopic data of the reference and object to provide projection data that are processed by conventional backprojection. Optical tomography images of a phantom consisting of axisymmetric double cylinders were experimentally obtained with the TEAM and time-gating and continuous-wave (CW) methods. The reconstructed TEAM images are compared with those obtained with the time-gating and CW methods and are found to have better spatial resolution.     

High-resolution full-field optical coherence tomography with a Linnik microscope

We describe an original microscope for high-resolution optical coherence tomography applications. Our system is based on a Linnik interference microscope with high-numerical-aperture objectives. Lock-in detection of the interference signal is achieved in parallel on a CCD by use of a photoelastic birefringence modulator and full-field stroboscopic illumination with an infrared LED. Transverse cross-section (en-face, or XY) images can be obtained in real time with better than 1-microm axial (Z) resolution and 0.5-microm transverse (XY) resolution. A sensitivity of approximately 80 dB is reached at a 1-image/s acquisition rate, which allows tomography in scattering media such as biological tissues.     

Natural light illumination system

In recent years, green energy has undergone a lot of development and has been the subject of many applications. Many research studies have focused on illumination with sunlight as a means of saving energy and creating healthy lighting. Natural light illumination systems have collecting, transmitting, and lighting elements. Today, most daylight collectors use dynamic concentrators; these include Sun tracking systems. However, this design is too expensive to be cost effective. To create a low-cost collector that can be easily installed on a large building, we have designed a static concentrator, which is prismatic and cascadable, to collect sunlight for indoor illumination. The transmission component uses a large number of optical fibers. Because optical fibers are expensive, this means that most of the cost for the system will be related to transmission. In this paper, we also use a prismatic structure to design an optical coupler for coupling n to 1. With the n-to-1 coupler, the number of optical fibers necessary can be greatly reduced. Although this new natural light illumination system can effectively guide collected sunlight and send it to the basement or to other indoor places for healthy lighting, previously there has been no way to manage the collected sunlight when lighting was not desired. To solve this problem, we have designed an optical switch and a beam splitter to control and separate the transmitted light. When replacing traditional sources, the lighting should have similar characteristics, such as intensity distribution and geometric parameters, to those of traditional artificial sources. We have designed, simulated, and optimized an illumination lightpipe with a dot pattern to redistribute the collected sunlight from the natural light illumination system such that it equals the qualities of a traditional lighting system. We also provide an active lighting module that provides lighting from the natural light illumination system or LED auxiliary sources, depending on circumstances. The system is controlled by a light detector. We used optical simulation tools to design and simulate the efficiency of the active module. Finally, we used the natural light illumination system to provide natural illumination for a traffic tunnel. This system will provide a great number of benefits for the people who use it.     

Effect of illumination in an integral imaging system with large depth of focus

We present the characteristics of integral imaging systems with large depth of focus (DOF) by use of two kinds of illumination: plane illumination and diffusing illumination. For each system, we perform ray analysis based on ray optics. To check the visual quality through optical experiments, we use an average image of observed images picked up at various positions within a large DOF. The synthesized elemental images for a three-dimensional (3-D) object with two character patterns were displayed in an optical system and its reconstruction experiments are performed. Experimental results show that use of diffusing illumination can improve visual quality of reconstruction 3-D images in depth-priority integral imaging.     

Three-dimensional optical-transfer-function analysis of fiber-optical two-photon fluorescence microscopy

The three-dimensional optical transfer function is derived for analyzing the imaging performance in fiber-optical two-photon fluorescence microscopy. Two types of fiber-optical geometry are considered: The first involves a single-mode fiber for delivering a laser beam for illumination, and the second is based on the use of a single-mode fiber coupler for both illumination delivery and signal collection. It is found that in the former case the transverse and axial cutoff spatial frequencies of the three-dimensional optical transfer function are the same as those in conventional two-photon fluorescence microscopy without the use of a pinhole.However, the transverse and axial cutoff spatial frequencies in the latter case are 1.7 times as large as those in the former case. Accordingly, this feature leads to an enhanced optical sectioning effect when a fiber coupler is used, which is consistent with our recent experimental observation.     

Geometric superresolution by code division multiplexing

In many highly resolved optical systems the resolution is limited not by the optics but by the CCD's nonzero pixel size. As a result, overall resolution is decreased. Here we propose a novel approach to enhancing resolution beyond the limit set by the CCD's pixels. This method does not involve additional mechanical elements, such as those used for microscans. In this scheme neither the CCD nor additional elements are moved. The geometric superresolving procedure is based on code-division multiplexing, with all its inherent benefits, such as relative noise immunity to single-tone interference. A setup is proposed for coherent and incoherent illumination, with slight modifications for the latter. A theoretical analysis of the setup is presented and compared with empirical results. This scheme is shown to enhance one-dimensional image resolution with the use of only a simple mask that doubles image resolution. This method can easily be expanded to two-dimensional images and to resolution-enhancement factors greater than 2.     

Optical transfer functions for specific-shaped apertures generated by illumination with a rectangular light pipe

We investigated the pupil functions of specific-shaped apertures generated by Lambertian illumination with a rectangular light pipe to derive the corresponding optical transfer functions (OTFs) in aberration-free and defocused optical systems. The semianalytical results indicate that the curves of the OTF of the optical system vary with the form of the shaped apertures that are generated by illumination with different geometric structures of light pipes and light sources. It was found that the OTF values of even-peak frequencies decrease when the Lambertian light source decreases. If there are a total of n x n individual apertures within a pupil, then n near-periodical peaks will appear on the OTF curve. It is evident that the values of the OTF remain almost unchanged even when the lengths of the light pipes are different. Furthermore, the geometric structure of the light pipe does not affect the resolution limit of the optical system, and under the condition of a larger defocused coefficient omega(20), the results of the defocused system can coincide with those of the aberration-free systems.     

Fundamental characteristics of a synthesized light source for optical coherence tomography

We describe the fundamental characteristics of a synthesized light source (SLS) consisting of two low-coherence light sources to enhance the spatial resolution for optical coherence tomography (OCT). The axial resolution of OCT is given by half the coherence length of the light source. We fabricated a SLS with a coherence length of 2.3 microm and a side-lobe intensity of 45% with an intensity ratio of LED1:LED2 = 1:0.5 by combining two light sources, LED1, with a central wavelength of 691 nm and a spectral bandwidth of 99 nm, and LED2, with a central wavelength of 882 nm and a spectral bandwidth of 76 nm. The coherence length of 2.3 microm was 56% of the shorter coherence length in the two LEDs, which indicates that the axial resolution is 1.2 microm. The lateral resolution was measured at less than 4.4 microm by use of the phase-shift method and with a test pattern as a sample. The measured rough surfaces of a coin are illustrated and discussed.     

Detection of ultrawide-band ultrasound pulses in optoacoustic tomography

A laser optoacoustic imaging system (LOIS) uses time-resolved detection of laser-induced pressure profiles in tissue in order to reconstruct images of the tissue based on distribution of acoustic sources. Laser illumination with short pulses generates distribution of acoustic sources that accurately replicates the distribution of absorbed optical energy. The complex spatial profile of heterogeneous distribution of acoustic sources can be represented in the frequency domain by a wide spectrum of ultrasound ranging from tens of kilohertz to tens of megahertz. Therefore, LOIS requires a unique acoustic detector operating simultaneously within a wide range of ultrasonic frequencies. Physical principles of an array of ultrawide-band ultrasonic transducers used in LOIS designed for imaging tumors in the depth of tissue are described. The performance characteristics of the transducer array were modeled and compared with experiments performed in gel phantoms resembling optical and acoustic properties of human tissue with small tumors. The amplitude and the spectrum of laser-induced ultrasound pulses were measured in order to determine the transducer sensitivity and the level of thermal noises within the entire ultrasonic band of detection. Spatial resolution of optoacoustic images obtained with an array of piezoelectric transducers and its transient directivity pattern within the field of view are described. The detector design considerations essential for obtaining high-quality optoacoustic images are presented.     

Full field of view super-resolution imaging based on two static gratings and white light illumination

The usage of two static gratings for obtaining super-resolved imaging dates back to the work by Bachl and Lukosz in 1967. However, in their approach a severe reduction in the field of view was the necessary condition for improving the resolution. We present an approach based on two static gratings without sacrificing the field of view. The key idea for not paying with the field of view is to use white light illumination to average the ghost images obtained outside the region of interest since the positions of those images are wavelength dependent. Moreover, large magnification is achieved by using a commercial microscope objective instead of a test system with a unity magnification as presented in previous works. Because of the large magnification, the second grating has a low spatial period. This allows us to create an intermediate image and use a second imaging lens with low resolution capability while still obtaining an imaging quality as good as that provided by the first imaging lens. This is an important improvement in comparison with the original super-resolving method with two fixed gratings.     

复谱频域OCT系统高分辨率图像重建

为了提高复杂多层样品层析图像的分辨率,构建了复谱频域光学相干层析成像(CSOCT)系统.由于其在低亮度和高速成像方面相对于时域OCT具有更高的灵敏度,因此在光学相干层析成像系统中具有重要的作用.本文对测试样品二层盖波片进行成像实验,基于光学相干层析基本理论,采用五帧相移算法,最终获得测试样品的复谱频域OCT图像.实验结果表明,该系统可以消除谱频域OCT图像中的寄生项和镜像,改善和提高层析图像的分辨率.     

Full-field birefringence imaging by thermal-light polarization-sensitive optical coherence tomography. I. Theory

A method for measuring birefringence by use of thermal-light polarization-sensitive optical coherence tomography is presented. The use of thermal light brings to polarization-sensitive optical coherence tomography a resolution in the micrometer range in three dimensions. The instrument is based on a Linnik interference microscope and makes use of achromatic quarter-wave plates. A mathematical representation of the instrument is presented here, and the detection scheme is described, together with a discussion of the validity domain of the equations used to evaluate the birefringence in the presence of white-light illumination.     

Frequency-domain sensitivity analysis for small imaging domains using the equation of radiative transfer

Optical tomography of small imaging domains holds great promise as the signal-to-noise ratio is usually high, and the achievable spatial resolution is much better than in large imaging domains. Emerging applications range from the imaging of joint diseases in human fingers to monitoring tumor growth or brain activity in small animals. In these cases, the diameter of the tissue under investigation is typically smaller than 3 cm, and the optical path length is only a few scattering mean-free paths. It is well known that under these conditions the widely applied diffusion approximation to the equation of radiative transfer (ERT) is of limited applicability. To accurately model light propagation in these small domains, the ERT has to be solved directly. We use the frequency-domain ERT to perform a sensitivity study for small imaging domains. We found optimal source-modulation frequencies for which variations in optical properties, size, and location of a tissue inhomogeneity lead to maximal changes in the amplitude and phase of the measured signal. These results will be useful in the design of experiments and optical tomographic imaging systems that probe small tissue volumes.     

Deconvolution-based spatial resolution in optical diffusion tomography

The role that deconvolution plays in the achievable spatial resolution in optical diffusion tomography is examined for the case of imaging an inhomogeneity in an otherwise homogeneous medium. It is shown that, in the measured data, it is the shape of the turbid medium modulation transfer function that determines the maximum spatial resolution. When the turbid medium transfer function is deconvolved from the measured data, it is the signal-to-noise ratio properties of the Fourier-transformed measured data that determine the maximum spatial resolution. It is shown that deconvolution-based methods can improve the spatial resolution in measured data up to a factor of 5 for realistic noise levels. These results are demonstrated with computer-simulated data.     

Wavelength-multiplexing system for single-mode image transmission

The expanding use of optical communication by means of optical fibers and the situation of drastically increasing amounts of data to be transmitted urge the exploration of novel systems permitting the transmission of large amounts of spatial information by fiber with smaller spatial resolution. An optical encoding and decoding system is suggested for transmitting one- or two-dimensional images by means of a single-mode fiber. The superresolving system is based on wavelength multiplexing of the input spatial information, which is achieved with diffractive optical elements. Preliminary experimental results demonstrate the capabilities of the suggested method for the one- and two-dimensional cases.     

Tomographic image reconstruction from optical projections in light-diffusing media

The recent developments in light generation and detection techniques have opened new possibilities for optical medical imaging, tomography, and diagnosis at tissue penetration depths of ~10 cm. However, because light scattering and diffusion in biological tissue are rather strong, the reconstruction of object images from optical projections needs special attention. We describe a simple reconstruction method for diffuse optical imaging, based on a modified backprojection approach for medical tomography. Specifically, we have modified the standard backprojection method commonly used in x-ray tomographic imaging to include the effects of both the diffusion and the scattering of light and the associated nonlinearities in projection image formation. These modifications are based primarily on the deconvolution of the broadened image by a spatially variant point-spread function that is dependent on the scattering of light in tissue. The spatial dependence of the deconvolution and nonlinearity corrections for the curved propagating ray paths in heterogeneous tissue are handled semiempirically by coordinate transformations. We have applied this method to both theoretical and experimental projections taken by parallel- and fan-beam tomography geometries. The experimental objects were biomedical phantoms with multiple objects, including in vitro animal tissue. The overall results presented demonstrate that image-resolution improvements by nearly an order of magnitude can be obtained. We believe that the tomographic method presented here can provide a basis for rapid, real-time medical monitoring by the use of optical projections. It is expected that such optical tomography techniques can be combined with the optical tissue diagnosis methods based on spectroscopic molecular signatures to result in a versatile optical diagnosis and imaging technology.     

Information loss and reconstruction in diffuse fluorescence tomography

This paper is a theoretical exploration of spatial resolution in diffuse fluorescence tomography. It is demonstrated that, given a fixed imaging geometry, one cannot-relative to standard techniques such as Tikhonov regularization and truncated singular value decomposition-improve the spatial resolution of the optical reconstructions via increasing the node density of the mesh considered for modeling light transport. Using techniques from linear algebra, it is shown that, as one increases the number of nodes beyond the number of measurements, information is lost by the forward model. It is demonstrated that this information cannot be recovered using various common reconstruction techniques. Evidence is provided showing that this phenomenon is related to the smoothing properties of the elliptic forward model that is used in the diffusion approximation to light transport in tissue. This argues for reconstruction techniques that are sensitive to boundaries, such as L1-reconstruction and the use of priors, as well as the natural approach of building a measurement geometry that reflects the desired image resolution.