## Thermal light edge enhancement ghost imaging of phase objects

Hanquan Song, Zhidan Yuan, Yuhang Ren, Dayu Zhao, Zili Zhang, Zhiyuan Zheng, and Lu Gao

China University of Geosciences, Beijing 100083, ∗ Corresponding author: gaolu@cugb.edu.cn (Dated: August 6, 2019)

∗ School of Science, China

We propose a experimental scenario of edge enhancement ghost imaging of phase objects with nonlocal orbital angular momentum (OAM) phase filters. Spatially incoherent thermal light is separated into two daughter beams, the test and reference beams, in which the detected objects and phase filters are symmetrically placed, respectively. The results of simulation experiment prove that the edge enhanced ghost images of phase objects can be achieved through the second-order light field intensity correlation measurement owing to the OAM correlation characteristics. Further simulation results demonstrate that the edge enhanced ghost imaging system dose not violate a Bell-type inequality for the OAM subspace, which reveals the classical nature of the thermal light correlation.

Usage:

Manuscript.

PACS numbers: (030.0030) Coherence and statistical optics; (110.0110) Imaging systems; (070.0070) Fourier optics and optical signal processing; (060.5060) Phase modulation.

## I. INTRODUCTION

Ghost imaging is a nonlocal image acquisition technique in which an image can be reconstructed by using a light beam that never interacts with the object through intensity correlation measurement. The phenomenon seems intriguing but can be physically explained by quantum mechanical or classical correlation with a two-photon entangled source or a thermal light source, respectively [1]-[6]. The spatial correlation properties of a quantum entangled source and a thermal light source are di ff erent. The quantum entangled source behaves as a mirror, whereas the classical thermal source looks like a phase-conjugate mirror in the correlated imaging [7]. The debate on the quantum vs classical nature of ghost imaging has lead to other interesting two-photon imaging e ff ects using classical sources [8], e.g., computational ghost imaging [9][10], sub-shot noise ghost imaging of weak absorbing objects [11][12], x-ray region ghost imaging [13][14] and even wavelength transforming ghost imaging of a single photon imaging system [15][16].

Recently, the nature of the correlations between di ff erent OAM components dependent on the phase of the fluctuations of pseudothermal light source have been studied, and the research results show that the thermal light correlation e ff ects share similarities with the quantum correlations in the azimuthal degree of freedom [17][18]. Furthermore, several investigations have addressed two-photon correlation in regard to the high dimensional OAM structure [19]-[21]. Ghost imaging and object identification explorations with OAMquantumcorrelations have been implemented [22]-[25]. Among these studies, B. Jack et al. of the University of Glasgow[23] applies the edge enhancement techniques in classical imaging to ghost imaging and proves that a phase filter, nonlocal with respect to the object, leads to enhanced coincidence images. They also interpret the high-contrast ghost images as a violation of a Bell inequality to demonstrate its quantum nature of the implementation of ghost imaging.

Here in this letter we propose an experimental scenario to accomplish the edge enhancement imaging of phase objects by making use of nonlocal phase filter. The intrinsic of the experiment is owing to the OAM correlation of the incoherent thermal light source. Although there are many studies on the thermal light correlation, we have not yet found the relevant exact proof of its intrinsic to our knowledge. Furthermore, we apply the proposed edge enhance ghost imaging system to measure the Bell-type violation. A circle phase object with π -phase step along the circumference and four phase filter with di ff erent orientations are used to achieve the result of non-violation of the Clauser-Home-Shinomy-Holt Bell-type inequality. The study results provide an explicit proof of the nature of the thermal light ghost imaging system embedded with edge enhancement imaging for phase objects.

## II. RESULTS

