Optimal illumination scheme for isotropic quantitative differential phase contrast microscopy

Yao Fan; Jiasong Sun; Qian Chen; Xiangpeng Pan; Lei Tian; Chao Zuo

doi:doi:10.1364/PRJ.7.000890

Photonics Research, 2019, 7 (8): 08000890, Published Online: Jul. 25, 2019

Optimal illumination scheme for isotropic quantitative differential phase contrast microscopy Download： 720次

论文大纲

Yao Fan ^1,2,3†Jiasong Sun ^1,2,3†Qian Chen ^1,2,5Xiangpeng Pan ^1,2,3Lei Tian ⁴Chao Zuo ^1,2,3,*

Author Affiliations

¹ School of Electronic and Optical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China

² Jiangsu Key Laboratory of Spectral Imaging & Intelligent Sense, Nanjing University of Science and Technology, Nanjing 210094, China

³ Smart Computational Imaging (SCI) Laboratory, Nanjing University of Science and Technology, Nanjing 210094, China

⁴ Department of Electrical and Computer Engineering, Boston University, Boston, Massachusetts 02215, USA

⁵ e-mail: chenqian@njust.edu.cn

Abstract

Differential phase contrast microscopy (DPC) provides high-resolution quantitative phase distribution of thin transparent samples under multi-axis asymmetric illuminations. Typically, illumination in DPC microscopic systems is designed with two-axis half-circle amplitude patterns, which, however, result in a non-isotropic phase contrast transfer function (PTF). Efforts have been made to achieve isotropic DPC by replacing the conventional half-circle illumination aperture with radially asymmetric patterns with three-axis illumination or gradient amplitude patterns with two-axis illumination. Nevertheless, the underlying theoretical mechanism of isotropic PTF has not been explored, and thus, the optimal illumination scheme cannot be determined. Furthermore, the frequency responses of the PTFs under these engineered illuminations have not been fully optimized, leading to suboptimal phase contrast and signal-to-noise ratio for phase reconstruction. In this paper, we provide a rigorous theoretical analysis about the necessary and sufficient conditions for DPC to achieve isotropic PTF. In addition, we derive the optimal illumination scheme to maximize the frequency response for both low and high frequencies (from 0 to

2 {NA}_{obj}

) and meanwhile achieve perfectly isotropic PTF with only two-axis intensity measurements. We present the derivation, implementation, simulation, and experimental results demonstrating the superiority of our method over existing illumination schemes in both the phase reconstruction accuracy and noise-robustness.

1. INTRODUCTION

