Fienup Algorithm With Sparsity Constraints: Application to Frequency-Domain Optical-Coherence Tomography

We address the problem of reconstructing a sparse signal from its DFT magnitude. We refer to this problem as the sparse phase retrieval (SPR) problem, which finds applications in tomography, digital holography, electron microscopy, etc. We develop a Fienup-type iterative algorithm, referred to as the Max-K algorithm, to enforce sparsity and successively refine the estimate of phase. We show that the Max-K algorithm possesses Cauchy convergence properties under certain conditions, that is, the MSE of reconstruction does not increase with iterations. We also formulate the problem of SPR as a feasibility problem, where the goal is to find a signal that is sparse in a known basis and whose Fourier transform magnitude is consistent with the measurement. Subsequently, we interpret the Max-K algorithm as alternating projections onto the object-domain and measurement-domain constraint sets and generalize it to a parameterized relaxation, known as the relaxed averaged alternating reflections (RAAR) algorithm. On the application front, we work with measurements acquired using a frequency-domain optical-coherence tomography (FDOCT) experimental setup. Experimental results on measured data show that the proposed algorithms exhibit good reconstruction performance compared with the direct inversion technique, homomorphic technique, and the classical Fienup algorithm without sparsity constraint; specifically, the autocorrelation artifacts and background noise are suppressed to a significant extent. We also demonstrate that the RAAR algorithm offers a broader framework for FDOCT reconstruction, of which the direct inversion technique and the proposed Max-K algorithm become special instances corresponding to specific values of the relaxation parameter.

Super-Resolution Reconstruction in Frequency-Domain Optical-Coherence Tomography Using the Finite-Rate-of-Innovation Principle