The scenario scheme is shown in Fig.1. A He-Ne laser is used to illuminate coherent light beam to a slowly rotation glass( RG ) to produce pseudo-thermal light beam. N is a telescope used to collimate and expand the speckle size of the laser beam, while polarizers P 1 and P 2 can modulate the intensity of the laser beam. The spatially incoherent light beam are generated after the RG and then it is divided into two twin beams, a test beam and a reference beam, by a 50:50 beam splitter ( BS ). The phase object and phase filter are implemented by two S LMs , which are placed in each beam symmetrically. Both S LMs are in the imaging plane of the thermal light source after the RG , and then re-imagined to the transverse planes of the CCDs at the end of each light arm. D 1 and D 2 are connected with a data acquisition card embedded in a personal computer to implement the intensity correlation calculation. The resolution of D 1 and D 2 is 8.3 µ m . Owing to the phase indistinguishability of the correlation characteristics of the incoherent thermal light, images of the detected

FIG. 1. The schematic diagram of the thermal light edge enhancement ghost imaging of phase objects. The phase object is stepped across the transverse plane of the S LM 2 and then imaged onto the CCD . At each position the local form of the object can be treated as a simple phase step of appropriate orientation and correlation measurements are made for di ff erent reference holograms. The details of the experimental schematic are explained in the main text.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01118v1-figure-1d8286a5-beba-4fa5-a52a-893cc9666747.jpg)

phase objects can not be observed from the intensity distribution in the transverse plane of each detector or even from the correlation measurement between the light field intensity of the two detectors. However, a special class of the optical vortices carries OAM, which is characterized by an azimuthal phase dependence of the form exp( il φ ). Here l is the OAM mode number and φ is the azimuthal angle. Recent studies [17][18] have demonstrated the nature of the correlations between di ff erent OAM components dependent on the phase of the fluctuations of pseudothermal light source. Thus the edge enhancement e ff ect of the pure phase object can be implemented through second-order light field intensity correlation measurement. The phase object and phase filter holograms are loaded on the S LMs symmetrically in the test and reference beams, which are shown in the inner pictures of Fig. 1. The size of the phase object is 500 ∗ 500 pixels and the phase filter hologram is 10 ∗ 10 pixels. The red parts of the inner pictures represent the phase of 0, and the blue part of the pictures represent the phase of π . The second-order correlation function between the two separated light beam can be given by ∆ G (2) ( l t , l r ) ∝ | ∫ d φ t exp i [ φ t ( l t -l r )] | 2 , where l t and l r characterize the helical phase modes of the test and reference beams, and φ t presents the phase distribution of the test beam. The second-order correlation exactly happens with l t = l r , which will decrease with the increase of the di ff erence between l t and l r . Any object, or part thereof, can be described by an appropriate superposition of the modes carrying a phase term of exp( il φ ). If the phase object contains the same OAM component as the phase filter, the edge enhancement e ff ect can be accomplished through correlation measurement.

The phase object 'ghost' is loaded on S LM 2 in the test beam, which is shown in Fig.2 (a). And the phase filter with topology charge of l r = 0, l r = 1 or a binary grating with a π -phase step in di ff erent directions is loaded on S LM 1 in the reference beam. The image is acquired by stepping the object in the transverse plane and recording the second-order correlation measurement value. The computational simulation results are shown in Fig. 2.

When a spiral phase filter, with index l r , is placed in the reference arm, the correlation measurement result is proportional to the model component of the object corresponding to l t = l r . The intensity distribution of the test beam is shown in Fig.2 (b) and no information of the object is contained owing to the incoherent characteristics of the thermal light beam. When the phase filter of l r = 0 is placed on S LM 1 , any part of the object with uniform phase of 0 gives a high correlation value owing to l t = l r and the correlation value comes to zero for other part with the non-zero phase distribution. The edge enhancement image is shown in Fig.2 (c), in which the correlation value at the edge of the phase object is lower than other parts with uniform phase distribution of 0. Then the phase filter on S LM 1 is replaced by l r = 1. For a part of the object with a π -phase step, an expansion in terms of exp( il obj φ ) contains non-zero components for l obj = ± 1 [23]. Such a phase step therefore gives a high second-order correlation measurement value for l re f = ± 1. The simulation result is shown in Fig.2 (d), in which the edge of the phase object gives higher secondorder correlation value than other parts of the object. Fig.2 (e) and (f) show the partial edge enhancement ghost images with oriented phase filters of horizontal or vertical π -phase step distributions, respectively. Correlation images appear where the contrast of the edge detection depends on the relative orientation of the edge with respect to the reference phase step.

