Translator Disclaimer
27 February 2018 Variable aperture-based ptychographical iterative engine method
Author Affiliations +
A variable aperture-based ptychographical iterative engine (vaPIE) is demonstrated both numerically and experimentally to reconstruct the sample phase and amplitude rapidly. By adjusting the size of a tiny aperture under the illumination of a parallel light beam to change the illumination on the sample step by step and recording the corresponding diffraction patterns sequentially, both the sample phase and amplitude can be faithfully reconstructed with a modified ptychographical iterative engine (PIE) algorithm. Since many fewer diffraction patterns are required than in common PIE and the shape, the size, and the position of the aperture need not to be known exactly, this proposed vaPIE method remarkably reduces the data acquisition time and makes PIE less dependent on the mechanical accuracy of the translation stage; therefore, the proposed technique can be potentially applied for various scientific researches.



High-resolution microscopic imaging presents fine details of the samples and becomes an indispensable tool for biological observation, medical diagnosis, and material analysis.1,2 To realize ideal imaging quality and spatial resolution, samples should always be sliced before being observed. Generally, these sample slices have the thickness of several micrometers and are roughly transparent objects, to obtain a high-contrast image with a common wide-field microscope, sample staining is required for most applications, though it not only requires complicated sample pretreating but also may damage the sample itself. Phase imaging is an alternative technique to realize high-contrast imaging without sample staining, where the spatial features of the samples are shown with the phase retardation of the transmitted light rather than the energy attenuation. Digital holography is a well-developed quantitative phase imaging technique, which can reconstruct both the amplitude and phase images of the sample simultaneously.37 Compared to the phase contrast microscopy and differential interference contrast microscopy, which are mainly for qualitative observations, digital holography can provide quantitative information on the sample amplitude and phase distributions. However, the additional reference beam makes the imaging system complicated and sensitive to environmental disturbances. Transport of intensity equation (TIE)-based phase microscopy is a noninterference phase sensing technique,815 which has a much lower requirement on the stability of the imaging system and working environment. However, like other lens-based imaging techniques, the image quality and accuracy of TIE heavily depend on the quality of imaging optics.

Coherent diffraction imaging (CDI), which is especially developed to solve the problems related to the lack of high-quality optics in the regime of x-ray, can retrieve both the amplitude and the phase of samples from recorded diffraction patterns via iterative computations. Since no lens or other complex optics is required, its resolving capability is only limited by the numerical aperture of the imaging system, and the diffraction limited resolution can be reached in theory. Unfortunately, the performance of traditional CDI algorithms, including the Gerchberg–Saxton algorithm,16 Fienup’s error reduction, hybrid input–output algorithms,17,18 and axial multi-intensity algorithm,19 are not ideal in terms of poor reliability, small field of view, and slow convergence speed. Ptychographical iterative engine (PIE),2026 which scans the object through a localized illuminating probe to a raster of positions and records all the diffraction patterns formed in the far field, can reconstruct the transmission function of the sample and the complex amplitude of the illumination accurately and rapidly with two counterpart updating formulas from the recorded diffraction patterns array. Compared to traditional CDI algorithms, PIE has outstanding advantages of fast convergence speed, high reconstruction quality, and theoretically infinite field of view, and accordingly it is regarded as a breakthrough of the CDI technique. PIE can also be adopted in the Fourier frequency domain to realize the so-called Fourier ptychographic microscopy,27,28 which does not require purely coherent illumination and extends the application range of PIE remarkably. Additionally, considering the reconstruction of other CDI techniques is only the exiting field from the sample, PIE can reconstruct the sample and the illumination separately and simultaneously, so the reconstructed image is free of the aberration from the illumination beam, and it lowers the requirement on the quality of the illumination optics. This is another outstanding advantage of PIE over other CDI techniques.

