## Abstract

We demonstrate a new approach to the transport of intensity equation (TIE) phase retrieval method which uses structured illumination to improve low-frequency noise performance. In the hybrid scheme, two phase images are acquired: one with uniform and one with sinusoidal grating illumination intensity. The former preserves the high spatial frequency features of the phase best, whereas the latter dramatically increase the response at low spatial frequencies (where traditional TIE notoriously suffers). We then theoretically prove the design of a spectral filter that optimally combines the two phase results while suppressing noise. The combination of uniformly and structured illuminated TIE (hybrid TIE) phase imaging is experimentally demonstrated optically with a calibrated pure phase object.

© 2014 Optical Society of America

## 1. Introduction

Phase information cannot be obtained directly from electromagnetic waves in the frequency range from THz and above; yet this information is often extremely valuable, as it records the strongest evidence of light interaction with materials of interest. Traditional imaging systems record the intensity of light, which maps to the imaginary part of the complex refractive index *n*; however, phase provides better contrast for many important classes of objects such as biological cells at visible frequencies, low–Z materials (including liquids and soft tissue at x–ray frequencies, etc.) [1]. The recovery of the real part of *n* is commonly referred to as “Quantitative Phase Retrieval,” even though the term is slightly misleading—it is not the phase itself that is of interest, but rather the optical path length through the specimen [2]. The latter interpretation has the benefit that it extends the problem naturally to partially coherent systems, where phase is not well defined [3]. Numerous techniques have been developed for the phase retrieval problem, including phase–shifting interferometry [4–6], intensity–based iterative [7, 8] or transport [9] methods, or methods based on Ptychography [10–12].

Among different quantitative phase retrieval methods, the transport of intensity equation (TIE) is convenient because it allows for a simple experimental setup with low demands on both the number and quality of the required optical elements [9, 13]. This feature makes TIE–based methods especially appealing for the X–ray regime; moreover, its demonstrated good tolerance to source incoherence [13–15] also relieves the necessity of using a synchrotron, opening up the possibility of eventual use with tabletop X-ray sources [16].

In the TIE method, the phase of the object is extracted from measurements of the intensity differentials along the propagation axis that occur due to propagation. For pure phase objects, these differentials act as sources and sinks in a Poisson equation which can then be inverted to the phase. However, because the Poisson inversion process significantly amplifies low spatial frequency information, the TIE has an intrinsic problem of amplifying low–frequency noise, which results in cloudy phase results [17,18]. Intuitively, this can be understood because small-scale phase features diffract light significantly, whereas large smooth objects bend the light only slightly, resulting in minimal intensity variations along the propagation direction. The problem is especially destructive in the case of smooth phase objects (e.g. cells, tissue) or low Fresnel numbers (small propagation distances). Many techniques have been proposed to tackle the problem of low frequency noise corruption in TIE, the most popular being computational regularization, linear (Tikhonov) [19] or nonlinear. The latter was recently demonstrated also with sparsity priors, specifically sparsity of hard edges in an otherwise smooth object, to fix the low frequencies [20]. These solutions effectively remove low frequencies, so produce nice results only when low frequency phase information does not exist or is not important. To properly recover low–frequency information, one must change the optical system in order to better capture such information. One way of doing this is to use a large defocus distance in the intensity image, which when combined with small defocus images can produce images with good fidelity across all spatial frequencies [17,18,21–23]. However, large propagation distances may be physically impossible or impractical. Another method places a diffuser in the sample plane in order to amplify the diffraction effects on propagation and thereby achieve better low frequency contrast [24]. Diffraction effects can also be enhanced by structuring the amplitude of the illumination with a spatial light modulator [25,26] and recovering the first gradient of phase for multiple directions, similar to methods in Xray grating interferometry [6, 27]. For example, lateral shear interferomery has been used for phase imaging pure low frequency objects whose high frequency diffraction can be ignored [28, 29]. Recently, Pang *et al* [30] achieved phase imaging using random illumination patterns with 10 intensity measurements under linearization assumptions.

In this paper, we introduce a novel hybrid TIE quantitative phase retrieval method with optimal response for pure phase objects at all spatial frequency bands. Our phase imaging method only requires three intensity measurements, one with uniform illumination and two with orthogonal sinusoidal (structured) illumination, all at the same defocused plane. By subtracting the diffraction signal with the uniform illumination from that with the structured illumination, we reveal a clean modulation induced signal that is solely depend on the phase gradient. This result is rigorously derived from the original TIE equation without any further approximation or linearization.