Bell inequality is a powerful method of discriminating whether a physical phenomenon is of quantum signature or not. Though the thermal light correlation e ff ects share similarities with the quantum correlations both in the spatial and azimuthal degree of freedom, however, the correlated images obtained with thermal light have a finite background. The thermal light correlation characteristics will not violate the Bell inequality, but no explicit proof has been given previously [23]. Here we provide the categorical demonstration through computational simulation based on the thermal light edge enhanced ghost imaging scenario. A circular phase object is chosen to replace the 'ghost' phase object in the test beam and the phase filters in the reference beam have the form of the π -phase step with di ff erent orientations. Both the phase object and the phase filter can be generated by holograms loaded on the S LMs , which are shown in Fig. 1. Owing to the π -phase step distribution of the circular phase object and phase filter, the measurement of the Bell-type inequality is based on the state space for the OAM states l re f = ± 1 with an equally weighted superposition. Only if the phase filter has the same azimuthal orientation as the phase circular object, the enhanced ghost image appears to be 'bright' through thermal light second-order correlation measurement. Since the circular phase object is much larger than the phase filter, the enhanced ghost image is achieved by stepping the phase object in its transverse plane and combining the results of each time correlation measurement. The orientations of the azimuthal

FIG. 2. Simulation results of the thermal light edge enhancement ghost imaging experiment. (a) Phase object. (b) The intensity distribution in the plane of D 2 . (c)-(f) Edge enhanced ghost images of the phase object shown in (a). The phase filter holograms in the reference beam are shown in insets.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01118v1-figure-75754086-0fdf-4541-ad72-f170c02e03dd.jpg)

phase filter in the reference beam are 0, π 2 , π 4 , and 3 π 4 , respectively. We record the thermal light second-order correlation value as a function of the relative angle between the phase step edge of the circular phase object and that of the phase filter in the reference beam. The simulation experimental results are shown in Fig. 3.

Such a Bell inequality measurement can be quantified by making use of a Clauser-Home-Shimony-Holt Bell type inequality with the | S | ≤ 2, and

$$S = E ( \theta _ { A }, \theta _ { B } ) - E ( \theta _ { A }, \theta _ { B } ^ { \prime } ) + E ( \theta _ { A } ^ { \prime }, \theta _ { B } ) + E ( \theta _ { A } ^ { \prime }, \theta _ { B } ^ { \prime } ), \quad \begin{array} { c } \overline { \Theta } \\ \Theta \end{array}$$

where θ A present the orientation of filter and θ B present azimuth angle of the phase object[27]. E ( θ A , θ B ) is calculated from the correlation measurement in particular orientations,

$$E ( \theta _ { A }, \theta _ { B } ) = \frac { C ( \theta _ { A }, \theta _ { B } ) + C ( \theta _ { A } ^ { * }, \theta _ { B } ^ { * } ) - C ( \theta _ { A } ^ { * }, \theta _ { B } ) - C ( \theta _ { A }, \theta _ { B } ^ { * } ) } { C ( \theta _ { A }, \theta _ { B } ) + C ( \theta _ { A } ^ { * }, \theta _ { B } ^ { * } ) + C ( \theta _ { A } ^ { * }, \theta _ { B } ) + C ( \theta _ { A }, \theta _ { B } ^ { * } ) }, \quad \begin{array} { c } \text{hologré} \\ \text{second} \\ \pi \text{ in } ( a ) \\ \text{the cur} \end{array}$$

