## Content adaptive sparse illumination for Fourier ptychography

Liheng Bian, 1 Jinli Suo, 1 Guohai Situ, 2 Guoan Zheng, 3 Feng Chen, 1 and Qionghai Dai 1 , ∗

1 Department of Automation, Tsinghua University, China

2 Shanghai Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, China

3 Biomedical Engineering &amp; ∗

Electrical and Computer Engineering, University of Connecticut, USA

Corresponding author: qhdai@tsinghua.edu.cn

## Compiled August 5, 2014

Fourier Ptychography (FP) is a recently proposed technique for large field of view and high resolution imaging. Specifically, FP captures a set of low resolution images under angularly varying illuminations and stitches them together in Fourier domain. One of FP's main disadvantages is its long capturing process due to the requisite large number of incident illumination angles. In this letter, utilizing the sparsity of natural images in Fourier domain, we propose a highly efficient method termed as AFP, which applies content adaptive sparse illumination for Fourier ptychography by capturing the most informative parts of the scene's spatial spectrum. We validate the effectiveness and efficiency of the reported framework with both simulations and real experiments. Results show that the proposed AFP could shorten the acquisition time of conventional FP by around 30%-60%. c © 2014 Optical Society of America

OCIS codes: (110.1758) Computational imaging; (110.3010) Image reconstruction techniques; (170.0180) Microscopy;.

Fourier Ptychography (FP) is a newly proposed technique for large field of view (FOV) and high resolution (HR) imaging [1-3]. This method sequentially captures a set of low resolution (LR) images describing different spatial spectrum bands of the light field, and these spectrum bands retrieved from captured LR images [4] are stitched together in Fourier domain to reconstruct the entire HR spatial spectrum and corresponding HR image. Recently, Zheng et al. [1] extend FP technique to microscopic imaging and propose Fourier Ptychography Microscopy (FPM). FPM assumes that the incident light of a microscope is plane wave, and thus different incident angles result in different shifted versions of the spatial spectrum of the light field through specimen. Since the whole light field is filtered by the microscope's objective and sampled by CCD, microscopy under changing incident illumination angles results in a set of images corresponding to different spectrum bands of the entire light field. According to [1], it takes the FPM setup about 3 minutes to capture enough images under various illumination angles for reconstructing a gigapixel grayscale image, whose optical resolution is ∼ 0.78 µ m and FOV as ∼ 120 mm . If one intends to obtain a color image, 2 the capturing time would be multipled. What's more, the calculation time will also increase linearly with the number of captured images.

There have been several studies improving FPM and FP from mainly two aspects. The first aspect is to correct various system aberrations. Bian et al. [5] derive an adaptive wavefront correction framework for FPM to minimize the error between the captured intensity image set and reconstructed image. This method can correct illumination fluctuation and pupil aberration effectively, but the global optimization process is time consuming. For efficiency improvement, a fast embedded pupil function recovery algorithm for FPM is proposed in [6] based on the extended ptychographic iterative engine [7], which

Fig. 1: Illustration of the sparsity prior of natural images in Fourier domain. The relationship between reserved ratio in Fourier domain and recovered RMSE in spatial domain of three kinds of public images including common natural scenes (USC-SIPI dataset), lung tissues and blood cells (example images and corresponding spatial spectra are shown on the left) is plotted in the right figure. Besides, reconstructed errors of the 'Lena' image with different reserved ratios are presented in the figure as well (amplified by 5 times for better visualization).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0577v1-figure-f39545c6-f3f7-4e8b-ad26-804e52e40d16.jpg)

is designed for fast probe function correction in conventional ptychography. The second kind of improvements is to reduce the hardware costs and system's running time. Multiplexing illuminations are applied to FP [8,9] and high quality HR images could be reconstructed using less shots. Dong et al. [10] investigate the requisite spatial data redundancy of FP and report an FP scheme that needs only one shot for each illumination angle instead of multiple differently exposed images. However, this study does not reduce the requisite large number of incident illumination angles, which also leads to long acquisition time of FPM and other FP setups.

To address the efficiency problem in current FPM, this letter investigates the spectrum redundancy of natural images in Fourier domain and proposes to stitch sparsely sampled but most informative subregions of spatial spectrum to reconstruct the latent HR image. Sparse sam-