Our approach to to low-frequency response enhancement is inspired by structured illumination microscopy [31]. Gustafsson’s insight was that modulating the object intensity shifts the high frequency components of the object down to low frequency region, enabling imaging resolution beyond the detector limit. Essentially, therefore, structured illumination microscopy solves the opposite problem of ours, making up for missing high spatial frequencies to achieve super resolution. Here, we propose to use the same principle, but in reverse: namely, to shift the *low* frequencies higher, thus making up for TIE’s inherent noise corruption problem. In that sense, it is also analogous to the HiLo reconstruction method [32], which has been shown to be effective in removing low frequency de-focused haze in tissue fluorescence microscopy. Even though in our case the origin of the low-frequency noise is different, we observe similar improvements in the hybrid reconstruction. We also show how to combine and optimize the phase retrieval results from the uniform illuminated TIE and the phase gradient signal to optimally mitigate the effects of the noise at all spatial frequency bands.

The remainder of the paper is structured as follows: in Section 2 we derive the forward and inverse operators for TIE with structured illumination. Based on these, in Section 3 we then design an optimal spectral filter that ideally selects and mixes signal from the uniformly and structure-illuminated TIE images, resulting in a hybrid quantitative phase image with minimum noise power spectrum over the entire spectral domain. Numerical simulation and experimental results with visible light are presented in Section 4.

## 2. Uniformly and structure-illuminated TIE phase imaging

Let *ϕ*(*x*, *y*) denote the phase of a purely transmitting object located at *z* = 0, and let *I*(*x*, *y; d*) denote the intensity after propagation by a (small) distance *d* behind the object or its conjugate in a suitable imaging system. When explicit reference to the lateral coordinates can be omitted without risk of confusion, we will abbreviate as *I*(*d*).