where θ ∗ ( A B , ) = θ ( A B , ) + π 2 ,the subscripts A and B denote the orientation of filter and azimuth angle of the phase object,and C is the second-order correlation value according to FIG.3(e). Fig. 3(a)-3(d) present images of the circular phase object with the reference hologram oriented at 0, π 2 , π 4 , and 3 π 4 , respectively. Fig. 3(e) shows the variation curves calculated from the enhanced images in each of Fig. 3(a)-3(d). The image data in Fig. 3(e) are averaged from the area unit of 8 radial pixels by 3 azimuthal degrees, and sinusoidal patterns are shown from the second-order correlation measurement. It should be stated that the maximum values of the curves (a) and (b) are larger than those of the (c) and (d). In the condition of (c) and (d) the azimuthal phase filter is inclined of 45 degree, which has the zigzag distribution owing to the square shape of each pixel point of the hologram loaded on the S LM . The zigzag edge phase-step distribution of the phase filter will destroy the value of the second-order correlation measurement.

By submitting θ A = 0, θ B = π 8 , θ ′ A = π 4 and θ ′ B = 3 π 8 into Eqs. (1) and (2), The maximal violation of the Bell inequality can be obtained [26]. In this condition we achieve the maximumvalue of S is 0.57623 through computational simulation, which is less than the local-hidden-variable bound of 2. The simulation result proves that the edge enhanced ghost imaging system is of the classical signature.

FIG. 3. (a)-(d) The simulation results of thermal light edge enhancement ghost imaging of the circular phase object with four di ff erent holograms orientated at 0, π 2 , π 4 ,and 3 π 4 , respectively. (e) present the second-order correlation values with the azimuthal degree from 0 to π in (a)-(d). The dashed lines show the scale in which each point of the curve in (e) are calculated from each local position.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01118v1-figure-8d3bb0b8-6f62-4009-8bfe-4d88aaf9a2a4.jpg)

## IV. CONCLUSION

In summary, we propose a thermal light edge enhanced ghost imaging scheme of phase object. The azimuthal phase

filter is placed non-locally to the phase object, which are in the spatially separated light beams. Isotropic and orientational edge enhanced ghost images of the π -phase step object have been achieved by the phase filter with di ff erent OAM topological charges. The Bell-type inequality has also been measured with a circular phase step object and oriented phase filter. The experimental results show that the maximal violation of the Bell inequality is less than 2, which means the classical signature of the thermal light ghost imaging system. Compared to the conventional optical edge enhanced imaging technology, the proposed scenario here can survive with incoherent light source and the edge enhanced image can even be obtained with non-localized phase filter. It may find potential applications in the fields of medical detection, microscopy imaging and remote sensing.

## ACKNOWLEDGEMENTS

This work was supported by the Fundamental Research Funds for the Central Universities and National Innovation and Entrepreneurship Training Program for College Student.

