3 April 2013 Three-dimensional deconvolution microfluidic microscopy using a tilted channel
Author Affiliations +
J. of Biomedical Optics, 18(4), 040503 (2013). doi:10.1117/1.JBO.18.4.040503
We have developed a microfluidic device that enables the computation of three-dimensional (3-D) images of flowing samples. Using a microfluidic channel that is tilted along the optical axis, we record several progressively defocused images of the flowing sample as it passes across the focal plane. The resulting focal stack is then deconvolved to generate 3-D images. Experimental results on flowing yeast cells reveal both volume and surface profile information. The microfluidic channel eliminates the need for a precise translation stage to control defocusing and enables high sample throughput in an insulated, nontoxic, liquid environment. The experimental device can be implemented in all existing microscopes as a modified slide stage and is ideally suited for 3-D profiling in flow cytometers.
Pégard and Fleischer: Three-dimensional deconvolution microfluidic microscopy using a tilted channel



Microfluidic microscopy is a relatively new imaging technique based on guiding samples under the observation window by flowing them through a microfluidic channel. This includes conventional microscopes (e.g., objectives) with modified slide stages as well as lensless optofluidic microscopes.12.3 The fluid displacement increases sample throughput, provides control over sample sorting,4 and gives an additional degree of freedom for imaging. To date, in-plane flow has been used to generate two-dimensional (2-D) images by displacing the target over a one-dimensional (1-D) array of detectors,1 or to slightly displace the sample and generate subpixel-resolution images2 and holograms.3 Three-dimensional (3-D) imaging microfluidic devices have also been explored, but they have relied on multiple observation axes3,5 or multiple illumination angles.6

Here, we keep the well-known configuration for white light microscopy, with a single illumination source and objective lens, and reduce the experimental device to a straight microfluidic channel tilted along the optical axis (a lensless device is possible, but using an objective simplifies the design, facilitates integration with existing systems, and significantly improves image quality7). By taking successive images as the sample flows along the channel axis, we collect a set of intensity measurements at different propagation planes. There are two complementary uses for such data: phase retrieval in the focal plane via a transport-of-intensity equation8,9 and volume imaging by assembly of the optical slices.10,11 In the former case, the tilt should be shallow, to reduce large variations in the phase and better approximate the diffraction derivatives using two spatially separated images. In the latter case, a steeper tilt is necessary, so that the bulk of the sample falls within the focal region.


Experimental Setup

The experimental setup is shown in Fig. 1. A 500-μm-wide, 50-μm-deep microfluidic channel is etched on a glass slide and located at the focal plane of a standard wide-field microscope. The slide is tilted at a 15-deg angle, α, with respect to the optical axis of the microscope objective, chosen to represent a good compromise between magnification and axial defocusing. The channel is illuminated with incoherent white light, and a 25× magnified image is recorded by a video camera operating at a constant frame rate of 30frames/s. As a test case, a suspension of 15 μm yeast cells in glycerol is flowed through the microfluidic channel.

Fig. 1

(a) Experimental setup. A microfluidic channel is located under the objective of a wide field microscope and tilted along the optical axis by an angle of 15 deg. The sample flows into the channel at a constant velocity. Multiple images are recorded by a camera operating at a constant frame rate. (b) Samples are observed as they flow along the channel axis and pass across the focal plane.


The frame rate and flow speed are adjusted so that each sample will be recorded in 100 consecutive frames as it passes from one end of the window of observation to the other. A constant flow is maintained using a fixed pressure difference (a 50-cm hydrostatic water column) between the channel input and output, and the exposure time is adjusted to reduce the flow-induced blur below the resolution limit of the imaging system. Narrow channel depth and high fluid viscosity guarantee a Hagen–Poiseuille type laminar flow.12 In addition, the concentration of particles is lower than 250μL1 for easy separation and to reduce interactions between flowing particles.


Principles of Operation