Fig. 2: Flow chart of the proposed APF framework. After initialization as the left part shows, final results can be exported by iteratively and adaptively updating the HR spectrum (as the blue hollow arrows show) using a small number of images captured in the first iteration (as the red solid arrows show) under different incident illuminations which correspond to the redlabeled spectrum bands shown in the central part.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0577v1-figure-7c61639f-328a-4e99-b4fa-af12deb1f1d7.jpg)

pling strategy will decrease the requisite illumination angles during acquisition, and correspondingly decrease the capturing time of FP technique.

To study the redundancy [11] in the spatial spectrum of a natural image statistically, we transform several public image data sets including both macroscopic (USC-SIPI common miscellaneous dataset [12]) and microscopic scenes (lung tissues and blood cells [13]) to Fourier domain, and then transform them back after truncating small values and calculate root-mean-square error (RMSE) of the reconstructed images. The relationship between reserved ratio of Fourier coefficients and RMSE is plotted in Fig. 1, where we also display the spatial spectrum of several exemplar images for further analysis. Besides the well known and widely existing sparsity of natural images' Fourier coefficients, they also exhibit some other properties. The spatial spectrum is apparently not uniform but exhibits inhomogeneous and monotonously decreasing trends along the radial direction. Besides, the spectrum is of some directionality whose distribution is closely related to the structure of the image. These inspiring observations motivate us to propose the adaptive Fourier ptychography (AFP) which does not exhaustively sample scene's Fourier domain as conventional FP [2], but utilizes sparse illuminations to adaptively sample the most important bands of the HR spectrum from central low frequency region outward towards the high frequency regions, i.e., in a circle-wise manner.

The entire framework of AFP is diagramed in Fig. 2. Before detailed explanation of the proposed AFP framework, we firstly define a term 'potential position', which defines the center of a subregion of a HR spatial spectrum that is potentially to be sampled. In addition, con- sidering different pixel numbers between the captured LR images and the recovered image, all the spectra in the entire framework are normalized by dividing corresponding pixel numbers. What's more, we respectively use x-y and u-v coordinate to describe an image in spatial domain and in Fourier domain. With the above notations, the framework of AFP is detailed as follows.

Firstly, the initial input of the reconstruction framework is set to be the zero-padding extension of spatial spectrum of the LR image shot under normal incident light. Besides, the initial potential position is set to be center of the HR spectrum.

Then we begin loops to update the HR spectrum. At the start of each loop, we determine potential positions of current loop in the HR Fourier domain by extending the positions of the previous loop towards the high frequency direction with a fixed step size (determined by the user-defined spectrum overlapping ratio), as the yellow dashed line shows in Fig. 2. After getting potential positions in current loop, we decide whether to update the subregion corresponding to each potential position. For each subregion, it will be updated only when it satisfies following two conditions simultaneously: 1) the maximum amplitude of the entries in this subregion is bigger than a manually set threshold, and 2) the closest distance from this subregion's center to other alreadyupdated centers is bigger than the step size. If one subregion is determined to be updated, the LR image under corresponding angular incident light should be captured as the red solid arrows show, and the subregion is then updated in the same way as FP. The correspondence between the subregion's center and the incident angle can be drawn from [14] as follows. If CCD's sampling rate, namely CCD's pixel size d , satisfies the sampling requirement d glyph[lessorequalslant] λ 2 NA [1] where λ is the wavelength of incident light and NA is the numerical aperture of the microscope's objective, then the correspondence is

$$\Delta u = \frac { d M } { \lambda } \cos ( \hat { \alpha } ) \text{ and } \Delta v = \frac { d N } { \lambda } \cos ( \hat { \beta } ), \quad ( 1 )$$

where ∆ u and ∆ v are the shifts of the subregion's center regarding to the center of the HR spectrum in Fourier domain; M and N are the pixel numbers of the captured LR images along two dimensions; ˆ and α ˆ are the inciβ dent direction angles regarding to the specimen plane. In addition, the filtering function of the objective is set to be the same as [1], i.e. a circular pupil whose radius is NA × 2 π λ in Fourier domain.