Quantitative phase imaging (QPI), which provides phase information about the refractive index distribution of transparent specimens, has drawn much attention in both optical and biomedical research ^{[13" target="_self" style="display: inline;">–3]}. The major advantage of QPI over conventional intensity imaging or fluorescence microscopy is that it requires no exogenous contrast agents (e.g., dyes or fluorescent protein) to enhance the contrast of the microscopic image, which enables label-free and stain-free optical imaging of live biological specimens in vitro ^[4,5]. The most common QPI methods are based on interferometry and holography with coherent illumination and a reference beam, making them expensive and sensitive to misalignment, vibrations, and speckle noise ^{[68" target="_self" style="display: inline;">–8]}. To overcome these limitations, non-interferometric QPI approaches using partially coherent illumination have been developed, such as transport-of-intensity equation (TIE) ^{[2,9–14" target="_self" style="display: inline;">–14]}, differential phase contrast microscopy (DPC) ^{[15–20" target="_self" style="display: inline;">–20]}, and Fourier ptychographic microscopy (FPM) ^{[21–25" target="_self" style="display: inline;">–25]}. TIE is a well-established non-interferometric phase retrieval approach, which enables the QPI of a transparent sample simply by measuring the intensities at multiple axially displaced planes ^{[2,9,11–14" target="_self" style="display: inline;">–14]}. The advantages of the TIE approach are that it is fully compatible with widely available bright-field microscopy hardware and able to offer an imaging resolution up to the incoherent resolution limit (two times better than the coherent diffraction limit) under matched annular illumination ^[26,27]. Without moving the position of a sample, DPC and FPM approaches retrieve the complex field of the sample by using asymmetric illuminations. In FPM, a set of low-resolution (LR) intensity images corresponding to different illumination angles, with the resolution determined by the numerical aperture (NA) of the objective lens, are acquired ^{[21,22,24,25]}. These LR intensity images are iteratively combined together in the Fourier domain, resulting in a wide-field, high-resolution complex image with the synthesized resolution determined by the sum of the objective lens and illumination NAs ^{[21–25" target="_self" style="display: inline;">–25,28,29]}. DPC, however, achieves phase recovery by using only four images with asymmetric illuminations in opposite directions ^[18,19]. It converts invisible sample phases into measurable intensity by shifting the sample’s spectrum in Fourier space to theoretically achieve a resolution of twice the coherent diffraction limit ^[18,19]. Assuming a linearized model for a weakly scattering sample, the DPC phase retrieval problem becomes a single-step deconvolution process using the phase contrast transfer function (PTF). By implementing DPC with a programable LED array ^[19,20,30] or a programable LCD panel ^[31], we are able to realize dynamic QPI along arbitrary axes of asymmetry, without any mechanical moving parts.

To recover the quantitative phase information of a weakly scattering sample, at least two complementary source patterns are required in DPC. However, the resultant PTF is anti-symmetric and zero at all spatial frequencies along the axis of asymmetry, which may lead to significant phase reconstruction artifacts if not properly handled. So in general, the illumination of DPC is designed with two-axis half-circle amplitude patterns, that is, four patterns (top, bottom, left, right half-circles) are used to avoid missing frequencies. However, artifacts of phase reconstruction still cannot be completely avoided, since DPC’s PTF is not circularly symmetric with only two-axis measurements. Recently, efforts have been made toward developing high-speed, or even single-shot QPI mechanisms, such as motion compensation and color-coded DPC based on wavelength multiplexing ^{[30,3234" target="_self" style="display: inline;">–34]}. Other efforts have been made to improve the phase reconstruction quality of DPC by optimizing the illumination scheme, phase reconstruction algorithm, and compensating for aberrations, and so on ^{[35–39" target="_self" style="display: inline;">–39]}. Among these methods, optimizing illumination becomes the most direct and effective approach to improve the image quality of DPC since the PTF of DPC is directly determined by illumination. By using three-axis radially asymmetric patterns ^[35] or two-axis gradient amplitude patterns ^[36] instead of more time-consuming intensity measurements along different illumination angles (12-axis) ^[20], the isotropic DPC is achieved with higher imaging efficiency. The machine learning is also introduced to the DPC illumination design, which improves the phase recovery accuracy, while reducing the raw intensity measurements at the same time ^[38]. Although the transfer response and isotropy of the PTF improved, the underlying theoretical mechanism of isotropic PTF has not been explored, so the optimal illumination scheme cannot be determined. Furthermore, the frequency responses of the PTFs under these engineered illuminations have not been fully optimized, leading to suboptimal phase contrast and signal-to-noise ratio (SNR) for phase reconstruction.

In this paper, we improve on these works by providing a rigorous theory for achieving isotropic DPC, where a new optimal illumination scheme is derived. The major advantages of the new illumination scheme are twofold. First, it is able to produce a circularly symmetrical PTF with only two-axis intensity measurements under partially coherent condition. Thus, it is expected to achieve high-quality phase reconstruction with isotropic transverse resolution and SNR by using only four intensity measurements. Second, the resulting PTF achieves a broadband frequency coverage for partially coherent imaging (from 0 to $2 {NA}_{obj}$ ) with a smooth and significantly enhanced response for both low and high frequencies, which alleviates the ill-posedness of the PTF inversion.

2. OPTIMAL ILLUMINATION SCHEME DESIGN

2.1 A. Derivation of Generalized Illumination Function for Isotropic Differential Phase Contrast

Consider a weak-scattering object with complex transmission function $t (r) = e^{- a (r) + i ϕ (r)}$ , illuminated by an oblique plane wave with a uniform intensity distribution $S (u_{j})$ ( $u_{j}$ denotes the spatial frequency of the tilted illumination). Invoking the weak object approximation $t (r) \approx 1 - a (r) + i ϕ (r)$ ^[40,41], the intensity spectrum of the bright-field image under oblique illumination can be separated into three terms ^[42], including the background, absorption contrast, and phase contrast terms (ignore systematic aberrations): $I_{j} (u) = S (u_{j}) δ (u) {| P (u_{j}) |}^{2} - S (u_{j}) A (u) [P^{*} (u_{j}) P (u + u_{j}) + P (u_{j}) P^{*} (u - u_{j})] + i S (u_{j}) Φ (u) [P^{*} (u_{j}) P (u + u_{j}) - P (u_{j}) P^{*} (u - u_{j})],$ (1)where $P (u)$ denotes the pupil function of the objective lens (assuming it is an ideal low-pass filter with a cutoff frequency of $\frac{{NA}_{obj}}{λ}$ ), and $A (u)$ and $Φ (u)$ are the Fourier transform of the sample’s absorption and phase function, respectively. To generate the phase contrast image $I_{l r}^{DPC}$ (e.g., the left–right DPC) along a specific direction of the phase gradient, a pair of images where each $S (u)$ has complementary gradient vector are used to calculate $I_{l r}^{DPC} = \frac{I_{r} - I_{l}}{I_{r} + I_{l}}$ ^[18,41]. Since two illumination patterns are symmetrical along the same axial direction (e.g., $x$ -axis in left–right DPC), the background term and absorption contrast term are cancelled, leaving only the phase contrast term. So, the corresponding PTF of the left–right DPC can be expressed as ${PTF}_{l r} (u) = \frac{\iint S_{l r} (u_{j}) [P^{*} (u_{j}) P (u + u_{j}) - P (u_{j}) P^{*} (u - u_{j})] d^{2} u_{j}}{\iint | S_{l r} (u_{j}) | {| P (u_{j}) |}^{2} d^{2} u_{j}} .$ (2)It can be found from Eq. (2) that once the optical configuration of the microscope is fixed (the pupil function of the objective lens is predefined), the PTF is fully determined by the illumination function ^[43].

To reconstruct the sample’s quantitative phase information from Eq. (2), we can solve the inverse problem with a single-step deconvolution ^[20]. Tikhonov regularization parameter $β$ is often introduced in the denominator to avoid singularity in the PTF inversion ^[44]: $ϕ (r) = F^{- 1} {\frac{\sum_{i} [{PTF}_{i}^{*} (u) \cdot I_{i}^{DPC} (u)]}{\sum_{i} {| {PTF}_{i}^{*} (u) |}^{2} + β}},$ (3)where ${PTF}_{i}^{*} (u)$ denotes complex conjugation of PTF along different axial directions. The denominator term $\sum_{i} {| {PTF}_{i}^{*} (u) |}^{2}$ represents the synthetic square of amplitude of the multi-axial PTFs, which can be used to indicate the degree of isotropy for DPC imaging. For simplicity, we use a shorthand notation $C (u)$ for this term in the following analysis.

The following derivation about isotropic DPC is performed in the polar coordinate instead of the Cartesian coordinate; that is, we use $(ρ, θ)$ instead of $(x, y)$ , where $ρ$ and $θ$ respectively represent the radius and the polar angle. This is because the optical systems have circular symmetry, and most circularly symmetric functions are separable in their polar coordinates, i.e., we can write $S_{l r} (ρ, θ)$ and $S_{u d} (ρ, θ)$ as a product of two one-dimensional functions about $ρ$ and $θ$ , respectively: $S_{l r} (ρ, θ) = L (ρ) M (θ), S_{u d} (ρ, θ) = L (ρ) N (θ) .$ (4)Since the phase contrast image is generated by complementary illumination in each axis measurement $S_{l r} = S_{r} - S_{l}$ , the illumination function $S_{l r} (ρ, θ)$ is an even function about $θ$ . Similarly, the illumination function $S_{u d} (ρ, θ)$ for the up–down axis is an odd function about $θ$ . Thus, $M (θ)$ and $N (θ)$ can be further expanded in the Fourier series defined on ( $- π$ , $π$ ]: $M (θ) = \sum_{n = 1}^{\infty} a_{n} \cos (n θ), N (θ) = \sum_{n = 1}^{\infty} b_{n} \sin (n θ) .$ (5)For general illumination and aperture function, the PTF can be calculated by integrating the overlapping areas between the objective pupil function (with its center at the origin $O$ ) and its off-axis version (with its center at the point $Q$ ), as illustrated by the red and blue regions in Fig. 1. This is because the illumination falling in the red regions can ensure the point $Q$ is within $P (u + u_{j}) = 1$ , and the illumination falling in the blue regions can ensure the point $Q (u)$ is within $P (u - u_{j}) = 1$ . However, it should be noted that when the illumination angle is close to the central axis of the objective lens, the two red and blue regions corresponding to $P (u + u_{j})$ and $P (u - u_{j})$ will be partially cancelled out by each other. So the integral interval should be partitioned according to the location of point $Q$ [see Figs. 1(a) and 1(b)], which can be represented as (taking left–right axis illumination for example) ${PTF}_{l r} (ρ, θ) = {\begin{matrix} \frac{2 \int_{ρ - {NA}_{obj}}^{{NA}_{obj}} \int_{θ - α}^{θ + α} S_{l r} (ξ, ε) d ε d ξ}{\int_{0}^{{NA}_{obj}} \int_{0}^{2 π} | S_{l r} (ξ, ε) | d ε d ξ} & {NA}_{obj} \leq ρ \leq 2 {NA}_{obj} \\ \frac{2 \int_{{NA}_{obj} - ρ}^{{NA}_{obj}} \int_{θ - α}^{θ + α} S_{l r} (ξ, ε) d ε d ξ}{\int_{0}^{{NA}_{obj}} \int_{0}^{2 π} | S_{l r} (ξ, ε) | d ε d ξ} & 0 \leq ρ < {NA}_{obj} \end{matrix} .$ (6)Substituting Eqs. (4) and (5) into Eq. (6) gives the PTF along the left–right axis direction. Although the mathematical formulas appear complicated, we should keep in mind that to achieve isotropic DPC, it is only required that the synthetic square of amplitude of the multi-axial PTFs $C (ρ, θ)$ [ $C (ρ, θ) = {| {PTF}_{l r} (ρ, θ) |}^{2} + {| {PTF}_{u d} (ρ, θ) |}^{2}$ ] is only a function of $ρ$ . Close inspection of Eq. (6) reveals that the isotropy cannot be achieved if there are cross-trigonometric terms in $C (ρ, θ)$ , suggesting there should be only one single-harmonic component in the Fourier series expansion of Eq. (5), that is (see Appendix A for detailed derivation): $\begin{matrix} S_{l r} (ρ, θ) = L (ρ) \cos (n θ) \\ S_{u d} (ρ, θ) = L (ρ) \sin (n θ) \end{matrix} (n = 1, 3, 5, \dots) .$ (7)Note that in Eq. (7), we neglect the unimportant constant factors. It should be further noted that when $n$ is even, the illumination pattern is centrosymmetric so that the resulting PTF will be completely canceled out. Therefore, $n$ should be an odd number $(n = 1, 3, 5, \dots)$ to guarantee a valid non-zero PTF. Equation (7) is the main result of our work, which provides the necessary and sufficient conditions for DPC to achieve isotropic PTF.

Schematic diagram of the integral for PTF along the left–right axis in the polar coordinate system. (a) The radius ρ of the point Q is in the range of NAobj≤ρ≤2NAobj. (b) The radius ρ of the point Q is in the range of 0≤ρ<NAobj.

Fig. 1. Schematic diagram of the integral for PTF along the left–right axis in the polar coordinate system. (a) The radius $ρ$ of the point $Q$ is in the range of ${NA}_{obj} \leq ρ \leq 2 {NA}_{obj}$ . (b) The radius $ρ$ of the point $Q$ is in the range of $0 \leq ρ < {NA}_{obj}$ .

下载图片查看所有图片

2.2 B. Optimal Illumination Scheme

Since no restrictions were imposed on the form of the function $L (ρ)$ , $L (ρ)$ can be any function of $ρ$ without affecting the isotropy of the DPC’s PTF. Thus, we can optimize the function $L (ρ)$ to improve the frequency coverage and response of the corresponding PTF. In our recent studies, an annular source has been demonstrated that optimizes the PTF of DPC ^[45], TIE ^[26,27], and FPM ^[28,42]. In this paper, the annular source is also used in the generalized illumination function of isotropic DPC to optimize the transfer response of PTF. In Appendix B, we compare the resultant PTFs of three different functions $L (ρ)$ : constant weight, linear weight, and Kronecker delta weight functions with a fixed $n = 1$ by simulation. The results suggest that the delta weight function (a thin annulus) produces the PTF with the highest response at almost all frequencies from 0 to $2 {NA}_{obj}$ .

To further analyze the influence of the thickness of the annulus on the PTF, we illustrate the corresponding PTFs of different annular thickness by fixing the NA of the outer circles to be ${NA}_{obj}$ and by only changing the thickness of the annulus. As might be expected, the phase contrast is gradually reduced as the annulus width increases, especially at low and high frequencies. This is because the paraxial illumination does not produce low-frequency phase contrast, and only illumination matching the objective ${NA}_{obj}$ can produce the strong response at all frequencies. Specific simulations and comparisons can be found in Appendix C. When the thickness of the annulus is extremely small ( $σ = \frac{{NA}_{ill}}{{NA}_{obj}} \to 0$ ), the phase response finally approaches that of the constant weight function $L (ρ) = 1$ . From the results in Appendix C, it can be deduced that we should choose the diameter of the annulus to be equal to that of the objective pupil and make its thickness as small as possible [ $δ (ρ - {NA}_{obj})$ , where $δ (ρ)$ is the delta function] to optimize the response of PTF.

Finally, we study the effect of the number $n$ in Eq. (7) on the PTF. Three odd $n$ numbers are selected to generate three illumination patterns and corresponding PTFs (see Appendix D). It is shown that when $n$ increases ( $n = 3, 5$ ), not only the PTF response is significantly attenuated, but also the number of zero crossings in $C (ρ, θ)$ increases. In such cases, the reconstruction phase can be severely distorted due to the ill-posedness of the PTF inversion. Thus, $n = 1$ provides the optimal PTF for DPC with the strongest response and no zero crossings from 0 to $2 {NA}_{obj}$ . Based on the above analysis, the optimal illumination scheme for isotropic DPC can be represented as $S_{l r} (ρ, θ) = δ (ρ - {NA}_{obj}) \cos θ, S_{u d} (ρ, θ) = δ (ρ - {NA}_{obj}) \sin θ .$ (8)Meanwhile, we can also give the analytical expressions of PTF under the optimal illumination scheme: ${PTF}_{l r} (ρ, θ) = \sin α \cos θ, {PTF}_{u d} (ρ, θ) = \sin α \sin θ .$ (9)Based on the geometric relationship of the isosceles triangle BOQ in Fig. 1, $α$ is determined by $\cos α = \frac{ρ}{2 {NA}_{obj}}$ . In this case, $C (ρ, θ)$ can be calculated as $C (ρ, θ) = 1 - \frac{ρ^{2}}{4 {NA}_{obj}^{2}} .$ (10)Equation (10) indicates that $C (ρ, θ)$ is only related to $ρ$ , which means that the corresponding PTF of DPC is isotropic.

3. IMAGING PERFORMANCE OF OPTIMAL ILLUMINATION SCHEME

To verify the isotropy of the optimal illumination scheme, we numerically simulated PTFs and $C (ρ, θ)$ for four DPC illumination patterns under two-axis illumination, including three states of the arts, namely uniform illumination ^[20], (two-axis) radial illumination ^[35], gradient amplitude illumination ^[36], and the optimal illumination proposed in this work, as shown in Fig. 2. In Figs. 2(b1)–2(b4), we show PTFs along the left–right axis under these different illumination patterns. It can be observed that the PTF under optimal illumination scheme has a smooth and significantly enhanced response at almost all frequencies of the theoretical bandwidth of the entire partially coherent imaging (from 0 to $2 {NA}_{obj}$ ). These can be seen more clearly from the amplitude of $C (ρ, θ)$ in Figs. 2(c1)–2(c4), where the $C (ρ, θ)$ under the optimal illumination scheme has obviously enhanced values at the low-frequency components near the zero frequency and the high-frequency components approaching $2 {NA}_{obj}$ .

PTF and C(ρ,θ) with four illumination patterns. (a1)–(a4) Four illumination patterns. (b1)–(b4) PTFs along the left–right axis. (c1)–(c4) C(ρ,θ) with two-axis illumination. (d) Quantitative curves of C(ρ,θ) along the black straight line under the four illumination patterns. (e1)–(e3) Quantitative curves of C(ρ,θ) under the four illumination patterns on three radii.

Fig. 2. PTF and $C (ρ, θ)$ with four illumination patterns. (a1)–(a4) Four illumination patterns. (b1)–(b4) PTFs along the left–right axis. (c1)–(c4) $C (ρ, θ)$ with two-axis illumination. (d) Quantitative curves of $C (ρ, θ)$ along the black straight line under the four illumination patterns. (e1)–(e3) Quantitative curves of $C (ρ, θ)$ under the four illumination patterns on three radii.

下载图片查看所有图片

To quantitatively compare these four PTFs, the black cross-section in Fig. 2(c1) is used to characterize the amplitude of $C (ρ, θ)$ . As shown in Fig. 2(d), the optimal illumination scheme has the maximum phase contrast at almost all frequencies from 0 to $2 {NA}_{obj}$ . Moreover, the peak values of $C (ρ, θ)$ for the other three illumination patterns are below 0.9, while they can reach 1 under the optimal illumination scheme. The strong phase contrast under the optimal illumination can be finally converted to the quantitative phase images by the PTF inversion, resulting in high-quality reconstructions with a uniform background and an improved resolution. To compare the degree of isotropy, we further plotted the values along three different concentric circles within $C (ρ, θ)$ in Figs. 2(e1)–2(e3). A constant frequency response can be obtained along the circle under the gradient amplitude illumination and the optimal illumination scheme, suggesting that $C (ρ, θ)$ obtained under these two illumination patterns are isotropic, while the uniform illumination ^[20] and the radial illumination ^[35] cannot generate isotropic PTFs under two-axis illumination due to their fluctuant frequency responses along the concentric circles.

The proposed optimal illumination scheme is then compared with different illumination patterns based on simulations. The Siemens star image [shown in Fig. 3(a)] ^[46] is used as an example phase object that is defined on a grid with $244 \times 244$ pixels with a pixel size of $0.2 μm \times 0.2 μm$ . The wavelength of the illumination is 525 nm, and the ${NA}_{obj}$ is 0.40. For such an imaging configuration, the ideal phase imaging resolution that can be achieved is 656 nm ( $λ / 2 {NA}_{obj}$ ), which is also shown in Fig. 3(b). For partially coherent image calculation, we use the Abbe’s method, in which each sub-image corresponding to point source in the aperture plane is superimposed at the image plane to generate captured images of DPC. To simulate the noise effect, each DPC image is corrupted by Gaussian noise with a standard deviation of 0.0003. In Fig. 3, we compare the phase retrieval results of different illumination patterns for different regularization conditions. For the case without regularization ( $β$ is an infinitesimal), the noise corresponding to the frequency components with extremely weak PTF responses was amplified, resulting in grainy artifacts superimposed on the reconstructed phases by using uniform illumination ^[20], (two-axis) radial illumination ^[35], and gradient amplitude illumination ^[36] [Figs. 3(c1)–3(c3)]. To stabilize the deconvolution process, a regularization parameter $β$ is generally introduced to the denominator of Eq. (3) to suppress the noise effect. As shown in Figs. 3(d1)–3(d3), when $β = 0.2$ , the grainy artifacts are significantly reduced under these three illumination patterns, but meanwhile, the low-frequency phase values are underestimated, and high-frequency features are significantly attenuated. Thus, a more properly chosen regularization parameter is required for these illumination patterns to achieve reliable phase reconstructions under noisy conditions. In contrast, the proposed optimal illumination scheme can always obtain accurate and high-quality reconstructed phases with any regularization parameters attributed to its significantly enhanced PTF response from the 0 to $2 {NA}_{obj}$ range [Fig. 2(b4)], as shown in Figs. 3(c4) and 3(d4). The root-mean-square error (RMSE) values for the optimal illumination DPC are compared with the other illumination patterns to quantitatively measure the SNR of these reconstructed phases. It can be found that the reconstructed phase under the optimal illumination always achieves the lowest RMSE values.

$Simulation results with different regularization parameters under four illumination patterns. (a) Original phase image. (b) Diffraction limit phase image of DPC (2NAobj). (c1)–(c4), (d1)–(d4) Phase results with regularization parameters of 0 and 0.2. (e) Phase values along three small circles corresponding to different spatial frequencies evenly distributed from 0 to 2NAobj. (f) Phase values along a small circle of the same radius in (d1)–(d4) under four illumination patterns.$

Fig. 3. Simulation results with different regularization parameters under four illumination patterns. (a) Original phase image. (b) Diffraction limit phase image of DPC ( $2 {NA}_{obj}$ ). (c1)–(c4), (d1)–(d4) Phase results with regularization parameters of 0 and 0.2. (e) Phase values along three small circles corresponding to different spatial frequencies evenly distributed from 0 to $2 {NA}_{obj}$ . (f) Phase values along a small circle of the same radius in (d1)–(d4) under four illumination patterns.

下载图片查看所有图片

To verify the isotropy of the DPC phase reconstruction under the optimal illumination scheme, phase values along three small circles corresponding to different spatial frequencies evenly distributed from 0 to $2 {NA}_{obj}$ are extracted and plotted in Fig. 3(e) (the black curve denotes the ground truth, and the red, green, and blue curves represent the low frequency, middle frequency, and high frequency, respectively). From these curves, the line widths corresponding to different frequency components can be clearly distinguished at all angles, which indicates that optimal illumination achieves an isotropic resolution at all spatial frequencies. In addition, it can be observed that the reconstructed phase contrast under optimal illumination is close to the ground truth phase. Furthermore, to quantitatively compare the resolution of DPC reconstruction results of these four illumination patterns, we then extracted and plotted the phase values along a small circle of the same radius in Figs. 3(d1)–3(d4). As shown in Fig. 3(f), the uniform, radial, and gradient amplitude illuminations all produce results with underestimated phase contrast (phase amplitude 0.1). In contrast, the phase amplitude is stabilized around 0.4 under our optimal illumination, which is closer to the ground truth. Moreover, it can be seen that the phase detail in the green marked regions, which corresponds to the weak PTF frequency responses under the uniform and radial illumination patterns, cannot be recovered correctly. However, the optimal illumination scheme achieves an isotropic resolution, as well as much more accurate phase values, due to its isotropic PTF with much stronger responses.

4. EXPERIMENTAL RESULTS

4.1 A. Phase Resolution Target

To verify the effectiveness of the optimal illumination scheme experimentally, we first measured a pure phase resolution target [Quantitative Phase Microscopy Target ( ${QPT}^{T M}$ ), Benchmark Technologies Corporation, USA]. Our setup was based on a commercial inverted microscope (IX83, Olympus), in which the original condenser diaphragm is replaced by a high-contrast an amorphous silicon (a-Si) thin-film transistor LCD screen [4.3 inch (1 inch = 2.54 cm), pixel resolution $480 \times 272$ ] ^[47]. In our experiments, the built-in halogen white light source with a green interference filter (central wavelength $λ = 550 nm$ , 45 nm bandwidth) was used for illumination, and the LCD screen was used to modulate the illumination in asymmetrical manners. The images were captured by an objective lens with a magnification of $10 \times$ and an NA of 0.25 (Olympus PLAN 10X/0.25), and finally digitalized by a charge coupled device (CCD) camera with the pixel size of 3.75 μm (the imaging source DMK 23U445). Figure 4(a) shows a bright-field image captured when the LCD screen is transparent, which has very little intensity contrast. To illustrate the imaging resolution more clearly, the small regions near the center of the image (green-boxed areas) of the reconstructed phases under the half-circular uniform illumination and optimal illumination schemes are shown in Figs. 4(c) and 4(d). In the reconstruction phase of Fig. 4(c), two more suitable regularization parameters are used to suppress low-frequency and high-frequency errors, respectively. To demonstrate the isotropy of the optimal illumination scheme, the phase values at three uniformly distributed frequencies of Fig. 4(d) are extracted to plot the quantitative curves. As shown in Fig. 4(e), the optimal illumination achieves isotropic resolution at different frequencies. Furthermore, line profiles along a small circle of the same radius (blue and red circles) are extracted and illustrated in Fig. 4(f) to quantitatively compare the highest achievable resolution. These results are generally consistent with the theoretical prediction, as well as our simulation results. It can be seen that the phase details in the green areas, which correspond to the frequency components with weak responses in the non-isotropic PTF of the uniform illumination pattern, cannot be recovered correctly, indicating that half-circular uniform illumination cannot provide phase reconstruction with isotropic transverse resolution and SNR. Once again, the phase reconstruction result obtained by using the optimal illumination scheme demonstrates much better isotropy.

Phase reconstruction results of a phase resolution target QPTTM. (a) A bright-field image. (b) A zoom-in of the interest region of the bright-field image. (c) Phase reconstruction result under the half-circular uniform illumination pattern. (d) Phase reconstruction result under the optimal illumination scheme. (d) Phase values along three small circles evenly distributed from 0 to 2NAobj under the optimal illumination pattern. (e) Phase values along a small circle of the same radius in (c), (d) under the half-circular uniform illumination pattern and optimal illumination pattern.

Fig. 4. Phase reconstruction results of a phase resolution target ${QPT}^{TM}$ . (a) A bright-field image. (b) A zoom-in of the interest region of the bright-field image. (c) Phase reconstruction result under the half-circular uniform illumination pattern. (d) Phase reconstruction result under the optimal illumination scheme. (d) Phase values along three small circles evenly distributed from 0 to $2 {NA}_{obj}$ under the optimal illumination pattern. (e) Phase values along a small circle of the same radius in (c), (d) under the half-circular uniform illumination pattern and optimal illumination pattern.

下载图片查看所有图片

4.2 B. Unstained HeLa Cells

The high-resolution QPI capability of the proposed optimal illumination scheme provides unique possibilities for the label-free imaging of cell growth in culture, using repeated imaging of cultures to assess the progression toward confluence over designated periods of time. In Fig. 5, we show the quantitative phase images of the human cervical adenocarcinoma epithelial (HeLa) cell division process over the course of 5 h. The experimental setup generally followed the parameters of the previous experiment, except that an objective lens with $10 \times$ , 0.4 NA (Olympus UPlanSApo 10X/0.4) and an additional $1.25 \times$ camera adapter were used (effective magnification $12.5 \times$ ). A time-lapse movie created with one-phase reconstruction per 22 s is provided in Visualization 1. We show the full-field-of-view phase reconstruction result in this video [one frame from the video is shown in Fig. 5(a)]. From two selected zoom-in regions (red-boxed and blue-boxed areas) in Figs. 5(b) and 5(c), subcellular features, such as cytoplasmic vesicles and pseudopodium, can be clearly observed. In Fig. 5(d), we further selected one cell [corresponding to the green-boxed region shown in Fig. 5(a)] to study its morphology during division, which spanned over about 1 h. These high-resolution phase images clearly reveal the cell morphological changes during different mitosis phases. In addition, since the optimal illumination scheme requires only two-axis illuminations, all these retracting, extending, reorganizing, migrating, and maturing processes of the cell(s) were recovered accurately without any motion blur. These results demonstrate that the optimal illumination scheme is capable of imaging unlabeled cells in a non-invasive manner, allowing for a high-resolution QPI over an extended period of time.

Phase reconstruction results of HeLa cells under the optimal illumination scheme. (a) Full-field-of-view phase distribution. (b), (c) Phase maps of two selected zooms. (d) Phase results at different time points.

Fig. 5. Phase reconstruction results of HeLa cells under the optimal illumination scheme. (a) Full-field-of-view phase distribution. (b), (c) Phase maps of two selected zooms. (d) Phase results at different time points.

下载图片查看所有图片

5. DISCUSSION AND CONCLUSIONS

In this work, a rigorous theoretical analysis about the necessary and sufficient conditions for DPC to achieve an isotropic PTF has been explored. We have derived an optimal illumination scheme to maximize the frequency response and meanwhile achieve isotropic DPC. Compared with a traditional DPC method that uses half-circle illumination, our optimal illumination scheme produces a perfectly circularly symmetrical PTF with only two-axis intensity measurements, avoids missing frequencies, and enhances the phase response, providing high-quality phase reconstruction with isotropic transverse resolution and SNR. The resultant PTF removes the ill-posedness of the PTF inversion so that artifacts in phase reconstruction results can be significantly reduced. A theoretical analysis, simulations, and experimental results have verified the superiority of our method over existing illuminations in both the phase reconstruction accuracy and noise-robustness. The investigation of live HeLa cell mitosis in vitro has demonstrated that our optimal DPC scheme is a simple, efficient, and stable approach for label-free quantitative cell imaging with subcellular resolution. Furthermore, the intrinsic advantages, such as being non-interferometric, having the compatibility with bright-field microscopic hardware, and incoherent diffraction-limited resolution up to $2 {NA}_{obj}$ , make it a competitive and promising technique for various microscopy applications in life sciences and biophotonics.

In this work, we derived the optimal illumination scheme in an ideal imaging system. It should be noted that the actual imaging system is often more complicated due to the aberrations and noise. Several studies have improved the imaging performance of DPC from other aspects, such as the optimized reconstruction algorithm, regularization algorithm, and aberration compensation ^[38,39]. In this paper, to make a fair comparison, we adopted the most common and the simplest algorithm to ensure that the improvement of the phase quality can only result from the illumination optimization. In Appendix E, it is demonstrated that the optimal illumination can still achieve stable high-quality imaging results in an imaging system with a small aberration. In essence, the optimal illumination scheme can be combined with other optimized compensation approaches to further improve the phase reconstruction quality of DPC. On the other hand, compared with the existing 12-axis and three-axis method, the optimal scheme improves the imaging speed while ensuring the isotropic, but does not exceed the imaging speed of the traditional two-axis DPC. When imaging rapidly varying samples, the optimal illumination cannot completely avoid motion artifacts. Combining our approach with color-multiplexing techniques to achieve single-shot isotropic DPC imaging can potentially solve this problem, which is an interesting direction for future work.

References

[1] PopescuG., Quantitative Phase Imaging of Cells and Tissues (McGraw Hill Professional, 2011).

[2] A. Barty, K. Nugent, D. Paganin, A. Roberts. Quantitative optical phase microscopy. Opt. Lett., 1998, 23: 817-819.

[3] E. Cuche, F. Bevilacqua, C. Depeursinge. Digital holography for quantitative phase-contrast imaging. Opt. Lett., 1999, 24: 291-293.

[4] Y. Kim, H. Shim, K. Kim, H. Park, S. Jang, Y. Park. Profiling individual human red blood cells using common-path diffraction optical tomography. Sci. Rep., 2014, 4: 6659.

[5] G. Popescu. Quantitative phase imaging of nanoscale cell structure and dynamics. Methods Cell Biol., 2008, 90: 87-115.

[6] C. J. Mann, L. Yu, C.-M. Lo, M. K. Kim. High-resolution quantitative phase-contrast microscopy by digital holography. Opt. Express, 2005, 13: 8693-8698.

[7] P. Marquet, B. Rappaz, P. J. Magistretti, E. Cuche, Y. Emery, T. Colomb, C. Depeursinge. Digital holographic microscopy: a noninvasive contrast imaging technique allowing quantitative visualization of living cells with subwavelength axial accuracy. Opt. Lett., 2005, 30: 468-470.

[8] B. Kemper, G. von Bally. Digital holographic microscopy for live cell applications and technical inspection. Appl. Opt., 2008, 47: A52-A61.

[9] M. R. Teague. Deterministic phase retrieval: a Green’s function solution. J. Opt. Soc. Am., 1983, 73: 1434-1441.

[10] S. S. Kou, L. Waller, G. Barbastathis, C. J. Sheppard. Transport-of-intensity approach to differential interference contrast (TI-DIC) microscopy for quantitative phase imaging. Opt. Lett., 2010, 35: 447-449.

[11] J. C. Petruccelli, L. Tian, G. Barbastathis. The transport of intensity equation for optical path length recovery using partially coherent illumination. Opt. Express, 2013, 21: 14430-14441.

[12] C. Zuo, Q. Chen, W. Qu, A. Asundi. High-speed transport-of-intensity phase microscopy with an electrically tunable lens. Opt. Express, 2013, 21: 24060-24075.

[13] C. Zuo, Q. Chen, W. Qu, A. Asundi. Noninterferometric single-shot quantitative phase microscopy. Opt. Lett., 2013, 38: 3538-3541.

[14] C. Zuo, Q. Chen, Y. Yu, A. Asundi. Transport-of-intensity phase imaging using Savitzky–Golay differentiation filter-theory and applications. Opt. Express, 2013, 21: 5346-5362.

[15] F. Pfeiffer, T. Weitkamp, O. Bunk, C. David. Phase retrieval and differential phase-contrast imaging with low-brilliance X-ray sources. Nat. Phys., 2006, 2: 258-261.

[16] D. Hamilton, C. Sheppard. Differential phase contrast in scanning optical microscopy. J. Microsc., 1984, 133: 27-39.

[17] B. Kachar. Asymmetric illumination contrast: a method of image formation for video light microscopy. Science, 1985, 227: 766-768.

[18] S. B. Mehta, C. J. Sheppard. Quantitative phase-gradient imaging at high resolution with asymmetric illumination-based differential phase contrast. Opt. Lett., 2009, 34: 1924-1926.

[19] L. Tian, J. Wang, L. Waller. 3d differential phase-contrast microscopy with computational illumination using an LED array. Opt. Lett., 2014, 39: 1326-1329.

[20] L. Tian, L. Waller. Quantitative differential phase contrast imaging in an LED array microscope. Opt. Express, 2015, 23: 11394-11403.

[21] G. Zheng, R. Horstmeyer, C. Yang. Wide-field, high-resolution Fourier ptychographic microscopy. Nat. Photonics, 2013, 7: 739-745.

[22] X. Ou, R. Horstmeyer, C. Yang, G. Zheng. Quantitative phase imaging via Fourier ptychographic microscopy. Opt. Lett., 2013, 38: 4845-4848.

[23] L. Tian, X. Li, K. Ramchandran, L. Waller. Multiplexed coded illumination for Fourier ptychography with an LED array microscope. Biomed. Opt. Express, 2014, 5: 2376-2389.

[24] C. Zuo, J. Sun, Q. Chen. Adaptive step-size strategy for noise-robust Fourier ptychographic microscopy. Opt. Express, 2016, 24: 20724-20744.

[25] J. Sun, C. Zuo, L. Zhang, Q. Chen. Resolution-enhanced Fourier ptychographic microscopy based on high-numerical-aperture illuminations. Sci. Rep., 2017, 7: 1187.

[26] C. Zuo, J. Sun, J. Li, J. Zhang, A. Asundi, Q. Chen. High-resolution transport-of-intensity quantitative phase microscopy with annular illumination. Sci. Rep., 2017, 7: 7654.

[27] J. Li, Q. Chen, J. Sun, J. Zhang, X. Pan, C. Zuo. Optimal illumination pattern for transport-of-intensity quantitative phase microscopy. Opt. Express, 2018, 26: 27599-27614.

[28] J. Sun, C. Zuo, J. Zhang, Y. Fan, Q. Chen. High-speed Fourier ptychographic microscopy based on programmable annular illuminations. Sci. Rep., 2018, 8: 7669.

[29] J. Sun, Q. Chen, Y. Zhang, C. Zuo. Efficient positional misalignment correction method for Fourier ptychographic microscopy. Biomed. Opt. Express, 2016, 7: 1336-1350.

[30] D. Lee, S. Ryu, U. Kim, D. Jung, C. Joo. Color-coded LED microscopy for multi-contrast and quantitative phase-gradient imaging. Biomed. Opt. Express, 2015, 6: 4912-4922.

[31] C. Zuo, J. Sun, S. Feng, M. Zhang, Q. Chen. Programmable aperture microscopy: a computational method for multi-modal phase contrast and light field imaging. Opt. Lasers Eng., 2016, 80: 24-31.

[32] W. Lee, D. Jung, S. Ryu, C. Joo. Single-exposure quantitative phase imaging in color-coded LED microscopy. Opt. Express, 2017, 25: 8398-8411.

[33] Z. F. Phillips, M. Chen, L. Waller. Single-shot quantitative phase microscopy with color-multiplexed differential phase contrast (cDPC). PLoS ONE, 2017, 12: e0171228.

[34] M. Kellman, M. Chen, Z. F. Phillips, M. Lustig, L. Waller. Motion-resolved quantitative phase imaging. Biomed. Opt. Express, 2018, 9: 5456-5466.

[35] Y.-Z. Lin, K.-Y. Huang, Y. Luo. Quantitative differential phase contrast imaging at high resolution with radially asymmetric illumination. Opt. Lett., 2018, 43: 2973-2976.

[36] H.-H. Chen, Y.-Z. Lin, Y. Luo. Isotropic differential phase contrast microscopy for quantitative phase bio-imaging. J. Biophoton., 2018, 11: e201700364.

[37] A. Robey, V. Ganapati. Optimal physical preprocessing for example-based super-resolution. Opt. Express, 2018, 26: 31333-31350.

[38] KellmanM.BostanE.RepinaN.WallerL., “Physics-based learned design: optimized coded-illumination for quantitative phase imaging,” IEEE Trans. Comput. Imaging (to be published).

[39] M. Chen, Z. F. Phillips, L. Waller. Quantitative differential phase contrast (DPC) microscopy with computational aberration correction. Opt. Express, 2018, 26: 32888-32899.

[40] H. Rose. Nonstandard imaging methods in electron microscopy. Ultramicroscopy, 1976, 2: 251-267.

[41] D. Hamilton, C. Sheppard, T. Wilson. Improved imaging of phase gradients in scanning optical microscopy. J. Microsc., 1984, 135: 275-286.

[42] J. Sun, Q. Chen, J. Zhang, Y. Fan, C. Zuo. Single-shot quantitative phase microscopy based on color-multiplexed Fourier ptychography. Opt. Lett., 2018, 43: 3365-3368.

[43] C. J. Sheppard, S. Roth, R. Heintzmann, M. Castello, G. Vicidomini, R. Chen, X. Chen, A. Diaspro. Interpretation of the optical transfer function: significance for image scanning microscopy. Opt. Express, 2016, 24: 27280-27287.

[44] BerteroM.BoccacciP., Introduction to Inverse Problem in Imaging (CRC Press, 1998).

[45] Y. Fan, J. Sun, Q. Chen, J. Zhang, C. Zuo. Wide-field anti-aliased quantitative differential phase contrast microscopy. Opt. Express, 2018, 26: 25129-25146.

[46] R. Horstmeyer, R. Heintzmann, G. Popescu, L. Waller, C. Yang. Standardizing the resolution claims for coherent microscopy. Nat. Photonics, 2016, 10: 68-71.

[47] C. Zuo, J. Sun, S. Feng, Y. Hu, Q. Chen. Programmable colored illumination microscopy (PCIM): a practical and flexible optical staining approach for microscopic contrast enhancement. Opt. Lasers Eng., 2016, 78: 35-47.

1. INTRODUCTION

2. OPTIMAL ILLUMINATION SCHEME DESIGN

2.1 A. Derivation of Generalized Illumination Function for Isotropic Differential Phase Contrast

2.2 B. Optimal Illumination Scheme

3. IMAGING PERFORMANCE OF OPTIMAL ILLUMINATION SCHEME

4. EXPERIMENTAL RESULTS

4.1 A. Phase Resolution Target

4.2 B. Unstained HeLa Cells

5. DISCUSSION AND CONCLUSIONS

Yao Fan, Jiasong Sun, Qian Chen, Xiangpeng Pan, Lei Tian, Chao Zuo. Optimal illumination scheme for isotropic quantitative differential phase contrast microscopy[J]. Photonics Research, 2019, 7(8): 08000890.

Optimal illumination scheme for isotropic quantitative differential phase contrast microscopy Download： 720次

1. INTRODUCTION

2. OPTIMAL ILLUMINATION SCHEME DESIGN

2.1 A. Derivation of Generalized Illumination Function for Isotropic Differential Phase Contrast

Fig. 1. Schematic diagram of the integral for PTF along the left–right axis in the polar coordinate system. (a) The radius $ρ$ of the point $Q$ is in the range of ${NA}_{obj} \leq ρ \leq 2 {NA}_{obj}$ . (b) The radius $ρ$ of the point $Q$ is in the range of $0 \leq ρ < {NA}_{obj}$ .

2.2 B. Optimal Illumination Scheme

3. IMAGING PERFORMANCE OF OPTIMAL ILLUMINATION SCHEME

4. EXPERIMENTAL RESULTS

4.1 A. Phase Resolution Target

4.2 B. Unstained HeLa Cells

Fig. 5. Phase reconstruction results of HeLa cells under the optimal illumination scheme. (a) Full-field-of-view phase distribution. (b), (c) Phase maps of two selected zooms. (d) Phase results at different time points.

5. DISCUSSION AND CONCLUSIONS

Article Outline

关于本站 Cookie 的使用提示

全站搜索

Optimal illumination scheme for isotropic quantitative differential phase contrast microscopy Download： 720次

1. INTRODUCTION

2. OPTIMAL ILLUMINATION SCHEME DESIGN

2.1 A. Derivation of Generalized Illumination Function for Isotropic Differential Phase Contrast

Fig. 1. Schematic diagram of the integral for PTF along the left–right axis in the polar coordinate system. (a) The radius ρ of the point Q is in the range of NAobj≤ρ≤2NAobj. (b) The radius ρ of the point Q is in the range of 0≤ρ<NAobj.

2.2 B. Optimal Illumination Scheme

3. IMAGING PERFORMANCE OF OPTIMAL ILLUMINATION SCHEME

4. EXPERIMENTAL RESULTS

4.1 A. Phase Resolution Target

4.2 B. Unstained HeLa Cells

Fig. 5. Phase reconstruction results of HeLa cells under the optimal illumination scheme. (a) Full-field-of-view phase distribution. (b), (c) Phase maps of two selected zooms. (d) Phase results at different time points.

5. DISCUSSION AND CONCLUSIONS

Article Outline

相关论文

相关资讯

关于本站 Cookie 的使用提示

全站搜索

Fig. 1. Schematic diagram of the integral for PTF along the left–right axis in the polar coordinate system. (a) The radius $ρ$ of the point $Q$ is in the range of ${NA}_{obj} \leq ρ \leq 2 {NA}_{obj}$ . (b) The radius $ρ$ of the point $Q$ is in the range of $0 \leq ρ < {NA}_{obj}$ .