The parabolic velocity flow ensures that particles flow at a constant velocity (v) along the channel axis but experience a shear-induced rotation at other points. As the acquisition of accurate focal stacks relies on the absence of rotation (or its compensation), the setup has been designed to minimize the effects of shear in all directions. Along the channel axis u, the absence of shear, v·u=0, is a property of the laminar flow. Along the y axis, shear-induced rotation effects are minimized by observing samples flowing in the middle part of a wide channel. Similarly, along the z axis, rotating objects are excluded by considering only particles flowing at the highest velocity in the middle of the channel, where the shear effects cancel. We note, however, that object rotation may be useful in other contexts, e.g., for multiple viewpoints, and is easily accessible by changing the injection point or imaging different parts of the flow.

Focal stacks are generated by tracking samples flowing into the channel, as shown in Fig. 2. The background noise is subtracted from the signal, and the zero value of the signal (in gray) corresponds to the nominal transparency of the free-running fluid. Let T be the frame recording period, and δ=Tv the object displacement along the channel axis between two frames. Let In(x,y) be the intensity of frame n. The focal stack, S, is given by



Fig. 2

Focal stacks are generated by observing samples as they flow in the tilted microfluidic channel. Constant flow velocity and frame rate enable the accumulation of 100 frames progressively defocused along the z axis. An object-tracking algorithm based on defocusing invariant properties of the center of gravity of the image is used for better accuracy. (a) Measured iso-intensity contours of yeast cells through focus (36μm). Insets show direct images of cells. (b) Normalized intensity of focal stack with background subtracted.


The extraction of focal stacks is based not only on the estimated position of the sample but also on the point spread function of the imaging system, P. In general, the function P can depend on position, which may vary due to fluctuations in the flow. We compensate for this by using a simple method of particle tracking. Assuming that P satisfies the symmetry condition P(x,y,z)=P(x,y,z), the first moment of the recorded image,


becomes a defocusing invariant.13 In our reconstruction algorithm, we compute the center of mass Gn for each frame in the focal stack. In the (x, y) plane, the position of Gn allows an accurate alignment of frames onto each other. Along the optical axis, the relation zGn·ex tan(α) allows a precise determination of the amount of defocusing.

The reconstruction of the volume absorption distribution of the object, O, relies on solving the well-known deconvolution problem using the experimentally measured focal stack, S, and point spread function, P, in 3-D space. O satisfies the volume integral


where r=x.x+y.y+z.z. In the ideal treatment of this problem, the optical transfer function, F[P], has only nonzero terms, so that the solution of Eq. (3) is given by



Unfortunately, this solution is known to be extremely sensitive to experimental noise.14 A variety of noise-reducing techniques have been developed to recover the object if the point spread function is not known, e.g., maximum likelihood estimation15 and blind deconvolution.16 In these methods, which in general can be quite complex, the point spread function is guessed instead of measured. Adaptive measures,17 which are particularly suited to flowing objects, or engineered point spread functions18 can be used as well. In the experiment here, we directly measure the point spread function of the microscope by acquiring the focal stack of sub-micron-based reference particles. For this calibration, we use a suspension of 800 nm dyed polystyrene beads in glycerol, with flow velocity, frame rate, and exposure conditions identical to those of the samples. We record a focal stack as we track one of the flowing reference particles, align the defocused images using the defocusing invariant Gn, and check that the field of observation is clear of other flowing objects. The resulting stack, centered in the window of computation and normalized to unitary absorption, represents the point spread function of the microfluidic microscope for this particular tilt angle. Although not done here, these particles can be embedded in the flow with the samples themselves, acting as real-time reference points for changing conditions and/or shear compensation.

Even with a known point spread function, the inversion (4) is sensitive to zeros and noise in the measurement. To compensate for this, we use a Wiener deconvolution filter.19 To implement the filter, we return to Eq. (3) and consider explicitly an additive noise term, N, which we assume to be independent from the signal. Equation (3) becomes:


