skip to main content

Title: Computational 3D microscopy with optical coherence refraction tomography

Optical coherence tomography (OCT) has seen widespread success as anin vivoclinical diagnostic 3D imaging modality, impacting areas including ophthalmology, cardiology, and gastroenterology. Despite its many advantages, such as high sensitivity, speed, and depth penetration, OCT suffers from several shortcomings that ultimately limit its utility as a 3D microscopy tool, such as its pervasive coherent speckle noise and poor lateral resolution required to maintain millimeter-scale imaging depths. Here, we present 3D optical coherence refraction tomography (OCRT), a computational extension of OCT that synthesizes an incoherent contrast mechanism by combining multiple OCT volumes, acquired across two rotation axes, to form a resolution-enhanced, speckle-reduced, refraction-corrected 3D reconstruction. Our label-free computational 3D microscope features a novel optical design incorporating a parabolic mirror to enable the capture of 5D plenoptic datasets, consisting of millimetric 3D fields of view over up to±<#comment/>75∘<#comment/>without moving the sample. We demonstrate that 3D OCRT reveals 3D features unobserved by conventional OCT in fruit fly, zebrafish, and mouse samples.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Optical Society of America
Date Published:
Journal Name:
Page Range / eLocation ID:
Article No. 593
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Illuminating or imaging samples from a broad angular range is essential in a wide variety of computational 3D imaging and resolution-enhancement techniques, such as optical projection tomography, optical diffraction tomography, synthetic aperture microscopy, Fourier ptychographic microscopy, structured illumination microscopy, photogrammetry, and optical coherence refraction tomography. The wider the angular coverage, the better the resolution enhancement or 3D-resolving capabilities. However, achieving such angular ranges is a practical challenge, especially when approaching±<#comment/>90∘<#comment/>or beyond. Often, researchers resort to expensive, proprietary high numerical aperture (NA) objectives or to rotating the sample or source-detector pair, which sacrifices temporal resolution or perturbs the sample. Here, we propose several new strategies for multiangle imaging approaching 4pi steradians using concave parabolic or ellipsoidal mirrors and fast, low rotational inertia scanners, such as galvanometers. We derive theoretically and empirically relations between a variety of system parameters (e.g.,  NA, wavelength, focal length, telecentricity) and achievable fields of view (FOVs) and importantly show that intrinsic tilt aberrations donotrestrict FOV for many multiview imaging applications, contrary to conventional wisdom. Finally, we present strategies for avoiding spherical aberrations at obliquely illuminated flat boundaries. Our simple designs allow for high-speed multiangle imaging for microscopic, mesoscopic, and macroscopic applications.

    more » « less
  2. One of the top priorities in observational astronomy is the direct imaging and characterization of extrasolar planets (exoplanets) and planetary systems. Direct images of rocky exoplanets are of particular interest in the search for life beyond the Earth, but they tend to be rather challenging targets since they are orders-of-magnitude dimmer than their host stars and are separated by small angular distances that are comparable to the classicalλ<#comment/>/Ddiffraction limit, even for the coming generation of 30 m class telescopes. Current and planned efforts for ground-based direct imaging of exoplanets combine high-order adaptive optics (AO) with a stellar coronagraph observing at wavelengths ranging from the visible to the mid-IR. The primary barrier to achieving high contrast with current direct imaging methods is quasi-static speckles, caused largely by non-common path aberrations (NCPAs) in the coronagraph optical train. Recent work has demonstrated that millisecond imaging, which effectively “freezes” the atmosphere’s turbulent phase screens, should allow the wavefront sensor (WFS) telemetry to be used as a probe of the optical system to measure NCPAs. Starting with a realistic model of a telescope with an AO system and a stellar coronagraph, this paper provides simulations of several closely related regression models that take advantage of millisecond telemetry from the WFS and coronagraph’s science camera. The simplest regression model, called the naïve estimator, does not treat the noise and other sources of information loss in the WFS. Despite its flaws, in one of the simulations presented herein, the naïve estimator provides a useful estimate of an NCPA of∼<#comment/>0.5radian RMS (≈<#comment/>λ<#comment/>/13), with an accuracy of∼<#comment/>0.06radian RMS in 1 min of simulated sky time on a magnitude 8 star. Thebias-corrected estimatorgeneralizes the regression model to account for the noise and information loss in the WFS. A simulation of the bias-corrected estimator with 4 min of sky time included an NCPA of∼<#comment/>0.05radian RMS (≈<#comment/>λ<#comment/>/130) and an extended exoplanet scene. The joint regression of the bias-corrected estimator simultaneously achieved an NCPA estimate with an accuracy of∼<#comment/>5×<#comment/>10−<#comment/>3radian RMS and an estimate of the exoplanet scene that was free of the self-subtraction artifacts typically associated with differential imaging. The5σ<#comment/>contrast achieved by imaging of the exoplanet scene was∼<#comment/>1.7×<#comment/>10−<#comment/>4at a distance of3λ<#comment/>/Dfrom the star and∼<#comment/>2.1×<#comment/>10−<#comment/>5at10λ<#comment/>/D. These contrast values are comparable to the very best on-sky results obtained from multi-wavelength observations that employ both angular differential imaging (ADI) and spectral differential imaging (SDI). This comparable performance is despite the fact that our simulations are quasi-monochromatic, which makes SDI impossible, nor do they have diurnal field rotation, which makes ADI impossible. The error covariance matrix of the joint regression shows substantial correlations in the exoplanet and NCPA estimation errors, indicating that exoplanet intensity and NCPA need to be estimated self-consistently to achieve high contrast.

    more » « less
  3. The optical phaseϕ<#comment/>is a key quantity in the physics of light propagating through a turbulent medium. In certain respects, however, the statistics of the phasefactor,ψ<#comment/>=exp⁡<#comment/>(iϕ<#comment/>), are more relevant than the statistics of the phase itself. Here, we present a theoretical analysis of the 2D phase-factor spectrumFψ<#comment/>(κ<#comment/>)of a random phase screen. We apply the theory to four types of phase screens, each characterized by a power-law phase structure function,Dϕ<#comment/>(r)=(r/rc)γ<#comment/>(wherercis the phase coherence length defined byDϕ<#comment/>(rc)=1rad2), and a probability density functionpα<#comment/>(α<#comment/>)of the phase increments for a given spatial lag. We analyze phase screens with turbulent (γ<#comment/>=5/3) and quadratic (γ<#comment/>=2) phase structure functions and with normally distributed (i.e., Gaussian) versus Laplacian phase increments. We find that there is a pronounced bump in each of the four phase-factor spectraFψ<#comment/>(κ<#comment/>). The precise location and shape of the bump are different for the four phase-screen types, but in each case it occurs atκ<#comment/>∼<#comment/>1/rc. The bump is unrelated to the well-known “Hill bump” and is not caused by diffraction effects. It is solely a characteristic of the refractive-index statistics represented by the respective phase screen. We show that the second-orderψ<#comment/>statistics (covariance function, structure function, and spectrum) characterize a random phase screen more completely than the second-orderϕ<#comment/>counterparts.

    more » « less
  4. Three-dimensional (3D) refractive index (RI) tomography has recently become an exciting new tool for biological studies. However, its limitation to (1) thin samples resulting from a need of transmissive illumination and (2) small fields of view (typically∼<#comment/>50µ<#comment/>m×<#comment/>50µ<#comment/>m) has hindered its utility in broader biomedical applications. In this work, we demonstrate 3D RI tomography with a large field of view in opaque, arbitrarily thick scattering samples (unsuitable for imaging with conventional transmissive tomographic techniques) with a penetration depth of ca. one mean free scattering path length (∼<#comment/>100µ<#comment/>min tissue) using a simple, low-cost microscope system with epi-illumination. This approach leverages a solution to the inverse scattering problem via the general non-paraxial 3D optical transfer function of our quantitative oblique back-illumination microscopy (qOBM) optical system. A theoretical analysis is presented along with simulations and experimental validations using polystyrene beads, and rat and human thick brain tissues. This work has significant implications for the investigation of optically thick, semi-infinite samples in a non-invasive and label-free manner. This unique 3D qOBM approach can extend the utility of 3D RI tomography for translational and clinical medicine.

    more » « less
  5. Capable of imaging blood perfusion, oxygenation, and flow simultaneously at the microscopic level, multi-parametric photoacoustic microscopy (PAM) has quickly emerged as a powerful tool for studying hemodynamic and metabolic changes due to physiological stimulations or pathological processes. However, the low scanning speed poised by the correlation-based blood flow measurement impedes its application in studying rapid microvascular responses. To address this challenge, we have developed a new, to the best of our knowledge, multi-parametric PAM system. By extending the optical scanning range with a cylindrically focused ultrasonic transducer (focal zone,76µ<#comment/>m×<#comment/>4.5mm) for simultaneous acquisition of 500 B-scans, the new system is 112 times faster than our previous multi-parametric system that uses a spherically focused transducer (focal diameter, 40 µm) and enables high-resolution imaging of blood perfusion, oxygenation, and flow over an area of4.5×<#comment/>1mm2at a frame rate of 1 Hz. We have demonstrated the feasibility of this system in the living mouse ear. Further development of this system into reflection mode will enable real-time cortex-wide imaging of hemodynamics and metabolism in the mouse brain.

    more » « less