The standard approach to signal reconstruction in frequency-domain optical-coherence tomography (FDOCT) is to apply the inverse Fourier transform to the measurements. This technique offers limited resolution (due to Heisenberg's uncertainty principle). We propose a new super-resolution reconstruction method based on a parametric representation. We consider multilayer specimens, wherein each layer has a constant refractive index and show that the backscattered signal from such a specimen fits accurately in to the framework of finite-rate-of-innovation (FRI) signal model and is represented by a finite number of free parameters. We deploy the high-resolution Prony method and show that high-quality, super-resolved reconstruction is possible with fewer measurements (about one-fourth of the number required for the standard Fourier technique). To further improve robustness to noise in practical scenarios, we take advantage of an iterated singular-value decomposition algorithm (Cadzow denoiser). We present results of Monte Carlo analyses, and assess statistical efficiency of the reconstruction techniques by comparing their performance against the Cramer-Rao bound. Reconstruction results on experimental data obtained from technical as well as biological specimens show a distinct improvement in resolution and signal-to-reconstruction noise offered by the proposed method in comparison with the standard approach.

Optimum Parameter Selection in Sparse Reconstruction of Frequency-Domain Optical-Coherence Tomography Signals

For a multilayered specimen, the back-scattered signal in frequency-domain optical-coherence tomography (FDOCT) is expressible as a sum of cosines, each corresponding to a change of refractive index in the specimen. Each of the cosines represent a peak in the reconstructed tomogram. We consider a truncated cosine series representation of the signal, with the constraint that the coefficients in the basis expansion be sparse. An l(2) (sum of squared errors) data error is considered with an l(1) (summation of absolute values) constraint on the coefficients. The optimization problem is solved using Weiszfeld's iteratively reweighted least squares (IRLS) algorithm. On real FDOCT data, improved results are obtained over the standard reconstruction technique with lower levels of background measurement noise and artifacts due to a strong l(1) penalty. The previous sparse tomogram reconstruction techniques in the literature proposed collecting sparse samples, necessitating a change in the data capturing process conventionally used in FDOCT. The IRLS-based method proposed in this paper does not suffer from this drawback.

Dynamic full-range Fourier-domain optical coherence tomography using sinusoidal phase-modulating interferometry

We propose a technique for dynamic full-range Fourier-domain optical coherence tomography by using sinusoidal phase-modulating interferometry, where both the full-range structural information and depth-resolved dynamic information are obtained. A novel frequency-domain filtering algorithm is proposed to reconstruct a time-dependent complex spectral interferogram from the sinusoidally phase-modulated interferogram detected with a high-rate CCD camera. By taking the amplitude and phase of the inverse Fourier transform of the complex spectral interferogram, a time-dependent full-range cross-sectional image and depth-resolved displacement are obtained. Displacement of a sinusoidally vibrating glass cover slip behind a fixed glass cover slip is measured with subwavelength sensitivity to demonstrate the depth-resolved dynamic imaging capability of our system. (c) 2007 Society of Photo-Optical Instrumentation Engineers.

Full-range parallel Fourier-domain optical coherence tomography using sinusoidal phase-modulating interferometry

We demonstrate a full-range parallel Fourier-domain optical coherence tomography (FD-OCT) in which a tomogram free of mirror images as well as DC and autocorrelation terms is obtained in parallel. The phase and amplitude of two-dimensional spectral interferograms are accurately detected by using sinusoidal phase-modulating interferometry and a two-dimensional CCD camera, which allows for the reconstruction of two-dimensional complex spectral interferograms. By line-by-line inverse Fourier transformation of the two-dimensional complex spectral interferogram, a full-range parallel FD-OCT is realized. Tomographic images of two separated glass coverslips obtained with our method are presented as a proof-of-principle experiment.

Dynamic full-range Fourier-domain optical coherence tomography using sinusoidal phase-modulating interferometry

We propose a technique for dynamic full-range Fourier-domain optical coherence tomography by using sinusoidal phase-modulating interferometry, where both the full-range structural information and depth-resolved dynamic information are obtained. A novel frequency-domain filtering algorithm is proposed to reconstruct a time-dependent complex spectral interferogram from the sinusoidally phase-modulated interferogram detected with a high-rate CCD camera. By taking the amplitude and phase of the inverse Fourier transform of the complex spectral interferogram, a time-dependent full-range cross-sectional image and depth-resolved displacement are obtained. Displacement of a sinusoidally vibrating glass cover slip behind a fixed glass cover slip is measured with subwavelength sensitivity to demonstrate the depth-resolved dynamic imaging capability of our system. (c) 2007 Society of Photo-Optical Instrumentation Engineers.

Advanced Monte Carlo Simulation and Machine Learning for Frequency Domain Optical Coherence Tomography

Optical Coherence Tomography(OCT) is a popular, rapidly growing imaging technique with an increasing number of bio-medical applications due to its noninvasive nature. However, there are three major challenges in understanding and improving an OCT system: (1) Obtaining an OCT image is not easy. It either takes a real medical experiment or requires days of computer simulation. Without much data, it is difficult to study the physical processes underlying OCT imaging of different objects simply because there aren't many imaged objects. (2) Interpretation of an OCT image is also hard. This challenge is more profound than it appears. For instance, it would require a trained expert to tell from an OCT image of human skin whether there is a lesion or not. This is expensive in its own right, but even the expert cannot be sure about the exact size of the lesion or the width of the various skin layers. The take-away message is that analyzing an OCT image even from a high level would usually require a trained expert, and pixel-level interpretation is simply unrealistic. The reason is simple: we have OCT images but not their underlying ground-truth structure, so there is nothing to learn from. (3) The imaging depth of OCT is very limited (millimeter or sub-millimeter on human tissues). While OCT utilizes infrared light for illumination to stay noninvasive, the downside of this is that photons at such long wavelengths can only penetrate a limited depth into the tissue before getting back-scattered. To image a particular region of a tissue, photons first need to reach that region. As a result, OCT signals from deeper regions of the tissue are both weak (since few photons reached there) and distorted (due to multiple scatterings of the contributing photons). This fact alone makes OCT images very hard to interpret.

This thesis addresses the above challenges by successfully developing an advanced Monte Carlo simulation platform which is 10000 times faster than the state-of-the-art simulator in the literature, bringing down the simulation time from 360 hours to a single minute. This powerful simulation tool not only enables us to efficiently generate as many OCT images of objects with arbitrary structure and shape as we want on a common desktop computer, but it also provides us the underlying ground-truth of the simulated images at the same time because we dictate them at the beginning of the simulation. This is one of the key contributions of this thesis. What allows us to build such a powerful simulation tool includes a thorough understanding of the signal formation process, clever implementation of the importance sampling/photon splitting procedure, efficient use of a voxel-based mesh system in determining photon-mesh interception, and a parallel computation of different A-scans that consist a full OCT image, among other programming and mathematical tricks, which will be explained in detail later in the thesis.

Next we aim at the inverse problem: given an OCT image, predict/reconstruct its ground-truth structure on a pixel level. By solving this problem we would be able to interpret an OCT image completely and precisely without the help from a trained expert. It turns out that we can do much better. For simple structures we are able to reconstruct the ground-truth of an OCT image more than 98% correctly, and for more complicated structures (e.g., a multi-layered brain structure) we are looking at 93%. We achieved this through extensive uses of Machine Learning. The success of the Monte Carlo simulation already puts us in a great position by providing us with a great deal of data (effectively unlimited), in the form of (image, truth) pairs. Through a transformation of the high-dimensional response variable, we convert the learning task into a multi-output multi-class classification problem and a multi-output regression problem. We then build a hierarchy architecture of machine learning models (committee of experts) and train different parts of the architecture with specifically designed data sets. In prediction, an unseen OCT image first goes through a classification model to determine its structure (e.g., the number and the types of layers present in the image); then the image is handed to a regression model that is trained specifically for that particular structure to predict the length of the different layers and by doing so reconstruct the ground-truth of the image. We also demonstrate that ideas from Deep Learning can be useful to further improve the performance.

It is worth pointing out that solving the inverse problem automatically improves the imaging depth, since previously the lower half of an OCT image (i.e., greater depth) can be hardly seen but now becomes fully resolved. Interestingly, although OCT signals consisting the lower half of the image are weak, messy, and uninterpretable to human eyes, they still carry enough information which when fed into a well-trained machine learning model spits out precisely the true structure of the object being imaged. This is just another case where Artificial Intelligence (AI) outperforms human. To the best knowledge of the author, this thesis is not only a success but also the first attempt to reconstruct an OCT image at a pixel level. To even give a try on this kind of task, it would require fully annotated OCT images and a lot of them (hundreds or even thousands). This is clearly impossible without a powerful simulation tool like the one developed in this thesis.

One-shot parallel complex Fourier-domain optical coherence tomography using a spatial carrier frequency

We propose a novel method of one-shot parallel complex Fourier-domain optical coherence tomography using a spatial carrier frequency for full range imaging. The spatial carrier frequency is introduced into the 2-D spectral interferogram in the lateral direction by using a tilted reference wavefront. This spatial-carrier- contained 2-D spectral interferogram is recorded with one shot of a 2-D CCD camera, and is Fourier-transformed in the lateral direction to obtain a 2-D complex spectral interferogram by a spatial-carrier technique. A full-range tomogram is reconstructed from the 2-D complex spectral interferogram. The principle of this method is confirmed by cross-sectional imaging of a glass slip object. (c) 2008 Society of Photo-Optical Instrumentation Engineers.

Intraoperative spectral domain optical coherence tomography for vitreoretinal surgery.

We demonstrate in vivo human retinal imaging using an intraoperative microscope-mounted optical coherence tomography system (MMOCT). Our optomechanical design adapts an Oculus Binocular Indirect Ophthalmo Microscope (BIOM3), suspended from a Leica ophthalmic surgical microscope, with spectral domain optical coherence tomography (SD-OCT) scanning and relay optics. The MMOCT enables wide-field noncontact real-time cross-sectional imaging of retinal structure, allowing for SD-OCT augmented intrasurgical microscopy for intraocular visualization. We experimentally quantify the axial and lateral resolution of the MMOCT and demonstrate fundus imaging at a 20Hz frame rate.

3D refraction correction and extraction of clinical parameters from spectral domain optical coherence tomography of the cornea.

Capable of three-dimensional imaging of the cornea with micrometer-scale resolution, spectral domain-optical coherence tomography (SDOCT) offers potential advantages over Placido ring and Scheimpflug photography based systems for accurate extraction of quantitative keratometric parameters. In this work, an SDOCT scanning protocol and motion correction algorithm were implemented to minimize the effects of patient motion during data acquisition. Procedures are described for correction of image data artifacts resulting from 3D refraction of SDOCT light in the cornea and from non-idealities of the scanning system geometry performed as a pre-requisite for accurate parameter extraction. Zernike polynomial 3D reconstruction and a recursive half searching algorithm (RHSA) were implemented to extract clinical keratometric parameters including anterior and posterior radii of curvature, central cornea optical power, central corneal thickness, and thickness maps of the cornea. Accuracy and repeatability of the extracted parameters obtained using a commercial 859nm SDOCT retinal imaging system with a corneal adapter were assessed using a rigid gas permeable (RGP) contact lens as a phantom target. Extraction of these parameters was performed in vivo in 3 patients and compared to commercial Placido topography and Scheimpflug photography systems. The repeatability of SDOCT central corneal power measured in vivo was 0.18 Diopters, and the difference observed between the systems averaged 0.1 Diopters between SDOCT and Scheimpflug photography, and 0.6 Diopters between SDOCT and Placido topography.

Optical coherence tomography is a valuable tool in the study of the effects of microneedle geometry on skin penetration characteristics and in-skin dissolution.

In this study, we used optical coherence tomography (OCT) to extensively investigate, for the first time, the effect that microneedle (MN) geometry (MN height, and MN interspacing) and force of application have upon penetration characteristics of soluble poly(methylvinylether-co-maleic anhydride, PMVE/MA) MN arrays into neonatal porcine skin in vitro. The results from OCT investigations were then used to design optimal and suboptimal MN-based drug delivery systems and evaluate their drug delivery profiles cross full thickness and dermatomed neonatal porcine skin in vitro. It was found that increasing the force used for MN application resulted in a significant increase in the depth of penetration achieved within neonatal porcine skin. For example, MN of 600 µm height penetrated to a depth of 330 µm when inserted at a force of 4.4 N/array, while the penetration increased significantly to a depth of 520 µm, when the force of application was increased to 16.4 N/array. At an application force of 11.0 N/array it was found that, in each case, increasing MN height from 350 to 600 µm to 900 µm led to a significant increase in the depth of MN penetration achieved. Moreover, alteration of MN interspacing had no effect upon depth of penetration achieved, at a constant MN height and force of application. With respect to MN dissolution, an approximate 34% reduction in MN height occurred in the first 15 min, with only 17% of the MN height remaining after a 3-hour period. Across both skin models, there was a significantly greater cumulative amount of theophylline delivered after 24 h from an MN array of 900 µm height (292.23 ± 16.77 µg), in comparison to an MN array of 350 µm height (242.62 ± 14.81 µg) (p < 0.001). Employing full thickness skin significantly reduced drug permeation in both cases. Importantly, this study has highlighted the effect that MN geometry and application force have upon the depth of penetration into skin. While it has been shown that MN height has an important role in the extent of drug delivered across neonatal porcine skin from a soluble MN array, further studies to evaluate the full significance of MN geometry on MN mediated drug delivery are now underway. The successful use of OCT in this study could prove to be a key development for polymeric MN research, accelerating their commercial exploitation.

The Quality of Reporting of Diagnostic Accuracy Studies of Optical Coherence Tomography in Glaucoma

Objective: To evaluate the quality of reporting of diagnostic accuracy studies using optical coherence tomography (OCT) in glaucoma. Design: Descriptive series of published studies. Participants: Published studies reporting a measure of the diagnostic accuracy of OCT for glaucoma. Methods: Review of English language papers reporting measures of diagnostic accuracy of OCT for glaucoma. Papers were identified from a Medline literature search performed in June 2006. Articles were appraised using the 25 items provided by the Standards for Reporting of Diagnostic Accuracy (STARD) initiative. Each item was recorded as full, partially, or not reported. Main Outcome Measures: Degree of compliance with the STARD guidelines. Results: Thirty papers were appraised. Eight papers (26.7%) fully reported more than half of the STARD items. The lowest number of fully reported items in a study was 5 and the highest was 17. Descriptions of key aspects of methodology frequently were missing. For example, details of participant sampling (e.g., consecutive or random selection) were described in only 8 (26.7%) of 30 publications. Measures of statistical uncertainty were reported in 18 (60%) of 30 publications. No single STARD item was fully reported by all the papers. Conclusions: The standard of reporting of diagnostic accuracy studies in glaucoma using OCT was suboptimal. It is hoped that adoption of the STARD guidelines will lead to an improvement in reporting of diagnostic accuracy studies, enabling clearer evidence to be produced for the usefulness of OCT for the diagnosis of glaucoma. © 2007 American Academy of Ophthalmology.