The principle of the Wiener filter is to find the best deconvolution operator, D, so that the retrieved object OR=[DS(r)] minimizes the RMS reconstruction error ER=|OOR|2.

The Wiener solution, representing an optimal compromise between noise and resolution, is given by


where ϵ=(|N|2/|S|2)=1.2*103 is a regularization constant corresponding to the inverse value of the signal-to-noise ratio. To implement this regularization, the signal intensity is normalized to 1 and the root-mean-square value of the noise (background intensity) is measured in an empty area near the flowing object.



Experimental results on flowing yeast cells are shown in Fig. 3. The focal stack, S, and the point spread function, P, are each processed using Eq. (6). 3-D data is retrieved, enabling the construction of iso-level surface contours. Figure 3(a) shows a projection of these contours along the optical axis, revealing small-scale surface features (12μm) which are clearly resolved (though smoothed somewhat by the regularization process). Most likely, they are due to early-stage budding, though many other factors can contribute to their morphology.20 Figure 3(b) gives contours from the side, showing that all the cells lie in the same vertical plane (a result of the controlled injection) and that each cell has flat side walls. Such deformation is common in flowing cells21 and is often used as a diagnostic.22 As stated above, both the observed shapes and surface profiles are minimally affected by noise in the deconvolution process, as Wiener regularization optimally smooths the reconstructed profiles. We emphasize that many details that are hidden in standard imaging using 2-D projections, such as cell orientation, 3-D shape, and surface roughness, are readily apparent in the volume images here.

Fig. 3

The three-dimensional (3-D) structure of flowing yeast cells is reconstructed digitally using Eq. (6). (a) An iso-level surface shows subcellular structures at the surface of the cellular membrane. (b) A three-quarter view of the samples, showing smooth structure of the cell walls and relative height in the channel.




In conclusion, we have introduced a new 3-D imaging technique based on recording a focal stack from samples flowing through a tilted microfluidic channel. As an example, we presented in-flow volume imaging of aggregated yeast cells, showing their size, position, orientation, and subcellular surface features. The method can be implemented on any existing microfluidic microscope and is ideally suited for 3-D profiling in flow cytometers.


The authors are grateful to their collaborators from the Complex Fluids Group at Princeton University.



X. Henget al., “Optofluidic microscopy: a method for implementing a high resolution optical microscope on a chip,” Lab Chip 6(10), 1274–1276 (2006).LCAHAM1473-0197http://dx.doi.org/10.1039/b604676bGoogle Scholar


G. Zhenget al., “Sub-pixel resolving optofluidic microscope for on-chip cell imaging,” Lab Chip 10(22), 3125–3129 (2010).LCAHAM1473-0197http://dx.doi.org/10.1039/c0lc00213eGoogle Scholar


W. BisharaH. ZhuA. Ozcan, “Holographic optofluidic microscopy,” Opt. Express 18(26), 27499–27510 (2010).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.18.027499Google Scholar


S. S. H. TsaiI. M. GriffithsH. A. Stone, “Microfluidic immunomagnetic multi-target sorting—a model for controlling deflection of paramagnetic beads,” Lab Chip 11(15), 2577–2582 (2011).LCAHAM1473-0197http://dx.doi.org/10.1039/c1lc20153kGoogle Scholar


S. Yanget al., “Stereoscopic optofluidic on-chip microscope,” in Proc. IEEE Winter Topicals, pp. 91–92, IEEE, Keystone, Colorado (2011).Google Scholar


S. Isikmanet al., “Lens-free optical tomographic microscope with a large imaging volume on a chip,” Proc. Natl. Acad. Sci. U. S. A. 108(18), 7296–7301 (2011).PNASA60027-8424http://dx.doi.org/10.1073/pnas.1015638108Google Scholar


N. PégardJ. Fleischer, “Optimizing holographic data storage using a fractional Fourier transform,” Opt. Lett. 36(13), 2551–2553 (2011).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.36.002551Google Scholar