In the data acquisition, the sample should be scanned to many positions relative to a static illumination probe with a proper overlapping ratio between two neighboring illuminated regions.29 High information redundancy introduced by a vast number of diffraction patterns is one prerequisite for high-quality reconstruction in PIE. Since about 10×10 diffraction patterns are recorded in most experiments to obtain satisfactory reconstruction, several minutes are often needed for a mechanical translation stage-based PIE system to finish the data acquisition, and high stability is required for both imaging system and working environment. Therefore, speedy PIE design using fast scanning devices to record as few diffractions patterns as possible is required. A kind of aperture scanning PIE was proposed recently,30 where a parallel beam is incident on a tiny aperture to form the illumination on the sample, and faithful information of the sample can be reconstructed from few diffraction intensities sequentially recorded when the aperture is changed to various sizes step by step. Since many fewer diffraction patterns are required, the data acquisition time can be remarkably reduced; however, various experimental problems should be solved before this method can be realized in practices. For example, it is impossible to obtain a perfectly planar illuminating beam on the variable aperture, and both the position and the shape of the aperture cannot be known exactly in experiments. In this paper, an improved variable aperture-based PIE (vaPIE) is designed, and the variable aperture is generated with a digital micromirror device (DMD), at the same time, the reconstruction algorithm is correspondingly modified to make it applicable for real experiments. To be specific, unlike the original algorithm processing the reconstruction between the detector plan and the sample plane,30 the modified algorithm carries out the reconstruction among the aperture plane, the sample plane, and the detector plane to retrieve the information of both sample and aperture, thus the illumination beam on the variable aperture needs not to be ideally planar, as well as both the shape and the center of the aperture. The underlying physics for this method to realize fast PIE imaging with less diffraction patterns is similar to that of the probe-diversity ptychography,31 which applies two probes to record the diffraction patterns and then can mathematically change the searching step more effectively than standard single-probe PIE.32,33 The vaPIE achieves the probe diversity with more probes using size changing apertures; therefore, it is more robust in avoiding the stagnation of local minimizations during the iterative computation.

In this paper, the feasibility of the vaPIE is demonstrated both numerically and experimentally, and it is found that the satisfactory reconstruction can be obtained from only 12 diffraction patterns recorded within 3 s. Moreover, both the shape and position of the apertures need not to be known exactly. The principle of the vaPIE is analyzed in Sec. 2, and the feasibility of this method is verified both numerically in Sec. 3 and experimentally with visible light in Sec. 4.