The TIE method exploits Fresnel propagation to relate the phase *ϕ*(*x*, *y*) with the intensity distribution differential *∂I*(*x*, *y*; *z*)/*∂z* along the propagation axis *z*. According to the TIE [9],

*x*and

*y*are coordinates in the plane perpendicular to the propagation

*z*axis,

*I*and

*ϕ*are the intensity and phase profiles at

*z*=

*z*

_{0}, respectively,

*k*is the wave number of the optical field, and ∇

_{⊥}≡ (

**i**

*∂/∂x*+

**j**

*∂/∂y*), is the gradient operator on the

*x*−

*y*plane, with

**i**and

**j**being the unit direction vector in the

*x*and

*y*directions, respectively.

The simplest case is that of a pure phase object and uniform illumination, which is well known but we repeat here for convenience and to establish notation. Let *I*_{0} denote the incident intensity, as shown in Fig. 1(a). The intensity distribution is assumed unchanged (still *I*_{0}) right after the thin-phase object. For a pure phase object, Eq. (1) yields

*d*) is the intensity differential signal. This is a Poisson equation and it can be solved in the Fourier domain as

*ℱ*is the Fourier transform operator,

*u*and

*v*are the spatial frequency coordinates, and is the forward transfer function of the Poisson operator. From this result, it is evident that the direct reconstruction described by Eq. (3) is ill-posed at low frequencies and singular at DC—the latter not surprising, since a constant phase factor is usually physically meaningless. Typically, a small regularization parameter is used to avoid the division-by-zero instability, but the spectrum of the intensity differential signal Δ(

*d*) and the noise it carries are still amplified by the inversion process, with an inverse quadratic dependence on spatial frequency. As a result, Eq. (3) typically produces phase estimates that are heavily corrupted by low frequency noise. Figure 2 illustrates this point by plotting the transfer function

*H*and the inverse transfer function

*H*

^{−1}as a function of frequency

*u*.

Here, we show that these low–frequency noise problems can be significantly improved by our hybrid structured-illumination TIE method. The idea is that low–frequencies can be shifted to high frequencies through the use of structured illumination. The procedure is shown schematically in Fig. 3. The pure phase object *ϕ*(*x*, *y*) is illuminated with a patterned intensity distribution *I*_{m}(*x*, *y*, *d*) as shown in Fig. 1(b). In the simplest case that we investigate here, the modulation is sinusoidal along an arbitrary lateral direction, e.g. *x*, as

*a*is the intensity gain/attenuation,

*η*is the modulation depth, and

*f*

_{m}is a suitably selected spatial frequency. (More on this later.)

After Fresnel propagation in the modulated illumination case, the intensity differential becomes

*I*(

_{m}*d*) is the intensity measurement after propagation by distance

*d*. Compared to Eq. (2), an additional cross term ∇

*I*

_{m}· ∇

*ϕ*contributes to the new intensity differential signal Δ

_{m}(

*d*). We can isolate the cross term by combining Eqs. (2) and (6), yielding

*s*

_{m}(

*d*) as the “modulation cross–term” signal. Note that no movement of the camera is necessary:

*I*

_{0}and

*I*

_{m}(0) are known (or can be calibrated with a measurement absent any object) while the two measurements

*I*(

*d*) and

*I*

_{m}(

*d*) are obtained with uniform and structured illumination, respectively, downstream at the (out–of–focus) location

*d*.

A derivation similar to Eqs. (5)–(7) has been reported for TIE under general coded illumination [25]; here, we proceed to exploit the equivalence of this result to a convolution (denoted as ⊗) in the spatial frequency domain. Taking the Fourier transform on both sides of this equation, we obtain

Equation (8) is easily solved when *I*_{m} is a sinusoidal pattern, as in Eq. (5). The result is

*s*

_{m,X}(

*d*) is the modulation cross–term signal after propagation by distance

*d*obtained by using the sinusoidal modulation pattern

*I*

_{m,X}(0). Note that the signal is proportional to the modulation frequency

*f*

_{m}and depth

*η*. The spectrum of the phase differential signal

*ϕ*is shifted back to the origin by convolution with the Dirac delta function

_{x}*δ*(

*u*+

*f*

_{m}) as

We can eliminate the double-modulated term *ℱ*(*ϕ _{x}*)(

*u*+ 2

*f*

_{m},

*v*) in two ways. A simple approach is to use a low pass filter ${T}_{x}(u,v)=\text{exp}\left[-{\left(u+{f}_{\text{m}}\right)}^{2}/{f}_{\text{m}}^{2}\right]$, which is centered at

*u*= 0 and cut off around

*u*=

*f*

_{m}. This approach is valid if the spectrum of the phase differential has negligible components beyond the modulation frequency

*f*

_{m}, in which case the two lobes of shifted spectrum do not overlap. The phase differential can then be obtained by

*linear*in the spatial frequency variable where the modulation took place, whereas in the uniform–illuminated case it was quadratic in both spatial frequency variables (compare Eq. (4) and Fig. 2). This explains why the noise influence on low spatial frequencies is expected to be milder in the structured–illumination case, as demonstrated in Fig. 2. The figure also shows that the quadratic response curve eventually catch up and surpass the linear response in the high frequency region beyond some critical frequency

*f*.

_{c}When the object contains substantial high frequency components beyond *f*_{m}, additional measurements can be used to eliminate the double-modulated term *ℱ*(*ϕ _{x}*)(

*u*+ 2

*f*

_{m}) [33]. For example, a second modulation cross–term measurement can be obtained by shifting the illumination by

*π*/2

*s*

_{q,X}(

*d*) relates to the spectrum of the phase differential signal as

The direct inversion described in principle in the two above cases (low–pass filtering and quadrature) is problematic, as we will also see in section 4. In practice, we include a second set of measurements using modulation in the perpendicular direction *y*. Thus, for the quadrature case, we obtain independently the phase differential in the *y* direction as

*s*

_{m,Y}(

*d*) and

*s*

_{q,Y}(

*d*) are modulation cross–term signals with intensity modulation patterns and respectively. (The low–pass filter approach can be similarly derived and we omit here.)

Having both *ϕ _{x}* and

*ϕ*independently available, the most robust way [34, 35] to recover

_{y}*ϕ*is by re–constructing the Poisson equation

*u*,

*v*. This is beneficial, but further improvement may be obtained by combining the estimate of

*ϕ*obtained from Eq. (23) with the standard Poisson solution (3). This is done in the next section.

## 3. Error statistics and optimized hybrid TIE reconstruction

Returning to Fig. 2, as we pointed out in the previous section, the inverse transfer function behaves better in the linear case (structured illumination) for spatial frequencies lower than the critical frequency *f _{c}* (Eq. (18)) whereas above