M. R. Teague, “Deterministic phase retrieval: a Greens function solution,” J. Opt. Soc. Am. 73(11), 1434–1441 (1983).JOSAAH0030-3941http://dx.doi.org/10.1364/JOSA.73.001434Google Scholar


S. GorthiE. Schonbrun, “Phase imaging flow cytometry using a focus-stack collecting microscope,” Opt. Lett. 37(4), 707–709 (2012).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.37.000707Google Scholar


A. Erhardtet al., “Reconstructing 3-D light-microscopic images by digital image processing,” Appl. Opt. 24(2), 194–200 (1985).APOPAI0003-6935http://dx.doi.org/10.1364/AO.24.000194Google Scholar


N. PégardJ. Fleischer, “3D microfluidic microscopy using a tilted channel,” in Biomedical Optics, Optical Society of America, Miami, Florida (2012).Google Scholar


H. StoneS. Kim, “Microfluidics: basic issues, applications, and challenges,” AIChE J. 47(6), 1250–1254 (2001).AICEAC0001-1541http://dx.doi.org/10.1002/aic.690470602Google Scholar


M. R. Teague, “Irradiance moments: their propagation and use for unique retrieval of phase,” J. Opt. Soc. Am. 72(9), 1199–1209 (1982).JOSAAH0030-3941http://dx.doi.org/10.1364/JOSA.72.001199Google Scholar


A. Tikhonovet al., Numerical Methods for the Solution of Ill-Posed Problems, Vol. 328, Springer, New York (1995).Google Scholar


W. Richardson, “Bayesian-based iterative method of image restoration,” J. Opt. Soc. Am. 62(1), 55–59 (1972).JOSAAH0030-3941http://dx.doi.org/10.1364/JOSA.62.000055Google Scholar


T. ChanC.-K. Wong, “Total variation blind deconvolution,” IEEE Trans. Image Process. 7(3), 370–375 (1998).IIPRE41057-7149http://dx.doi.org/10.1109/83.661187Google Scholar


W. Donget al., “Image deblurring and super-resolution by adaptive sparse domain selection and adaptive regularization,” IEEE Trans. Image Process. 20(7), 1838–1857 (2011).IIPRE41057-7149http://dx.doi.org/10.1109/TIP.2011.2108306Google Scholar


S. Pavaniet al., “Three-dimensional, single-molecule fluorescence imaging beyond the diffraction limit by using a double-helix point spread function,” Proc. Natl. Acad. Sci. U. S. A. 106(9), 2995–2999 (2009).PNASA60027-8424http://dx.doi.org/10.1073/pnas.0900245106Google Scholar


C. ChatwinR. Wang, Frequency Domain Filtering Strategies for Hybrid Optical Information Processing, Research Studies Press, Baldock, United Kingdom (1996).Google Scholar


Y. Ohyaet al., “High-dimensional and large-scale phenotyping of yeast mutants,” Proc. Natl. Acad. Sci. U. S. A. 102(52), 19015–19020 (2005).PNASA60027-8424http://dx.doi.org/10.1073/pnas.0509436102Google Scholar


M. Abkarianet al., “Cellular-scale hydrodynamics,” Biomed. Mater. 3(3), 034011 (2008).BMATEM1748-6041http://dx.doi.org/10.1088/1748-6041/3/3/034011Google Scholar


C. Westendorfet al., “Live cell flattening—traditional and novel approaches,” BMC Biophys. 3(1), 9 (2010).http://dx.doi.org/10.1186/1757-5036-3-9Google Scholar

Nicolas C. Pegard, Jason W. Fleischer, "Three-dimensional deconvolution microfluidic microscopy using a tilted channel," Journal of Biomedical Optics 18(4), 040503 (3 April 2013). http://dx.doi.org/10.1117/1.JBO.18.4.040503
Submission: Received ; Accepted


3D image processing



Point spread functions


3D metrology

Back to Top