- [1] T. B. Pittman, Y. H. Shih, D. V . Strekalov, and A. V . Sergienko, Phys. Rev. A 52, R3429 1995.
- [2] D. V. Strekalov, A. V . Sergienko, D. N. Klyshko, and Y . H. Shih, Phys. Rev. Lett. 74, 3600 1995.
- [3] A. Gatti, E. Brambilla, and L. A. Lugiato, Progress in Optics 51, 251 2008.
- [4] R. S. Bennink, S. J. Bentley, and R. W. Boyd, Phys. Rev. Lett. 89, 113601 2002.
- [5] A. Gatti, E. Brambilla, M. Bache, and L. A. Lugiato, Phys. Rev. A 70, 013802 2004.
- [6] J. Cheng and S. S. Han, Phys. Rev. Lett. 92, 093903 2004.
- [7] D. Zh. Cao, J. Xiong, and K. G. Wang, Phys. Rev. A 71, 013801 2005.
- [8] R. Kaltenbaek, J. Lavoie, D. N. Biggersta ff , and K. J. Resch, Nature Phys. 4, 864 (2008).
- [9] J. H. Shapiro, Computational ghost imaging, Phys. Rev. A 78, 061802(R) (2008).
- [10] Y. Altmann, S. McLaughlin, M. J. Padgett, V. K Goyal, A. O. Hero, and D. Faccio, Quantum-inspired computational imag-
- ing, Science 361, eaat2298 (2018).
- [11] G. Brida, M. Genovese, and I. Ruo Berchera, Experimental realization of sub-shot-noise quantum imaging, Nat. Photonics 4, 227 (2010).
- [12] N. Samantaray, I. Ruo Berchera, A. Meda, and M. Genovese, Realization of the first sub-shot-noise wide field microscope, Light Sci. Appl. 6, e17005 (2017).
- [13] D. Pelliccia, A. Rack, M. Scheel, V. Cantelli, and D. M. Paganin, Experimental X-ray Ghost Imaging, Phys. Rev. Lett. 117, 113902 (2016).
- [14] H. Yu, R. Lu, S. Han, H. Xie, G. Du, T. Xiao, and D. Zhu, Fourier-Transform Ghost Imaging with Hard X Rays, Phys. Rev. Lett. 117, 113901 (2016).
- [15] P. A. Morris, R. S. Aspden, J. E. C. Bell, R.W. Boyd, and M. J. Padgett, Imaging with a small number of photons, Nat. Commun. 6, 5913 (2015).
- [16] R. S. Aspden, D. S. Tasca, R.W. Boyd, and M. J. Padgett, EPRbased ghost imaging using a single-photon-sensitive camera, New J. Phys. 15, 073032 (2013).
- [17] O. S. Maga na-Loaiza, M. Mirhosseini, R. M. Cross, S. M. H. Rafsanjani, and R. W. Boyd, Sci. Adv. 2, e1501143 (2016).
- [18] Lu Gao, Seyed Mohammad Hashemi Rafsanjani, Yiyu Zhou, Zhe Yang, Omar S. Maga-Loaiza, Mohammad Mirhosseini, Jiapeng Zhao, Boshen Gao, and Robert W. Boyd, Appl. Phys. Lett. 110, 071107 (2017)
- [19] A. Mair, A. Vaziri, G. Weihs, and A. Zeilinger, Nature (London) 412, 313 (2001).
- [20] S. Franke-Arnold, S. M. Barnett, M. J. Padgett, and L. Allen, Phys. Rev. A 65, 033823 (2002).
- [21] J. P. Torres, A. Alexandrescu, and L. Torner, Phys. Rev. A 68, 050301 (2003).
- [22] N. Uribe-Patarroyo, A. Fraine, D. S. Simon, O. Minaeva, and A. V. Sergienko, Phys. Rev. Lett. 110, 043601 (2013).
- [23] B. Jack, J. Leach, J. Romero, S. Franke-Arnold, M. RitschMarte, S. M. Barnett, and M. J. Padgett, Phys. Rev. Lett. 103, 083602 (2009).
- [24] L. Chen, J. Lei, and J. Romero, Light: Sci. Appl. 3, e153 (2014).
- [25] Z. Yang, O. S. Maga na-Loaiza, M. Mirhosseini, Y. Zhou, B. Gao, L. Gao, S. M. Hashemi Rafsanjani, G. Long, and R. W. Boyd, Light Sci Appl 6, e17013 (2017).
- [26] Leach, J. et al. Violation of a Bell inequality in two-dimensional orbital angular momentum state-spaces. Optics Express 17, 8287-8293 (2009).
- [27] Lebedev, A. V., Lesovik, G. B.,and Blatter, G. Entanglement in a noninteracting mesoscopic structure. Physical Review B 71, 045306 (2005).
Hanquan Song Zhidan Yuan Yuhang Ren Dayu Zhao Zili Zhang Zhiyuan Zheng Lu Gao
Thermal light edge enhancement ghost imaging of phase objects
We propose a experimental scenario of edge enhancement ghost imaging of phase
objects with nonlocal orbital angular momentum (OAM) phase filters. Spatially
incoherent thermal light is separated into two daughter beams, the test and
reference beams, in which the detected objects and phase filters are
symmetrically placed,respectively. The results of simulation experiment prove
that the edge enhanced ghost images of phase objects can be achieved through
the second-order light field intensity correlation measurement owing to the OAM
correlation characteristics. Further simulation results demonstrate that the
edge enhanced ghost imaging system dose not violate a Bell-type inequality for
the OAM subspace, which reveals the classical nature of the thermal light
correlation.