After traversing all the potential positions in current loop, the construction process moves to next loop until the potential positions reach the highest frequency of the HR spectrum, i.e., the spectrum boundary.

To improve the recovery accuracy further, users can repeat the above updating iteration in Fourier domain multiple times (experimentally ∼ 3 times are enough for a convergent result, which is the same as FP), as the blue hollow arrows in Fig. 2 illustrate. There are two worth noting difference between the repetitive iterations

Fig. 3: Comparison of simulation results by AFP and FP on images of different spectrum sparsities and directivities. (a) plots reconstructed errors of AFP (solid lines) and FP (dashed lines) at different spectrum reserved ratios. The blue, green and red lines respectively correspond to the experiment results conduced on Lena, lung tissue and blood cell image shown in Fig. 1. (b) shows example normalized reconstruction error images of AFP (left column) and FP (right column) with the same reserved ratio (25%) on the three images.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0577v1-figure-f6842e38-95bd-4c2d-aacc-bd0fba661ab2.jpg)

and the first one: (i) The previously captured LR images are utilized to avoid repetitive acquisition, i.e., we only adaptively capture LR images under angularly varying illuminations in the first iteration (as the red solid arrows show). (ii) The reconstruction result of the previous iteration instead of the initial image serves as the input of current iteration.

To measure the advantage of AFP over conventional FP quantitatively, here we respectively apply FP and the proposed APF on simulated inputs, and show pros and cons of the two reconstruction methods. The captured image volume is synthesized in following three steps: we do FFT to the original HR image, and select subregions corresponding to calculated incident angles by multiplying the HR spectrum with an ideal pupil function (all ones in the pupil circle and zeros outside the circle). Then by shifting these subregions to spectrum center and doing inverse FFT, we get complex images in spatial domain. Finally, we remove the phase of these images and retain only their amplitude as LR input images.

In the simulation, LR images' pixel numbers are set to be one tenth of the original HR image along each dimension, and overlapping ratio of spectrum subregions is 65%, which is the same as in [1]. Besides, the initial HR image is set to be all-one matrix, and the iteration number is set to be 5 times. The user-defined normalized updating threshold ranges from 10 256 2 to 100 256 2 with a step as 10 256 2 , which leads to the reserved ration ranging from around 70% down to around 5%. The simulation results are shown in Fig. 3. In (a), we compare AFP's reconstruction errors with that of FP on different images and different spectrum reservation ratios. The plots on three kinds of images consistently show that AFP is of less reconstruction errors than FP under the same reserved ratio. In other words, AFP needs less shots than FP to obtain the same recovery precision. Specifically, the amount of saved time is related to the sparsity and directionality of the scene's Fourier coefficients. For Lena and the lung tissue images which contain much high frequency and of strong directionality in Fourier domain, APF could save more than 50% shots, while for the blood cell image which owns less high frequency and directionality, AFP could save more than 30% shots. In Fig. 3(b), three groups of reconstructed error images are represented. The spectrum reserved ratio is set to be the same as 25%, and the errors are normalized by the maximum error value of FP for visualization. From the results we can clearly see that the recovered errors of AFP are much less than that of FP, which means that AFP preserves more crucial image details while of lower noise than FP under the same spectrum reserved ratio.

To further validate the proposed AFP, we build a microscopy setup to capture raw images, as shown in Fig. 4(a). The setup parameters are almost the same as the FPM setup in [1] except that Olympus BX43 microscope and 2 × objective (NA = 0.08) are utilized for our commercial availability. The LED plate is placed around 8cm under the specimen, and the lateral distance between two adjacent LEDs is 4 mm which is determined by the use-defined spectrum overlapping ratio (35%) and the objective. The central wavelength of incident light is 632nm. The pixel size in the final shot raw images is calculated by dividing the detector's pixel size ( ∼ 5.8 µ m) with the objective's magnification factor (2) as ∼ 2.9 µ m, which is smaller than the maximum pixel size of raw images required by the sampling requirement as λ 2 NA obj = 3 95 . µm . We capture three LR images for each illumination angle whose exposure time are respectively 0.005s, 0.05s and 0.5s, and then combine them to obtain a high-dynamic range image for subsequent reconstruction. In addition, the iteration is still set to be 5 times and the embedded pupil function correction technique [6] is utilized in the reconstruction process for smaller reconstruction error.