The principle of vaPIE is shown in Fig. 1, where Fig. 1(a) shows the data collection procedure and Fig. 1(b) shows the reconstruction process. In Fig. 1(a), the roughly parallel beam is incident on the n’th (n=1,,N) aperture Sn(r)={1,r<rn0,rrn, and the exiting light from the aperture forms the probe Pn(r) incident on the sample after propagating a distance of d1. The diffraction intensity In(r) is recorded by an image recorder with the distance d2 behind the sample. Assuming the transmitting function of the sample and the illumination on the aperture as O(r) and E(r), respectively, the illumination on the sample and the diffraction at the recording plane can be calculated as Pn(r)=Fd1[E(r)Sn(r)] and Un(r)=Fd2[Pn(r)O(r)]=|Un(r)|·exp[iφn(r)], where Fd represents the Fresnel propagation of distance d.

Fig. 1

(a) Scheme of vaPIE and (b) flowchart of the retrieval algorithm.


The iterative reconstructing computation is carried out with the following steps after two uniform initial guesses are given to the illumination E(r) and the sample function O(r), respectively.

  • 1. In the m’th round iterative reconstruction, the n’th illuminating probe Pm,n(r) on the sample is calculated as

    Eq. (1)


  • 2. The transmitting field of the sample Tm,n(r) and the complex field at the recording plane are computed as

    Eq. (2)


  • 3. Replace the modulus of Um,n(r) with the square root of the recorded diffraction intensity In(r) and keep its phase unchanged as

    Eq. (3)


  • 4. Backpropagate the updated complex field Um,n(r) to the sample plane as

    Eq. (4)


  • 5. Update the probe light Pm,n(r) and the sample function with Eq. (5), where both α and β are constants with values between 0 and 1

    Eq. (5)


  • 6. Propagate Pm,n(r) to the aperture plane and update the illuminating field Em(r) as

    Eq. (6)


  • 7. Jump to step 1 until all the aperture sizes are addressed.

  • 8. Calculate the residual error ϵ with Eq. (7), if the accuracy meets the requirement, the iteration stops, or else jump to step 1 to start another round of iterative computation

    Eq. (7)


Since the shape of the aperture forming the illuminating probes P(r) cannot be known exactly in practice, in the iterative computations, the aperture size used in Eq. (6) should be slightly larger than the real value, and it will be discussed in detail in Sec. 3.


Numerical Simulations

To check the performance of the vaPIE, numerical simulations were used to determine the main parameters for real experiments. In simulations, the distances of the sample to the aperture plane and that to the recording plane were set as 315 and 18 mm, respectively. The wavelength was 632.8 nm, and the pixel size of the detector was 7.4  μm. The diameter of the circular aperture was scanned from 0.74 to 2.54 mm with an increment step of 0.15 mm, and 12 frames of diffraction patterns were captured for reconstructing the complex transmitting function of the sample. Figure 2(a) shows the setting amplitude and phase distributions of the sample, and Fig. 2(b) exhibits all the diffraction patterns formed at the recording plane. With the above described reconstruction algorithm, the reconstructed phase and amplitude distributions with 50 iterations are shown in Fig. 2(c), where the residual error calculated by Eq. (7) was only about 1.0%. Since only 12 diffraction patterns are required, compared to the common PIE, which often needs 100 or more diffraction patterns and about 10 min for data acquisition, the data acquisition of this method is remarkably accelerated, making the vaPIE more applicable in practice.

Fig. 2

Numerical simulation of vaPIE: (a) presetting sample amplitude and phase distributions, (b) diffraction patterns with different variable apertures, and (c) retrieved sample amplitude and phase distributions. The color bar of the amplitude indicates the normalized magnitude, and that of the phase presents the value in radians.


In the simulations shown in Fig. 2, both the shape and the position of each tiny aperture are assumed to be ideally accurate; however, in real experiments, both of them cannot be known exactly. Thus, the influence of the shape and position error on the reconstruction accuracy should be quantitatively evaluated. First, the influence of the aperture deviation was discussed in detail. Figure 3(a) shows the comparison between the apertures used for diffraction patterns calculation and sample reconstruction, indicating the definition of the aperture deviation. Figure 3(b) shows the convergence speed changes with different aperture deviations, where the reconstruction error increases with the inaccuracy of the aperture size. When the size of the reconstructing aperture is 20% larger than that of the real aperture, the reconstruction error increases from 1% to about 2%, and the reconstructed image shown in Fig. 3(c) is still of high quality; but when the size of the reconstructing aperture is 33.3% larger than that of the real aperture, the reconstruction error increases to 6%, and the reconstructed image is seriously blurred. This means that in real experiments, good reconstruction can still be available when the size and the shape of the aperture are only approximately known.

Fig. 3

Numerical simulation on aperture deviation: (a) definition of the aperture deviation, (b) retrieved sample amplitude and phase distributions with different aperture deviations, and (c) quantitative evaluation on retrieval quality according to the residual error. The color bar of the amplitude indicates the normalized magnitude, and that of the phase presents the value in radians.


The positions or the centers of all the apertures are difficult to keep unchanged in the whole data acquisition in experiments, in other words, the positions of the apertures cannot be known exactly in practices. To find out the influence of the center deviation on the reconstruction quality, another set of simulations was carried out in Fig. 4, where Fig. 4(a) shows the definition of the center deviation Δd, and Fig. 4(b) shows the convergence and the residual reconstruction error according to various center deviations. It is found that with the increasing center deviation, the convergence speed obviously slows down and the residual error in the reconstructions increases remarkably. However, for the case of center deviation as 6% shown in Fig. 4(c), the residual error in the final reconstruction is only about 3.5%, proving that the proposed vaPIE is still robust in sample retrieval even with existing center deviations.

Fig. 4

Numerical simulation on center deviation: (a) definition of the center deviation, (b) retrieved sample amplitude and phase distributions with different center deviations, and (c) quantitative evaluation on retrieval quality according to the residual error. The color bar of the amplitude indicates the normalized magnitude, and that of the phase presents the value in radians.


During the data acquisition, the size of the aperture forming the illumination on the sample increases step by step, Fig. 5 shows the simulation results reflecting the influence of the step size Δr on the reconstruction quality, where Fig. 5(a) shows the definition of Δr, and Fig. 5(b) shows the convergence speeds and the residual error of the reconstruction when the aperture size increases from 0.74 to 2.54 mm with Δr=0.05, 0.10, 0.15, 0.20, and 0.25 mm, respectively. It is found that the reconstruction error increases slowly with the step size, for the step size of 0.05 mm, the residual error is <1%, and for the step size of 0.25 mm, the residual error is about 6%. The underlying physics for the influence of step size on the reconstruction quality is that, with the increasing step size, fewer diffraction patterns and accordingly less information on the object are recorded and then the reliability of the reconstruction degrades.

Fig. 5

Numerical simulation on increment of the variable-aperture size: (a) definition of the increment of the variable-aperture size, (b) retrieved sample amplitude and phase distributions with different increments of the variable-aperture size, and (c) quantitative evaluation on retrieval quality according to the residual error. The color bar of the amplitude indicates the normalized magnitude, and that of the phase presents the value in radians.


The numerical simulations illustrated above show that the proposed vaPIE can reconstruct both the phase and amplitude images simply from diffraction patterns obtained using a variable aperture to form the changing illumination on the sample. Moreover, according to the numerical simulation, the size, the shape, and the position need not to be known exactly in the reconstruction, and the step size of the changing aperture should be as small as possible to record as many diffraction patterns as needed to achieve faithful reconstruction.



To verify the feasibility of the proposed method experimentally, a series of experiments were proposed with visible light. Figure 6(a) shows the experimental setup, where the illumination on the object was a parallel He-Ne laser beam reflected from a DMD (DLP Discovery D4100, Texas Instruments). Since DMD can be programmed into reflecting mirrors of any shape to reflect the laser beam to the sample, it acted as the variable aperture in experiments. The top row of Fig. 6(b) shows three apertures of different sizes, where the white regions within the red circles are the effective reflecting regions and the bottom row shows the corresponding recorded diffraction patterns. The sample was placed 315 mm behind the DMD, and a CCD camera (AVT Pike F421B, Germany) positioned 18 mm behind the sample was used to record the diffraction patterns. During the data acquisition, the radius of the aperture increased from 54 pixels (0.74 mm) to 185 pixels (2.54 mm), and the radius increment was 11 pixels (0.15 mm). The exposure time of the CCD was set to 200 ms, and the whole data acquisition was finished within 3 s to record 12 diffraction patterns. The final field of view is roughly estimated as 22π12  mm2 determined by the central region of the largest aperture. Note that the largest field of view achievable for this method is mainly limited by the size of the DMD, and this can be a shortcoming compared to common PIE, which has infinite field of view in theory. However, since DMD has a size of 10  mm×14  mm, this size is large enough for most applications.

Fig. 6

Experimental setup of the vaPIE: (a) vaPIE system and (b) variable-aperture scanning tactic realized by DMD.


A 1951 USAF resolution target (Edmund Optics) was first used to test the resolution of the imaging system. Figure 7(a) lists partial diffraction patterns recorded with different apertures. Figure 7(b1) is the reconstructed image and Fig. 7(b2) is the zoomed-in image, where we can find that the third element of group six can be clearly resolved as shown in Fig. 7(c), corresponding to a resolution of 6.20  μm. Since the probe generated by a small aperture is remarkably smaller than that of a large aperture, the image will be updated more times at the center than at the edge in the iterative computation from steps 1 to 7, the central region becomes clear first in the iterative computation, but after many rounds of iterations, the region close to the edge will also be updated for enough times and then becomes clear at last.

Fig. 7

Experimental measurements on a standard 1951 USAF resolution test chart using vaPIE: (a) diffraction patterns on the detector plane and the corresponding illumination beams on the sample plane, (b1) and (b2) retrieved amplitudes of 1951 USAF resolution test chart in different field of views, and (c) intensity distributions along the red and blue lines in (b2).


Next, two fixed biological samples were imaged to check the performance of the vaPIE: one was the corn section and another was the monocotyledon section. Figure 8 shows the reconstructed phase and amplitude images, where Figs. 8(a1) and 8(b1) are the two amplitude images, and Figs. 8(a2) and 8(b2) are the two corresponding phase images. It is found that the structure of each individual cell can be clearly identified, proving the high reconstruction quality of the proposed vaPIE. The structure of the deemed concentric circles can be found in these reconstructed images, and according to the simulations shown Fig. 5, this is the reconstruction noise related to the big increment step in changing the aperture size and can disappear when the finer increment step is used.

Fig. 8

Experimental measurements on biological samples as (a) corn section and (b) monocotyledon section. (1) and (2) Retrieved amplitude and phase distributions via vaPIE. The color bar of the amplitude indicates the normalized magnitude, and that of the phase presents the value in radians. The black bar in (a1) indicates 0.3 mm.




An improved vaPIE is proposed to reconstruct the sample amplitude and phase images, and its feasibility is verified both numerically and experimentally. Using a DMD to generate the aperture with variable sizes, the data acquisition of 12 frames of diffraction patterns can be finished within 3 s. Moreover, combining with the modified reconstruction algorithm to carry out the iterative computation among the aperture plane, the sample plane and the detector plane, satisfactory reconstruction can be obtained within 50 iterations. In other words, the vaPIE can realize speedy microscopic quantitative images by remarkably reducing the data acquisition and reconstruction time. At the same time, since the shape, the size, and the position of the variable aperture need not to be known exactly, the vaPIE is less dependent on the mechanical accuracy of the imaging system. Considering its advantages of fast convergence speed, rapid data acquisition, high reconstruction quality, and simplified setup, the proposed vaPIE can also be applied in many applications for scientific research.


All of the authors have no relevant financial interests in this article and no potential conflicts of interest to disclose.


This research was supported by the National Natural Science Foundation of China (Nos. U1730132, 61705092, 31522056, and 11647144), the Natural Science Foundation of Jiangsu Province of China (Nos. BK20130162 and BK20170194), the Shanghai Sailing Program (No. 17YF1407000), and the Fundamental Research Funds for the Central Universities (Nos. JUSRP115A14 and JUSRP51721B).



M. Mir et al., “Quantitative phase imaging,” Prog. Opt., 57 133 –217 (2012). POPTAN 0079-6638 Google Scholar


H. Majeed et al., “Breast cancer diagnosis using spatial light interference microscopy,” J. Biomed. Opt., 20 (11), 111210 (2015). JBOPFO 1083-3668 Google Scholar


V. Bianco et al., “Quasi noise-free digital holography,” Light Sci. Appl., 5 e16142 (2016). Google Scholar


F. Merola et al., “Tomographic flow cytometry by digital holography,” Light Sci. Appl., 6 e16241 (2017). Google Scholar


A. Ozcan and E. McLeod, “Lensless imaging and sensing,” Annu. Rev. Biomed. Eng., 18 (18), 77 –102 (2016). ARBEF7 1523-9829 Google Scholar


D. Li et al., “Accurate retrieval algorithm of amplitude from radial-shearing interferogram,” Opt. Lett., 35 (18), 3054 –3056 (2010). OPLEDP 0146-9592 Google Scholar


X. Wu et al., “Deviation influences on sectional image reconstruction in optical scanning holography using a random-phase pupil,” Appl. Opt., 52 (1), A360 –A366 (2013). Google Scholar


X. Tian et al., “Real-time quantitative phase imaging based on transport of intensity equation with dual simultaneously recorded field of view,” Opt. Lett., 41 (7), 1427 –1430 (2016). OPLEDP 0146-9592 Google Scholar


X. Tian et al., “In-focus quantitative intensity and phase imaging with numerical focusing transport of intensity equation method,” J. Opt., 18 (10), 105302 (2016). Google Scholar


W. Yu et al., “Real time quantitative phase microscopy based on single-shot transport of intensity equation (ssTIE) method,” Appl. Phys. Lett., 109 (7), 071112 (2016). APPLAB 0003-6951 Google Scholar


X. Meng et al., “Smartphone based hand-held quantitative phase microscope using the transport of intensity equation method,” Lab Chip, 17 104 –109 (2017). LCAHAM 1473-0197 Google Scholar


L. Waller, L. Tian and G. Barbastathis, “Transport of intensity phase-amplitude imaging with higher order intensity derivatives,” Opt. Express, 18 (18), 12552 –12561 (2010). OPEXFF 1094-4087 Google Scholar


L. Tian, J. Petruccelli and G. Barbastathis, “Nonlinear diffusion regularization for transport of intensity phase imaging,” Opt. Lett., 37 (19), 4131 –4133 (2012). OPLEDP 0146-9592 Google Scholar


S. Kou et al., “Transport-of–intensity approach to differential interference contrast (TI-DIC) microscopy for quantitative phase imaging,” Opt. Lett., 35 (3), 447 –449 (2010). OPLEDP 0146-9592 Google Scholar


L. Waller et al., “Transport of intensity phase imaging in a volume holographic microscope,” Opt. Lett., 35 (17), 2961 –2963 (2010). OPLEDP 0146-9592 Google Scholar


R. W. Gerchberg and W. O. Saxton, “A practical algorithm for the determination of phase from image and diffraction plane pictures,” Optik, 35 227 –246 (1972). OTIKAJ 0030-4026 Google Scholar


J. Fienup, “Reconstruction of an object from the modulus of its Fourier transform,” Opt. Lett., 3 (1), 27 –29 (1978). OPLEDP 0146-9592 Google Scholar


J. Fienup, “Phase retrieval algorithms: a comparison,” Appl. Opt., 21 (15), 2758 –2769 (1982). Google Scholar


B. Li et al., “An enhanced reconstruction algorithm to extend CT scan field-of-view with z-axis consistency constraint,” Med. Phys., 39 (10), 6028 –6034 (2012). MPHYA6 0094-2405 Google Scholar


J. Rodenburg and H. Faulkner, “A phase retrieval algorithm for shifting illumination,” Appl. Phys. Lett., 85 (20), 4795 –4797 (2004). APPLAB 0003-6951 Google Scholar


H. Faulkner and J. Rodenburg, “Movable aperture lensless transmission microscopy: a novel phase retrieval algorithm,” Phys. Rev. Lett., 93 (2), 023903 (2004). PRLTAO 0031-9007 Google Scholar


H. Faulkner and J. Rodenburg, “Error tolerance of an iterative phase retrieval algorithm for moveable illumination microscopy,” Ultramicroscopy, 103 (2), 153 –164 (2005). ULTRD6 0304-3991 Google Scholar


J. Rodenburg, “Transmission microscopy without lenses for objects of unlimited size,” Ultramicroscopy, 107 (2–3), 227 –231 (2007). Google Scholar


P. Thibault et al., “Probe retrieval in ptychographic coherent diffractive imaging,” Ultramicroscopy, 109 (4), 338 –343 (2009). ULTRD6 0304-3991 Google Scholar


A. Maiden and J. Rodenburg, “An improved ptychographical phase retrieval algorithm for diffractive imaging,” Ultramicroscopy, 109 (10), 1256 –1262 (2009). ULTRD6 0304-3991 Google Scholar


W. Yu et al., “High-quality image reconstruction method for ptychography with partially coherent illumination,” Phys. Rev. B, 93 (24), 241105 (2016). Google Scholar


G. Zheng, R. Horstmeyer and C. Yang, “Wide-field, high-resolution Fourier ptychographic microscopy,” Nat. Photonics, 7 (9), 739 –745 (2013). NPAHBY 1749-4885 Google Scholar


C. Kuang et al., “Digital micromirror device-based laser-illumination Fourier ptychographic microscopy,” Opt. Express, 23 (21), 26999 –27010 (2015). OPEXFF 1094-4087 Google Scholar


O. Bunk et al., “Influence of the overlap parameter on the convergence of the ptychographical iterative engine,” Ultramicroscopy, 108 (5), 481 –487 (2008). ULTRD6 0304-3991 Google Scholar


T. Shimobaba et al., “Ptychography by changing the area of probe light and scaled ptychography,” Opt. Commun., 331 (22), 189 –193 (2014). OPCOB8 0030-4018 Google Scholar


I. Peterson, R. Harder and I. Robinson, “Probe-diverse ptychography,” Ultramicroscopy, 171 77 –81 (2016). ULTRD6 0304-3991 Google Scholar


C. Yang et al., “Iterative algorithms for ptychographic phase retrieval,” (2011). Google Scholar


J. Qian et al., “Efficient algorithms for ptychographic phase retrieval,” Inverse Problems and Applications, 615 261 –280 American Mathematical Society(2014). Google Scholar

Biographies for the authors are not available.

© 2018 Society of Photo-Optical Instrumentation Engineers (SPIE) 1083-3668/2018/$25.00 © 2018 SPIE
Aihui Sun, Yan Kong, Xin Meng, Xiaoliang He, Ruijun Du, Zhilong Jiang, Fei Liu, Liang Xue, Shouyu Wang, and Cheng Liu "Variable aperture-based ptychographical iterative engine method," Journal of Biomedical Optics 23(2), 026007 (27 February 2018).
Received: 6 November 2017; Accepted: 26 January 2018; Published: 27 February 2018

Back to Top