Super-resolution iaging using the spatial-frequency filtered intensity fluctuation correlation Jane Sprigg 1,*, Tao Peng 1, and Yanhua Shih 1 arxiv:1409.134v3 [physics.optics] 1 Nov 016 1 University of Maryland, Baltiore County - 1000 Hilltop Circle, Baltiore MD, 150, United States of Aerica * jane.sprigg@gail.co ABSTRACT We report an experiental deonstration of a nonclassical iaging echanis with super-resolving power beyond the Rayleigh liit. When the classical iage is copletely blurred out due to the use of a sall iaging lens, by taking advantage of the intensity fluctuation correlation of theral light, the deonstrated caera recovered the iage of the resolution testing gauge. This ethod could be adapted to long distance iaging, such as satellite iaging, which requires large diaeter caera lenses to achieve high iage resolution. Introduction Iproving the resolution of optical iaging has been a popular research topic in recent years. 1 6 A coonly used siple approach is to easure the autocorrelation of two identical classical iages, effectively squaring the classical iage, I 1 (ρ 1 I 1 (ρ 1, where ρ 1 is the transverse coordinate of the detector. This autocorrelation produces a axiu gain of the spatial resolution. However, the iaging resolution of such a setup can be further iproved by changing the easureent fro I 1 (ρ 1 I 1 (ρ 1, in ters of intensity, or n 1 (ρ 1 n 1 (ρ 1, in ters of photon nuber counting, to the intensity fluctuation correlation I 1 (ρ 1 I (ρ, or n 1 (ρ 1 n (ρ, where ρ 1 and ρ are the transverse coordinates of two spatially separated detectors. Then, if only those fluctuation correlations due to the higher spatial frequencies fro I (ρ are selected, a super-resolving iage can be observed fro the joint detection of the intensity fluctuations at the two detectors. The physics behind this super-resolution is siilar to the original theral light ghost iaging, 7,8 and is quite different fro an autocorrelation easureent. It should be ephasized that the reported result is also different than that of Oh et. al.; 6 while the authors easured the intensity fluctuations, it was the intensity fluctuation autocorrelation I 1 (ρ 1, which was still liited by the resolution iproveent of an autocorrelation easureent. In this Report, we deonstrate a caera with resolution beyond the classical Rayleigh liit. Siilar to the original theral ghost iaging experients, 7,8 the caera produces an iage by the easureent of I 1 (ρ 1 I (ρ ; the caera consists of a typical iaging setup, except it has two sets of independent and spatially separated detectors: D 1 placed on the iage plane, and D placed on the Fourier transfor plane. Crucially, D integrates (sus only the higher spatial frequencies, or transverse wavevectors, by blocking the central area of the Fourier transfor plane. The iage is calculated fro the intensity fluctuations of D 1, at each transverse position ρ 1, and the bucket detector D. The easureent can be forulated as R c (ρ 1 = I 1 (ρ 1 dκ F(κ I (κ, where F(κ is a filter function which selects the higher spatial frequencies. Experiental setup Fig. 1 illustrates the laboratory-deonstrated caera setup. We used a standard narrow spectral bandwidth pseudo-theral light source consisting of a 10 diaeter 53 n wavelength laser bea scattered by illions of tiny diffusers on the surface of a rotating ground glass. The object iaged was the 5-3 eleent of a 1951 USAF Resolution Testing Gauge. Like a traditional caera, the iaging lens, L I, had an aperture liited by an adjustable pinhole to approxiately 1.36 diaeter. However, this iaging device has two optical ars behind its iaging lens L I. The light transitted by the object falls on a single-ode 0.005 diaeter fiber tip that scans the iage plane of ar one and is then interfaced with a photon counting detector D 1. The cobination of the scanning fiber tip, D 1, and photon counting detector acts as a CCD array. In ar two, the second lens L B is placed behind the iage plane and perfors a Fourier transfor of the field distribution of the iage plane of ar two. D, a scannable ultiode 0.105 diaeter fiber interfaced with a photon counting detector, is placed in the Fourier transfor plane of L B and integrates only the higher spatial frequencies while filtering out the lower spatial frequency
Figure 1. Experiental setup: a 10 diaeter 53 n wavelength laser bea scatters fro a rotating ground glass (i and strikes a 1951 USAF Resolution Testing Gauge (ii, then the iaging lens L I (iii and a pinhole of 1.36 diaeter (iv. The light splits to ar one, where it is collected by a scanning point detector D 1, and ar two, where it is collected by a spatially filtered bucket detector, consisting of a lens L B located on the iage plane and a ultiode fiber tip in the Fourier transfor plane. The iage is then calculated using the Photon Nuber Fluctuation Correlation (PNFC circuit. Figure. Resolution coparison for different iaging ethods of three 0.0141 wide slits iaged by a 10 diaeter source: (a unresolved first-order classical iage, where the gray shading arks the location of the slits; (b unresolved iages fro the fluctuation autocorrelation. The black dots indicates[i 1 (x 1 ] and the blue triangles show [ I 1 (x 1 ], as seen in Oh et. al.; 6 (c copletely resolved iage observed fro R c (x 1 where the solid line is a Gaussian data fitting. /7
odes; we ephasize that the placeent of this spatial filter does not depend on knowledge of the Fourier transfor of the iage. The intensity fluctuations in this experient were recorded by a Photon Nuber Fluctuation Correlation (PNFC circuit, 9,10 which independently records the arrival tie of each photo-detection event at D 1 or D. The intensity, easured by the nuber of photons detected per second, is divided into a sequence of short tie windows, each of which needs to be less than the second-order coherence tie of the light; it is iportant that the width of the tie window not be too long. The software first calculates the average intensity per short tie window, n s, where s=1, indicates the detector, and then the difference or fluctuation ter for each tie window: n j,s = n j,s n s. The corresponding statistical average of n 1 n is thus n 1 n = 1 [ ] N ( n j,1 n j,. (1 j It should be ephasized that when the fluctuation correlation is calculated between the two detectors, n 1 and n have not yet been tie-averaged. The tie averaging is perfored after the correlation, as indicated by the ultiplication appearing inside the su over j. Experiental results Typical experiental results are presented in Fig.. In this easureent, the 5-3 eleent of a 1951 USAF Resolution Test Target was iaged in one diension by scanning D 1 in the x-direction along the slits. Fig. (a shows a copletely unresolved classical iage of the three slits, I 1 (x 1, that was directly easured by the scanning detector D 1. For reference the gray shading indicates the location of the slits. Fig. (b shows two results: the black dots plot the autocorrelation [I 1 (x 1 ], while the blue triangles show the fluctuation autocorrelation [ I 1 (x 1 ] at each point. [ I 1 (x 1 ] was calculated fro the intensity fluctuation autocorrelation of D 1. The easureents in Fig. (b have a resolution gain, siilar to that of Oh et. al. 6 Using the Rayleigh liit 11,1 δx=1.λ s I /D, 11,1 where s I is the distance fro lens L I to the iage plane and λ is the wavelength of illuination, the expected resolution of the autocorrelation in the iage plane is approxiately δx/ =0.13 which, as seen in (b, is not enough to resolve the three slits which have a slit-to-slit separation of about 0.13. However, by spatially filtering ar two, the three slits of the 5-3 eleent of the gauge were clearly separated when correlated with ar 1, as seen in Fig. (c. The error bars in Figs. (a and (b are quite sall, especially when copared to those in Fig. (c. This is a typical negative feature of second-order easureents; copared to first-order classical iaging, in order to achieve the sae level of statistics the reported iaging echanis needs a longer exposure tie. How uch longer the easureent takes depends on the power of the light source and other experiental paraeters. Fig. (c is the su of two easureents obtained by placing the bucket fiber tip at two points in the Fourier transfor plane: x + = 0. and x = 0.4 ; this selects the higher spatial frequency odes which fall onto the two fiber tips and blocks all other spatial frequency odes. We represent this atheatically with the filter function F(x =Π(x x +,D F +Π(x x,d F to siulate the physical spatial filtering, where Π(x,w is a rectangle function of width w, D F = 0.105 is the fiber diaeter, and x is easured fro the central axiu. Then in one diension for κ kx / f B, where f B is the focal length of the bucket lens, R c (x 1 = I 1 (x 1 dx F(x I (x, and only integrates the higher spatial frequencies collected by the bucket detector in the neighborhood of x = x +,x. Again, this spatial filtering does not require any knowledge of the object or its Fourier transfor function. One way to iprove these results is to replace D with a CCD array; the CCD would still be in the Fourier transfor plane, but with the central pixels blocked. This would allow D to collect ore light of higher spatial frequencies. Although the liits of our equipent, software data storage, and tie constraints prevented the authors fro aking such iproveents, in this reported easureent all three slits of the resolution gauge are certainly well-resolved, while both the classical iaging and the autocorrelation echaniss could not resolve it. Discussion and theory In the experient, we use the spatial correlation of the noise, I 1 (r 1,t 1 I (r,t, to produce an iage fro the joint photo-detection of two independent and spatially separated photodetectors, D 1 and D. In the following, we outline the theory behind our experient. First we briefly consider how a first-order or classical caera produces an iage in its iage plane, I 1 (ρ 1. The experient was perfored using a pseudotheral light source, created by placing a rotating ground glass in the path of a laser bea. The ground glass contains a large nuber of tiny scattering diffusers which act as sub-sources. The wavepackets of scattered light play the role of subfields; each diffuser scatters a subfield to all possible directions, resulting in the subfields acquiring rando phases. Each sub-field propagates fro the source plane to the iage plane by eans of 3/7
a propagator or Green s function, E (ρ 1 = E g (ρ 1, where E is the initial phase and aplitude of the field eitted by sub-source and g (ρ 1 is the Green s function which propagates the light fro the th sub-source located at ρ to the point ρ 1 at soe distance z fro the source plane. To siplify the proble, we assue the fields are onochroatic and ignore the teporal part of the propagator. Then the light easured at coordinate (r,t is the result of the superposition of a large nuber of subfields, =1 E (r,t, each eitted fro a point sub-source, I(r,t= E (r,t n E n (r,t= E (r,t + E(r,tE n (r,t= I(r,t + I(r,t ( n where I(r, t, the ean intensity, is the result of the th subfield interfering with itself; I(r, t, the intensity fluctuation, is the result of the th subfield interfering with the nth subfield, n, and is usually considered noise because I(r,t =0 when taking into account all possible rando phases of the subfields. A classical iaging syste easures the ean intensity distribution on the iage plane, I 1 (ρ 1, where we have assued a point detector D 1 is placed at coordinate ρ 1, the transverse coordinate of the iage plane. In an ideal iaging syste, the self-interference of subfields produces a perfect point-to-point iage-foring function. The ideal classical iage assuing an infinite lens is the convolution between the aperture function of the object A(ρ O and the iage-foring δ-function which characterizes the point-to-point relationship between the object plane and the iage plane. 11 13 I(ρ 1 = E dρ O g (ρ O A(ρ O g O (ρ 1 ob j dρ O A(ρ O δ(ρ O + ρ 1 µ = A( ρ 1 /µ (3 where µ = s I /s O is the agnification factor, g (ρ O is a Green s function propagating the th subfield fro the source plane to the object plane over a distance z, and g O (ρ 1 is a function propagating the subfield fro the object plane to the detection plane over a distance s O + s I, and including the iaging lens. A(ρ O is an arbitrary function describing the object aperture. In reality, due to the finite size of the iaging syste, we rarely have a perfect point-to-point correspondence. Incoplete constructive-destructive interference blurs the point-to-point correspondence to point-to-spot correspondence. The δ-function in the convolution of Eq. 3 is then replaced by a point-to-spot iage-foring function, or a point-spread function which is deterined by the shape and size of the lens. For a lens with a finite diaeter, one coon odel describes the shape or pupil of the lens as a disk of diaeter D: I 1 (ρ 1 = dρ O A(ρ O sob [ π D ρ ob j ] (4 where the sobrero-like point-spread function is defined as sob(x J 1 (x/x; J 1 (x is the first-order Bessel function. The iage resolution is deterined by the width of the sob-function: the narrower the higher. A larger diaeter lens results in a narrower sob-function and thus produces iages with higher spatial resolution. To siplify the atheatics, it is coon to approxiate a finite lens as a Gaussian e (ρ L/(D/ with diaeter D, but a soother falloff than the disk approxiation. This leads to a Gaussian iaging-foring function: I 1 (ρ 1 = dρ O A(ρ O exp [ ( π D ρ ob j ]. (5 This Gaussian version of the iaging equation will be used later in nuerical calculations to siplify the atheatical evaluation. It is clear fro Eqs. 4 and 5 that for a chosen value of distance s O, a larger iaging lens and shorter wavelength will result in a narrower point-spread function, and thus a higher spatial resolution of the iage. Now we consider the noise produced iage that is observed fro the easureent of Fig. 1 by eans of I 1 (ρ 1 I (ρ. To ake the explanation of the experiental results easier to follow, first we exaine the case where two point scanning detectors D 1 and D are placed in the iage planes of ar one and ar two: I 1 (ρ 1 I (ρ = E (ρ 1 E n (ρ 1 n p q E p(ρ E q (ρ = E (ρ 1 E (ρ E n (ρ 1 E n(ρ =q n=p E (ρ 1 E (ρ (6 4/7
The calculation of E (ρ 1 E (ρ is straightforward: [ E (ρ 1 E (ρ = E dρ O g (ρ O dκa (κ, ρ O g O (κ, ρ 1 ][E dρ O g (ρ O ] dκ A(κ, ρ O g O (κ, ρ [ = E dρ O dρ O g (ρ O E g (ρ O dκ A (κ, ρ O e iκ ρ O sob [ π D ρ ]] [ dκ A(κ, ρ O e iκ ρ O sob [ π D ρ λ s O + ρ ] e ik(z 0+s O /(z 0 s O (ρo ρ O e ik/(s I(ρ1 ρ] Next, we coplete the suation over in ters of the subfields, or the sub-sources, by eans of an integral over the entire source plane. This integral results in the well-known Hanbury-Brown Twiss (HBT correlation: sob [(π θ/λ ρ O ρ O ], where θ is the angular diaeter of the light source relative to the object plane. To siplify further calculations, we assue a large value of θ and approxiate the sob-function to a δ-function evaluated at ρ O = ρ O, κ = κ. I 1 (ρ 1 I (ρ is therefore approxiately equal to: I 1 (ρ 1 I (ρ dρ O A(ρ O sob [ π D ρ ] sob [ π D ρ λ s O + ρ ] e ik/(s I(ρ1 ρ (8 It is clear that when ρ 1 = ρ in Eq. 8, the easureent of I 1 (ρ 1 I (ρ 1 produces an iage with a resolution gain, with an iaging resolution due to the iage-foring sob-functions, i.e., ρ 1 = ρ µ ρ O. When the lens is large enough to resolve the object, the result is a point-to-point reproduction of the iage only when ρ 1 = ρ ; otherwise for sall lens apertures Eq. 8 fors a point-to-spot iage when ρ O + ρ 1 /µ <λ s O/D and ρ O + ρ /µ <λ s O/D. Now we ove D to the Fourier transfor plane of L B of ar two, i.e., to its focal plane, effectively perforing a Fourier transfor of the field distribution of the iage plane. In addition, D is placed off-center relative to the optic axis of the lens to select part of the spatial frequencies on the Fourier transfor plane, acting as a spatial frequency filter. Matheatically, R c (ρ 1 = I 1 (ρ 1 dρ I (ρ ( dκ F(κ dρ dρ O A(ρ O sob[ π D ρ π D ] sob[ ρ λ s O + ρ ]e ik/(s I (ρ1 ρ (7 e iκ ρ. (9 As a result of the spatial filter function F(κ, the iaging resolution of Eq. 9 is uch narrower than that of a first-order iage; however, it is difficult to siplify this equation further in this for. To get a better understanding of the physics behind Eq. 9, instead of odeling the finite radius of the lens as a disk, which results in the sob-function, we approxiate the finite radius of the lens as the Gaussian function e x L /(D/ with a halfwidth D/, and evaluate in one diension. This leads to a Gaussian iaging-foring function instead of the sob-function. Working in one diension, we change ρ O to x O ; ρ 1 to x 1 ; ρ to x, etc. Then R c (x 1 siplifies to: R c (x 1 dx O A(x O e i k s x kd I 1 ( 4s e (x O +x 1 /µ O kd ( 4s dκ dx F(κ e (x O +x /µ O e i s k x I e iκ x. (10 Corresponding to the experiental easureent, where D was placed at two off-center points in the Fourier transfor plane, we odel the filter function in one diension by two rectangle functions: F(κ =Π(x F x +,D F +Π(x F x,d F, where κ = kx F, while A(x O = Π(x O,w+Π(x O w,w+π(x O + w,w; assuing x, x O, and κ (x F are integrated fro ±, the resulting equation is an analytic expression. Define [ ( kd ( x f(x 1,x O =exp 1 4s + x 1x O µ + 8µ +(kd x O /(4s O ] kx exp[ i( 1 k 3 D 4 x ] O 4µ +(kd /(4s O µs O k D 4 + 18µ s 3 A(x O, O (11 which contains the iaging equation in Gaussian for for x 1 and x, in addition to soe phase ters which are not observable in a first-order iage, and a= µ s O k(k(d/ + iµs O, b(x O= ikµd x O kd + 8iµs O ; (1 5/7
Figure 3. All plots were calculated fro Eq. 13. (a Theoretical coparison of the resolution of a 1.36 lens for three slits with a width and separation of 0.0141. The transparent gold shading shows the first-order unresolved iage, while the light gray indicates the ideal iage; the lines plot Eq. 13 with spatial filter F(κ for several cutoff values κ = kx f B, x ± =±x F, where k is the wavenuber and f B is the focal length of lens L B, versus the position of D 1 in the iage plane. As the cutoff frequency increases the ability to resolve the three slits also increases. (b Coparison of the first-order Gaussian iaging function (black to the second-order iaging function (teal for x F =0.18. The y-axis is in arbitrary units of noralized intensity. (c Theoretical coparison of the resolution of Eq. 13 for a 1.36 lens as the cutoff frequency x F increases. The dot arks the second-order iaging resolution of x F = 0.18 fro (a; the first-order iaging resolution corresponds to x F = 0 Then Eq. 10 is, after evaluating the κ and x integrals, ( R c (x 1 = dx O f(x 1,x O e b(x O (Erfi( b(x O +ak( D F f 4a B x a + ( Erfi( b(x O+ak( D F f B x + a Erfi( b(x O ak( D F f B + x + a Erfi( b(x O ak( D F f B + x a. (13 Then it is easy to see that restricting the allowed spatial frequencies of the Erfi functions constrains the values x O is allowed to take, which, together with f(x 1,x O, iproves the ability to resolve different points on the object plane. However, without evaluation Eq. 13 ay still not be clear enough to show exactly how the resolution is affected, so we have included the following figures which plot soe inforative values to support our experiental observation. Fig. 3(a copares the theoretical first-order unresolved iage of three slits with the second-order fluctuation correlation iage calculated using Eq. 13. The filter function is calculated for a fiber diaeter D F of 0.105 and varying distances x F fro the center of the Fourier transfor plane. Note that the plot for x F = 0.18 deonstrates siilar behavior to the observed experiental data, including the shift of the left and right peaks away fro center. It is clear that, for a lens diaeter of 1.36, the gold transparent plot is copletely unresolved. However, as the cutoff frequency κ = kx F / f B increases, the second-order resolution also increases as seen in the increasing separation of the peaks in Fig. 3. This is ore clearly illustrated in Figs. 3(b and (c; in (b the iaging function of the first-order iage (black is plotted with the second-order iaging function (teal for x F = 0.18. It is clear that the second-order iaging function is uch narrower. In (c the half-width of the iaging function in Eq. 13 is calculated at x 1 =0 for increasing values of x F. Using the estiated experiental paraeters, the second-order iaging resolution starts equal to the first order at x F = 0 and increases to the experiental setup s liit of about 0.00. It is evident fro the experient and theoretical calculations that the increase in spatial resolution is strongly dependent on the chosen spatial filter. It is, in effect, applying a high-pass spatial filter to one ar, producing an edge-sharpening effect. 11,1 The interesting part is that the correlation of the spatially filtered intensity fluctuations with ar one produces a resolved iage, especially since neither ar sees a resolved iage. This correlation filters out the lower spatial frequencies of the unresolved iage of ar one, yielding a resolved iage in the intensity fluctuation correlation rather than the intensity. Conclusion In suary, by using a high-pass spatial filter in the non-resolving side of a two-ar caera, the easureent of the intensity fluctuation correlation I 1 (ρ 1 I (ρ was able to resolve an object that could not be resolved by a traditional caera. This iaging ethod would be particularly useful for long-distance iaging in situations where it is ipractical to have large 6/7
lenses but high resolution is still desired, as it could take advantage of the large angular size of the sun, 0.5 relative to the earth, and the correspondingly sall coherence length, on the order of 0.. In addition, since the theral light iage in I 1 (ρ 1 I (ρ is in general turbulence-free, 14 this ethod would be particularly attractive for satellite caeras taking high resolution iages of objects on the ground. Technically ore coplicated optics or electronics for practical sunlight iaging will be discussed separately. References 1. Ferri, F. et al. High-resolution ghost iage and ghost diffraction experients with theral light. Phys. Rev. Lett. 94, 18360 (005.. Giovannetti, V., Lloyd, S., Maccone, L. & Shapiro, J. H. Sub-Rayleigh-diffraction-bound quantu iaging. Phys. Rev. A 79, 01387 (009. 3. Guerrieri, F. et al. Sub-Rayleigh iaging via n-photon detection. Phys. Rev. Lett. 105, 16360 (010. 4. Gong, W. & Han, S. Multiple-input ghost iaging via sparsity constraints. J. Opt. Soc. A. A 9, 1571 1579 (01. 5. Chen, J., Gong, W. & Han, S. Sub-rayleigh ghost iaging via sparsity constraints based on a digital icro-irror device. Physics Letters A 377, 1844 1847 (013. 6. Oh, J.-E., Cho, Y.-W., Scarcelli, G. & Ki, Y.-H. Sub-Rayleigh iaging via speckle illuination. Opt. Lett. 38, 68 684 (013. 7. Valencia, A., Scarcelli, G., D Angelo, M. & Shih, Y. Two-photon iaging with theral light. Phys. Rev. Lett. 94, 063601 (005. 8. Scarcelli, G., Berardi, V. & Shih, Y. Can two-photon correlation of chaotic light be considered as correlation of intensity fluctuations? Phys. Rev. Lett. 96, 06360 (006. 9. Chen, H., Peng, T. & Shih, Y. 100% correlation of chaotic theral light. Phys. Rev. A 88, 03808 (013. 10. Chen, H. The high-order quantu coherence of theral light. Ph.D. thesis, University of Maryland Baltiore County (014. 11. Hecht, E. Optics (Addison Wesley, 00. 1. Goodan, J. W. Introduction to Fourier Optics (McGraw-Hill, Inc., 1968. 13. Shih, Y. An Introduction to Quantu Optics: Photon and Biphoton Physics (Taylor & Francis, 011. 14. Meyers, R., Deacon, K. & Shih, Y. A new two-photon ghost iaging experient with distortion study. Journal of Modern Optics 54, 381 39 (007. Acknowledgents The authors wish to thank J. Sion and Hui Chen for their helpful discussions. This aterial is based upon work partially supported by the National Science Foundation and the Maryland Innovation Initiative (MII. Author contributions stateent Y.H. conceived the experient(s, J.S. conducted the experient(s, J.S. and T.P. analyzed the results. All authors reviewed the anuscript. Additional Inforation Copeting financial interests The authors declare no copeting financial interests. 7/7