Results recovered by FP and the proposed AFP are respectively shown in Fig. 4(b) and Fig. 4(c), where the resolution of LR shot images are 100 × 100, and reconstruction magnification factor is 15. From the figure we can see that AFP works successfully on real captured images including both of the USAF chart and the mouse brain slice (fiber tissue). Specifically, the USAF recovery results show that with 85 LR images AFP could resolve the feature of group 9, element 3 on the USAF target, which achieves the same resolution obtained by FP with 225 LR images. Namely, around 60% shoots are saved by the proposed AFP method.

Both the above simulation and real-shot experiments indicate that AFP could significantly shorten the acquisition time compared to FP. To explain the source of AFP's efficiency advantages, we illustrate its sampling strategy in parallel with that of conventional FP in Fig. 5. As shown in (a), FP captures all the images corresponding to the central part of HR spectrum exhaustively and then updates the spectrum sequentially. There

Fig. 4: The prototype setup and real reconstruction results. (a) shows the prototype microscopy system, as well as two examples of LR images captured under normal illumination including USAF chart and mouse brain slice. (b) presents reconstructed spatial spectra (amplitude) and corresponding HR images by FP utilizing sequential 225 LR images, while the results by the proposed AFP are shown in (c), where the user-defined normalized threshold is set to be 1 5 . × 10 -3 (85 LR input images used for the USAF chart and 72 images used for the brain slice).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0577v1-figure-e1234551-33fe-41ab-90f7-10c9dee5c064.jpg)

Fig. 5: Comparison of Fourier domain sampling sequence beween FP (top row) and AFP (bottom row) on the Lena image. (a) shows the theoretical sampling sequences of FP and AFP. (b) and (c) respectively present the simulated sampling central positions and recovered spectra of Lena image by FP and AFP under the same sampling ratio as 25%.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0577v1-figure-44e6c69c-e1b5-4e61-9116-01ce05165e8c.jpg)

are two main disadvantages of FP's sampling method: 1) Some areas of small amplitudes such as the '1' labelled region in (a) are updated, but this improves little reconstruction performance, so capturing their corresponding images is inefficient. 2) Because FP only samples the central part of the whole spectrum with a light NA as ∼ 0.5 [1], some important high frequencies owning large amplitudes are neglected, such as the '2' labeled area.

Instead, our proposed AFP adaptively updates the spectrum subregions in a circle-wise manner from low frequency to high frequency. In each circle it only updates areas of large spectrum amplitudes. For example, as Fig. 5(b) and (c) show, the horizonal outward regions of bigger amplitudes are updated selectively, while the regions of smaller amplitudes (such as the top and bottom high frequency regions, as well as some sporadic low frequency regions such as the red box labeled regions) are neglected. To sum up, benefiting from the adaptive updating scheme with highly efficient implementation, AFP can capture the most important information in Fourier domain and thus can save the acquisition time significantly, especially when object images own stronger directivity and thus more sparse spatial spectra.

We thank Xuemei Hu, Dongsheng An, Jingyu Lin, Weixin Jiang and Yulun Zhang for valuable discussions. This work is supported by the National Nature Science Foundation of China under grant 61327902, 61377005 and 61120106003, as well as the Recruitment Program of Global Youth Experts.

## References