*f*the reconstruction from uniform illumination should be better. To quantitively describe the “goodness” of imaging strategies, we introduce the power spectral density (PSD) error as a statistical measurement of phase noise corruption. Based on this scheme, we show how to optimally combine the two reconstructions under the assumption of additive white Gaussian noise. The result is our final hybrid TIE reconstruction filter.

_{c}To carry out the derivation, we begin by evaluating the noise statistics for four reconstruction strategies, namely:

- uniform illumination alone (standard TIE), which we will denote as
*ϕ̂*_{1}; - sinusoidal illumination along one direction with quadrature recovery, denoted as
*ϕ̂*_{2}; - sinusoidal illumination alone with two orthogonal direction measurements and quadrature recovery, denoted as
*ϕ̂*_{3}; - combined uniform and structured illumination, the latter being sinusoidal with two orthogonal direction measurements and quadrature recovery; the exact optimal way of doing the combination will emerge at the end of this section from the preceding analysis, and we will denote the result as
*ϕ̂*_{h}.

*ϕ̂*

_{1},

*ϕ̂*

_{3}will be used in the optimal filter derivation. Other reconstruction strategies (e.g., low-pass filter recovery) can also be carried out similarly, but we have omitted them for brevity and because they are sub–optimal in terms of noise performance.

Let the noise–corrupted measurement of the differential signal Δ̂(*d*) be denoted as

*n*(

*x*,

*y*) is a wide–sense stationary random process of zero mean. This measurement noise originates from the optical intensity measurement, normally consisting of dark noise of the detector and shot noise from photon statistics. The detector dark noise is uncorrelated for each CCD pixel, and is not dependent on the measured intensity. As a result, the dark noise has a flat spatial power spectrum within the bandwidth of the CCD. The shot noise power scales with the intensity and its spatial spectrum is related to the input intensity profile. However, in typical TIE experiments, the intensity differential Δ(

*d*) and Δ

_{m}(

*d*) are small compared to the initial intensity

*I*

_{0}. For uniform and weakly modulated input fields, the shot noise level of each pixel is also similar to others, resulting in spatially white noise. Under all the above assumptions, the power spectral density of the intensity differential noise

*n*is constant and we define it as where

*σ*

^{2}is the aggregated noise variance.

Starting with the uniformly illuminated case, let *ϕ̂*_{1} denote the noisy estimate of the phase that we obtain from the noisy measurement Eq. (24). Using Eq. (3) and the Wiener–Khintchine theorem, it is straightforward to show that the power–spectral densities *S _{ϕ}*(

*u*,

*v*),

*S*

_{ϕ̂1}(

*u*,

*v*) are related as

*S*(

_{n}*u*,

*v*) is the PSD of

*n*(

*x*,

*y*), and

*E*(

*u*,

*v*) is the PSD error we wish to minimize. Substituting the noise PSD from Eq. (25), the well–known PSD error of the phase retrieval is calculated for uniform TIE in our notation as

Similarly, we can calculate the noisy phase reconstruction and extract the PSD error of the phase for the sinusoidal-patterned-illuminated TIE. In this case, the PSD error of the phase depends on both the way one recovers phase differentials from Eq. (10) (low–pass filtering or quadrature) and the approach one uses for the final retrieval (direct integral or Poisson reconstruction).

We first consider the PSD error *E*_{ϕ̂x} of the phase differential *ϕ _{x}* itself. When

*ϕ*is obtained by the quadrature method as in Eq. (16), the PSD error becomes

_{x}*d*) and Δ̂

_{m}(

*d*) are identical. The low–pass filter recovery case yields a similar result, but larger by a factor of 2 (compare Eqs. (13) and (17)).

We next consider the PSD error when relying on a single sinusoidal structured illumination along the *x* direction, as in Eq. (12), and we denote the reconstruction result as *ϕ̂*_{2}. Then

*ϕ̂*

_{3}as

*ϕ*has the same statistics as that of

_{y}*ϕ*, as is true for systems with orientation symmetry.

_{x}Inserting the expression of PSD error of phase differentials from Eq. (28) into the equations above, we finally obtain the corresponding PSD errors of two approaches as

The power spectral densities *E*_{ϕ̂1} (uniform-illuminated TIE), *E*_{ϕ̂2} (structured illuminated TIE, forward integral) and *E*_{ϕ̂2} (structured illuminated TIE, Poisson reconstruction) are shown in Fig. 4 in log scale. We can see, firstly, the clear improvement of noise suppression using two orthogonal sinusoidal illumination steps as in Eq. (23) compared to one-directional integration. This is because the one-directional integration yields an anisotropic noise spectrum—inversely quadratic along the *u* direction only, and constant along the *v* direction. As a result, *ϕ̂*_{2} is heavily corrupted by noise on the *v* axis. The isotropic power spectral density error of *E*_{ϕ̂3}, on the other hand, is inversely linear in all directions. In fact, *ϕ*_{2} is the least-square phase retrieval solution given phase differential profiles in both directions. These conclusions are further demonstrated by the numerical simulations further along in Fig. 5.

On the other hand, the isotropic profiles of PSD error spectra *E*_{ϕ̂1} and *E*_{ϕ̂3} are consistent with the inverse transfer function behavior depicted in Fig. 2. As predicted in the previous section, noise is relatively suppressed in the low frequency region for the structured illumination TIE but decays relatedly faster in the high frequency region for the uniform illuminated TIE. To obtain the best reconstruction independent of spatial frequency, an optimized phase estimate can be constructed by weighting each of the two signals in the hybrid set of measurements according to its frequency response in different spectral regions. In other words, we can combine high frequency components from the uniformly-illuminated TIE and the low frequency components from the structured-illumination TIE to minimize the overall effect of the noise.

This is formally done as follows: first, we form a convex linear combination of the two reconstructions *ϕ̂*_{1} and *ϕ̂*_{3}, filtered respectively as

*t*(

*x*,

*y*) is the impulse response of the shift–invariant filter that we wish to optimize. We assume that

*t*is square integrable and denote its Fourier transform (transfer function) as

*T*(

*u*,

*v*). Forming the autocorrelation of the random process

*ϕ̂*

_{h}and taking its Fourier transform, it is straightforward to show that the PSD error for the new estimate

*ϕ̂*

_{h}is

*u*,

*v*) = (0, 0). The latter exclusion is justified because the DC component of a phase signal has no physical meaning and so we are justified in ignoring it. Under these assumptions, the total noise power

*P*can be computed as

_{N}_{★}denotes the integral over all spatial frequencies excluding an appropriately small region around the origin, as noted just before.

To minimize *P _{n}*, we use the variational method. By setting the variation of

*P*to equal zero,

_{n}*f*, and we refer to it as the hybrid TIE (hybrid uniformly and structured illuminated TIE) filter. It yields the optimal (in the sense of minimal overall noise power) PSD error

_{c}The minimized noise spectrum is shown Fig. 4(d). We can see that, indeed, the overall noise spectrum is reduced compared to all the power spectra *E*_{ϕ̂1}, *E*_{ϕ̂2} and *E*_{ϕ̂3}.

## 4. Simulation and experimental results

Here we report simulation and experimental results to confirm and characterize the performance of various structured illumination approaches for TIE phase retrieval. We first simulate the phase imaging in the x-ray region. A pure phase profile of a 2×2 mm^{2} Lena’s image, as shown in Fig. 5(a), is used as the object. The largest phase value is set to ±2.6 rad, while the background is 0 rad. We assume that the illumination wavelength is in the hard x-ray range of 0.0633 nm (∼20 keV), propagation distance *d* = 0.4 m and the sinusoidally patterned illumination has a single spatial frequency *f*_{m} at 14.2 mm^{−1} with 50% modulation depth. Pixel size is assumed to be 4.4 *μ*m and space–bandwidth product is set to 1024×1024. Additive white noise in the intensity differential measurement is assumed to be *σ*^{2} = 10^{−4}, normalized to background intensity power
${I}_{0}^{2}=1$. All these numbers are within practical range of realistic x-ray experiment. Note that the propagation distance of 0.4 m is far smaller than the Talbot distance of the grid, thereby is insufficient for obtaining diffraction signal from low frequency components (< *f _{m}*) using uniform illumination TIE. However, multiplication of the propagation distance quickly becomes impractical for in-door operations. We show that structure illumination can revive low frequency signals under propagation distance constrain.

Figure 5(b) shows the phase reconstruction using uniformly-illuminated TIE, corrupted with low-frequency noise. The structured illuminated TIE phase retrieval *ϕ̂*_{2} shown in Fig. 5(c) contains substantial strip-like high–frequency noise along the *x* direction, as predicted (see Fig. 4(b).) The low frequency background noise is substantially suppressed for *ϕ̂*_{3} as shown in Fig. 5(d), where we used the modulation in both orthogonal directions. However, the image is not as sharp as in *ϕ̂*_{1}. The hybrid phase reconstruction *ϕ̂*_{h} shown in Fig. 5(e) both removes noise effectively and preserves sharp and smooth features alike, resulting in a high fidelity phase retrieval. In Figs. 5(f)–5(i) we also show the corresponding PSD error, as computed from the reconstructions, showing good agreement with the theoretical predictions given in Fig. 4.