- 1. G. Zheng, R.Horstmeyer, and C. Yang, Nat. Photonics 7 , 739 (2013).
- 2. X. Ou, R.Horstmeyer, C. Yang, and G.Zheng, Opt. Lett. 38 , 4845 (2013).
- 3. G.Zheng, IEEE Photonics J. 6 , 1 (2014).
- 4. J. R. Fienup, Appl. Opt. 21 , 2758 (1982).
- 5. Z. Bian, S. Dong, and G. Zheng, Opt. Express 21 , 32400 (2013).
- 6. X. Ou, G.Zheng, and C. Yang, Opt. Express 22 , 4960 (2014).
- 7. A. M. Maiden, and J. M. Rodenburg, Ultramicroscopy 109 , 1256 (2009).
- 8. S. Dong, R. Shiradkar, R. Nanda, and G. Zheng, Biomed. Opt. Express 5 , 1757 (2014).
- 9. L. Tian, X. Li, K. Ramchandran, and L. Waller, Biomed. Opt. Express 5 , 2376 (2014).
- 10. S. Dong, Z. Bian, R. Shiradkar and G. Zheng, Opt. Express 22 , 5455 (2014).
- 11. M. W. Marcellin, Springer (2002).
- 12. Univeisity of Southern California, 'SIPI Image Database,' http://sipi.usc.edu/database/.
- 13. D. Chuckor, and H. Louisa, 'Electron Microscopy Images,' http://www.dartmouth.edu/%7Eemlab/gallery/.
- 14. J. W. Goodman, Introduction to Fourier optics (Roberts and Company Publishers, 2005).

## Informational Fourth Page

## References

- 1. G. Zheng, R. Horstmeyer, and C. Yang, 'Wide-field, high-resolution Fourier ptychographic microscopy,' Nat. Photonics 7 , 739-745 (2013).
- 2. X. Ou, R. Horstmeyer, C. Yang, and G. Zheng, 'Quantitive phase imaging via Fourier ptychographic microscopy,' Opt. Lett. 38 , 4845-4848 (2013).
- 3. G. Zheng, 'Breakthroughs in photonics 2013: Fourier ptychographic imaging,' IEEE Photonics J. 6 , 1-7 (2014).
- 4. J. R. Fienup, 'Phase retrieval algorithms: a comparison,' Appl. Opt. 21 , 2758-2769 (1982).
- 5. Z. Bian, S. Dong, and G. Zheng, 'Adaptive system correction for robust Fourier ptychographic imaging,' Opt. Express 21 , 32400-32410 (2013).
- 6. X. Ou, G. Zheng, and C. Yang, 'Embedded pupil function recovery for Fourier ptychographic microscopy,' Opt. Express 22 , 4960-4972 (2014).
- 7. A. M. Maiden, and J. M. Rodenburg, 'An improved ptychographical phase retrieval algorithm for diffractive imaging,' Ultramicroscopy 109 , 1256-1262 (2009).
- 8. S. Dong, R. Shiradkar, R. Nanda, and G. Zheng, 'Spectrum multiplexing and coherent-state decomposition in Fourier ptychographic imaging,' Biomed. Opt. Express 5 , 1757-1767 (2014).
- 9. L. Tian, X. Li, K. Ramchandran, and L. Waller, 'Multiplexed coded illumination for Fourier Ptychography with an LED array microscope,' Biomed. Opt. Express 5 , 2376-2389 (2014).
- 10. S. Dong, Z. Bian, R. Shiradkar and G. Zheng, 'Sparsely sampled Fourier ptychography,' Opt. Express 22 , 54555464 (2014).
- 11. M. W. Marcellin, 'JPEG2000: image compression fundamentals, standards, and practice,' Springer (2002).
- 12. Univeisity of Southern California, 'SIPI Image Database,' http://sipi.usc.edu/database/.
- 13. D. Chuckor, and H. Louisa, 'Electron Microscopy Images,' http://www.dartmouth.edu/%7Eemlab/gallery/.
- 14. J. W. Goodman, Introduction to Fourier optics (Roberts and Company Publishers, 2005).
Liheng Bian Jinli Suo Guohai Situ Guoan Zheng Feng Chen Qionghai Dai
Content adaptive sparse illumination for Fourier ptychography
Fourier Ptychography (FP) is a recently proposed technique for large field of
view and high resolution imaging. Specifically, FP captures a set of low
resolution images under angularly varying illuminations and stitches them
together in Fourier domain. One of FP's main disadvantages is its long
capturing process due to the requisite large number of incident illumination
angles. In this letter, utilizing the sparsity of natural images in Fourier
domain, we propose a highly efficient method termed as AFP, which applies
content adaptive sparse illumination for Fourier ptychography by capturing the
most informative parts of the scene's spatial spectrum. We validate the
effectiveness and efficiency of the reported framework with both simulations
and real experiments. Results show that the proposed AFP could shorten the
acquisition time of conventional FP by around 30%-60%.