We next built a TIE phase imaging experiment in the optical range. The experiment is designed to mimic the x-ray simulation presented above. We use the operation wavelength of 633 nm and a propagation distance of 0.4 mm, while keeping the pixel size at 4.4 *μ*m. By keeping the Fresnel number (*dx*)^{2}/*λd* fixed, the TIE signal stays the same. We use a lens (N-BK7 Half-Ball Lens 6.0 mm Diameter from Edmund Optics) as the phase object. To reduce the total phase difference thereby improving the condition of the inverse problem, the lens object was immersed in index-matching oil. At the operation wavelength of 633 nm, the refractive index difference between the lens and the immersion oil is effectively reduced to 0.0036±0.0001, yielding a nominal maximum phase variation of 106.6 rad. We chose the lens object deliberately because *(i)* its shape is calibrated so it provides ground truth and *(ii)* it contains mostly low spatial frequency features, thus it accentuates TIE’s problems at that band.

The experimental arrangement is shown in Fig. 6. A multimode fiber delivers the light beam from a broadband LED illuminator to a focusing lens. The beam is focused down to a pinhole (diameter 50 *μ*m) and then collimated to provide quasi–plane wave illumination. A narrow band (1 nm @ 633 nm) spectral filter is inserted in the collimator. In this way, we obtain uniform illumination with high spatial and temporal coherence. The object is imaged with a magnification ×0.5 via a 4f telescope system onto the detector. The lenses used in the telescope system are 2-inch aspheric lenses with focal length of 20 cm and 10 cm, respectively. The detector has pixel size of 4.4 *μ*m. It is placed on a rail that can be translated along the optical propagation axis *z*. In the uniformly-illuminated TIE experiment, we move the detector to obtain images both in focus and at a defocus distance *d* = 0.1 mm, corresponding to a propagation distance of *d* = 0.4 mm for the object. The additional in-focal imaging is taken for background normalization. For the structured illuminated TIE, we only take the defocused images with the patterned illumination, which is generated by inserting a sinusoidal amplitude grating (spatial frequency 10 mm^{−1}, modulation depth of 40%, Applied Imaging Inc.) immediately before the lens, but within the depth of focus such that it is effectively collocated with the sample. Both horizontal and vertical structured illumination were implemented in the experiment and one defocused image is taken for each. The phase differential is obtained via the filter in Eq. (11) and the phase profile is retrieved via the Poisson Eq. (23). Note that the 4-f system is only used to separate the object and image plane while keeping the effective propagation distance to 0.4 mm. No operation on the optical Fourier plane is implemented, as no 4-f system is available in the x-ray region.

Figure 7 shows the experimental phase retrieval results. As expected due to the low frequency content, the uniformly-illuminated TIE result (Fig. 7(a)) is gravely corrupted by noise. On the other hand, the structured–illumination approach on the other hand, successfully captures the phase profile (Fig. 7(b)), but misses the sharp boundary. Finally, the hybrid reconstruction *ϕ*_{h} shown in Fig. 7(c), also obtains the low–frequency phase profile successfully. Moreover, from Fig. 7(c) and the cross–sections shown in Fig. 7(d), we observe that a small bubble in the immersion oil was accidentally trapped in the vicinity of the lens. This bubble is not visible in Fig. 7(c) where, again as expected, high frequencies/small features are filtered out. The overall measured SNR (bubbles ignored) for the uniform, structure illumination TIE and the hybrid reconstruction are 0.79, 4.51 and 4.77, respectively. It is noted that the noise spectrum for the uniform TIE result shown in Fig. 7(a) is heavily concentrated in the low frequency domain. This could be the result of intensity noise due to propagation of illumination non-uniformities. Further optimization of the hybrid reconstruction can be made based on the noise spectrum obtained in the experiment.

## 5. Conclusion

We have demonstrated both theoretically and experimentally that substantial noise suppression is obtained, especially in the low–frequency region, by a hybrid reconstruction of TIE phase retrieval using both uniform and structured illumination. Future work will focus on designing the illumination pattern optimally for a given situation [25], as well as exploring the possibility of also achieving enhanced resolution through proper design on of the pattern. The latter has been done recently for phase imaging with interferometry [36] and can be incorporated into our algorithm in the future. It is noted that so far we have not utilized sparsity or other compressive priors. These should be expected to further suppress noise and improve reconstruction quality.

This research was funded by the US Department of Homeland Security. We would like to thank Zhengyun Zhang for helpful discussions.

## References and links

**1. **P. Bach, J. Jett, U. Pastorino, M. Tockman, S. Swensen, and C. Begg, “Computed tomography screening and lung cancer outcomes,” Jama **297**(9), 953–961 (2007). [CrossRef] [PubMed]

**2. **K. A. Nugent, D. Paganin, and T. E. Gureyev, “A phase odyssey,” Phys. Today **54**(8), 27–32 (2001). [CrossRef]

**3. **J. W. Goodman, *Statistical Optics* (Wiley-Interscience, 1985).

**4. **R. K. Crane, “Interference phase measurement,” Appl. Opt. **8**, 538 (1969).

**5. **O. Y. Kwon, “Multichannel phase-shifted interferometer,” Opt. Lett. **9**(2), 59–61 (1984). [CrossRef] [PubMed]

**6. **T. Weitkamp, A. Diaz, C. David, F. Pfeiffer, M. Stampanoni, P. Cloetens, and E. Ziegler, “X-ray phase imaging with a grating interferometer,” Opt. Express **13**(16), 6296–6304 (2005). [CrossRef] [PubMed]

**7. **R. Gerchberg and W. Saxton, “Phase determination for image and diffraction plane pictures in the electron microscope,” Optik **34**(3), 275–284 (1971).

**8. **J. R. Fienup, “Reconstruction of an object from the modulus of its Fourier transform,” Opt. Lett. **3**(1), 27–29 (1978). [CrossRef] [PubMed]

**9. **M. R. Teague, “Deterministic phase retrieval: a Green’s function solution,” J. Opt. Soc. Am. **73**(11), 1434–1441 (1983). [CrossRef]

**10. **J. M. Rodenburg and H. M. Faulkner, “A phase retrieval algorithm for shifting illumination,” Appl. Phys. Lett. **85**(20), 4795–4797 (2004). [CrossRef]

**11. **P. Thibault, M. Dierolf, O. Bunk, A. Menzel, and F. Pfeiffer, “Probe retrieval in ptychographic coherent diffractive imaging,” Ultramicroscopy **109**(4), 338–343 (2009). [CrossRef] [PubMed]

**12. **G. Zheng, R. Horstmeyer, and C. Yang, “Wide-field, high-resolution Fourier ptychographic microscopy,” Nat. Photon. **7**(9), 739–745 (2013). [CrossRef]

**13. **D. Paganin and K. A. Nugent, “Noninterferometric phase imaging with partially coherent light,” Phys. Rev. Lett. **80**(12), 2586 (1998). [CrossRef]

**14. **T. E. Gureyev, A. Roberts, and K. A. Nugent, “Partially coherent fields, the transport-of-intensity equation, and phase uniqueness,” J. Opt. Soc. Am. **12**(9), 1942–1946 (1995). [CrossRef]

**15. **T. E. Gureyev, Y. I. Nesterets, D. M. Paganin, A. Pogany, and S. W. Wilkins, “Linear algorithms for phase retrieval in the Fresnel region. 2. Partially coherent illumination,” Opt. Commun. **259**(2), 569–580 (2006). [CrossRef]

**16. **R. Lewis, “Medical phase contrast x-ray imaging: current status and future prospects,” Phys. Med. Biol. **49**(16), 3573 (2004). [CrossRef] [PubMed]

**17. **L. Waller, L. Tian, and G. Barbastathis, “Transport of intensity phase-amplitude imaging with higher order intensity derivatives,” Opt. Express **18**(12), 12552–12561 (2010). [CrossRef] [PubMed]

**18. **D. Paganin, A. Barty, P. J. McMahon, and K. A. Nugent, “Quantitative phase–amplitude microscopy. III. The effects of noise,” J. Microsc. **214**(1), 51–61 (2004). [CrossRef] [PubMed]

**19. **L. Turner, B. Dhal, J. Hayes, A. Mancuso, K. Nugent, D. Paterson, R. Scholten, C. Tran, and A. Peele, “X-ray phase imaging: Demonstration of extended conditions for homogeneous objects,” Opt. Express **12**(13), 2960–2965 (2004). [CrossRef] [PubMed]

**20. **L. Tian, J. C. Petruccelli, and G. Barbastathis, “Nonlinear diffusion regularization for transport of intensity phase imaging,” Opt. Lett. **37**(19), 4131–4133 (2012). [CrossRef] [PubMed]

**21. **E. D. Barone-Nugent, A. N. T. O. N. Barty, and K. A. Nugent, “Quantitative phase–amplitude microscopy I: optical microscopy,” J. Microsc. **206**(3), 194–203 (2002). [CrossRef]

**22. **C. Zuo, Q. Chen, Y. Yu, and A. Asundi, “Transport-of-intensity phase imaging using Savitzky-Golay differentiation filter-theory and applications,” Opt. Express **21**, 5346–5362 (2013). [CrossRef] [PubMed]

**23. **Z. Jingshan, R. A. Claus, J. Dauwels, L. Tian, and L. Waller, “Transport of Intensity phase imaging by intensity spectrum fitting of exponentially spaced defocus planes,” Opt. Express **22**(9), 10661–10674 (2014) [CrossRef] [PubMed]

**24. **P. F. Almoro, L. Waller, M. Agour, C. Falldorf, G. Pedrini, W. Osten, and S. G. Hanson, “Enhanced deterministic phase retrieval using a partially developed speckle field,” Opt. Lett. **37**(11), 2088–2090 (2012). [CrossRef] [PubMed]

**25. **A. Shanker, L. Tian, and L. Waller, “Defocus-based quantitative phase imaging by coded illumination,” Proc. SPIE **8949**, 89490R (2014).

**26. **Y. Zhu, A. Pan, and G. Barbastathis, “Low-noise TIE phase imaging by structured illumination,” in *Computational Optical Sensing and Imaging*, Optical Society of America (2014), CTh3C-5. [CrossRef]

**27. **F. Pfeiffer, T. Weitkamp, O. Bunk, and C. David, “Phase retrieval and differential phase-contrast imaging with low-brilliance X-ray sources,” Nat. Phys. **2**(4), 258–261 (2006). [CrossRef]

**28. **K. Ichikawa, A. W. Lohmann, and M. Takeda, “Phase retrieval based on the irradiance transport equation and the Fourier transform method: experiments,” Appl. Opt. **27**, 3433–3436 (1988). [CrossRef] [PubMed]

**29. **F. Roddier, “Wavefront sensing and the irradiance transport equation,” Appl. Opt. **29**, 1402–1403 (1990). [CrossRef] [PubMed]

**30. **P. Gao, G. Pedrini, C. Zuo, and W. Osten, “Phase retrieval using spatially modulated illumination,” Opt. Lett. **39**, 3615–3618 (2014). [CrossRef] [PubMed]

**31. **M. G. Gustafsson, “Nonlinear structured-illumination microscopy: wide-field fluorescence imaging with theoretically unlimited resolution,” Proc. Natl. Acad. Sci. U.S.A. **102**(37), 13081–13086 (2005). [CrossRef] [PubMed]

**32. **J. Mertz and J. Kim, “Scanning light-sheet microscopy in the whole mouse brain with HiLo background rejection,” J. Biomed Opt. **15**(1), 016027 (2010). [CrossRef] [PubMed]

**33. **M. R. Arnison, K. G. Larkin, C. J. Sheppard, N. I. Smith, and C. J. Cogswell, “Linear phase imaging using differential interference contrast microscopy,” J. Microsc. **214**(1), 7–12 (2004). [CrossRef] [PubMed]

**34. **S. Velghe, J. Primot, N. Guérineau, M. Cohen, and B. Wattellier, “Wave-front reconstruction from multidirectional phase derivatives generated by multilateral shearing interferometers,” Opt. Lett. **30**(3), 245–247 (2005). [CrossRef] [PubMed]

**35. **S. Ettl, J. Kaminski, M. C. Knauer, and G. Häusler, “Shape reconstruction from gradient data,” Appl. Opt. **47**(12), 2091–2097 (2008). [CrossRef] [PubMed]

**36. **S. Chowdhury and J. Izatt, “Structured illumination quantitative phase microscopy for enhanced resolution amplitude and phase imaging,” Biomed. Opt. Express **4**(10), 1795–1805 (2013). [CrossRef] [PubMed]