Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Thresholded single-photon underwater imaging and detection

Open Access Open Access

Abstract

Optical underwater target imaging and detection have been a tough but significant challenge in deep-sea exploration. Distant reflected signals drown in various underwater noises due to strong absorption and scattering, resulting in degraded image contrast and reduced detection range. Single-photon feature operating at the fundamental limit of the classical electromagnetic waves can broaden the realm of quantum technologies. Here we experimentally demonstrate a thresholded single-photon imaging and detection scheme to extract photon signals from the noisy underwater environment. We reconstruct the images obtained in a high-loss underwater environment by using photon-limited computational algorithms. Furthermore, we achieve a capability of underwater detection down to 0.8 photons per pulse at Jerlov type III water up to 50 meters, which is equivalent to more than 9 attenuation lengths. The results break the limits of classical underwater imaging and detection and may lead to many quantum-enhanced applications, like air-to-sea target tracking and deep-sea optical exploration.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Over the past few decades, the interest in air-to-sea optical exploration has rised [1,2] because of its potential scalability to different platforms like aircrafts or artificial satellites, and its wide applications in diverse fields, including underwater target detection [3], seafloor topography and marine science. By adding a laser source to the systems, the imaging methods can significantly improve the reduced image contrast under typical underwater conditions [4]. Moreover, the detrimental effects of scattering can be dramatically suppressed by using spatial, temporal or modulation discrimination techniques [5].

Range-gated imaging technology is an effective spatial-temporal discrimination method to overcome the backscattering of water media, which can improve the working distance and the signal-to-noise ratio [611]. The range-gated imaging systems equipped with a high-repetition-frequency laser have gradually become dominant due to its advantages such as small size, low power consumption, and commercially available [1215]. As the distance increases, the signal photons returned from the object become sparse and submerged in scattered photons, which cannot be effectively extracted. After the accurate time-gate operation, most backscattered photons can be filtered out. However, some scattered photons still randomly enter the camera with a certain probability, which is the critical factor that limits the detection range and further improves the image quality.

Recent advances of sensors are transforming target detection to the single-photon regime, especially in low light conditions [1624]. At the same time, various emerging imaging principles and computational algorithms enhance the capabilities of reconstructing objects [2528]. Photon has the unique features of robust to environmental noise, low decoherence properties and low detectable energy level. They play a critical role in quantum computing [2931], quantum simulation [3235], quantum communication [3638], and quantum imaging [3942]. Experiments with photons have been employed to test the foundations of quantum mechanics and the current developments in quantum optics enable the exploration of many new directions. Especially, we witness the recent progresses in underwater quantum communication [4346] and quantum state transmission [47]. However, it has been proven that loss eliminates most of the quantum effects in sensing and imaging, which means that the use of quantum states in these applications is invalid [4851]. Therefore, it is a challenge to realize the potential of quantum technologies and exhibit quantum advantages in underwater conditions.

2. Theory and experiment

Here, we present an experimental demonstration of the quantum-enhanced underwater imaging and detection scheme with the single-photon threshold strategy, which does not detect quantum states directly but separates signals from noises. As is shown in Fig. 1(a), the whole experiment is implemented in a largescale marine test platform with 300 m long, 16 m wide, and 10 m deep, which is the biggest multiple function deep-water towing tank in Asia. The semiopen environment is similar to the natural field situation, which is a benefit to testing the imaging system in a real field condition. The schematic diagram of the experimental setup is shown in Fig. 1(b). The air-to-sea range-gated imaging system contains two sections together located on the same side of the water. We guide photons from the sender section into the water by a cloud stage and collect photons back with the receiver section. A 532 nm pulsed laser beam illuminates the target, provided by a solid-state AOM Q-switched laser operating at a repetition rate of 10 kHz with a tunable average power range of 1-8 W and the full width of half-maximum (FWHM) pulse duration of 50 ns. A small part of the pulse energy ($\sim$ 6 mW) is reflected by the polarization beam splitter onto the Si-PIN photodiode and converted into electrical pulse signals. By comparing the power between the sender and receiver, we can derive the underwater transmission channel loss. The measured attenuation coefficient is about $\alpha$=0.18 $\mathrm {m}^{-1}$, equivalent to an attenuation length of 5.5 m. This attenuation coefficient is also close to the one in coastal seawater Jerlov type III ($\alpha$=0.179 $\mathrm {m}^{-1}$). Our optical setup is a bistatic system in which the separate illuminating and collecting sections can avoid back reflections from the optical components causing damage to the sensitive CCD chip. By the spatial discrimination, a height difference of about 15 cm between the two underwater mirrors can minimize the overlap between the transmitted laser and the field of view of the receiver, which can reduce the effect of backscattering.

 figure: Fig. 1.

Fig. 1. Experimental implementation. A 532 nm pulsed laser beam illuminates the target, which is generated by a Q-switched laser operated at a repetition rate of 10 kHz with the full width of half-maximum pulse duration of 50 ns. A part of the laser output is sent to a linear optical detector to trigger the ICCD camera. The reflected light from the target is gathered by a telescope with an optical aperture of 150 mm. After spectral filtering with a bandpass filter at 532 nm (FWHM $\pm$ 10 nm), photons reflected by the target are focused on the ICCD camera. The focus modulator (FM) is utilized to adjust the size of laser spot by two lens with focal lengths of 30 mm and 75 mm, and the collimation modulator (CM) can guide the laser to illuminate the center of targets. The intensifier is triggered by a Si-PIN detector. (a) Real field-test environment. (b) Sketch of the experimental setup. HWP: half wave plate; QWP: quarter wave plate; PBS: polarization beam splitter; BF: band-pass filter ; ICCD: intensified charge coupled device. (c) Timing diagram of synchronization control signals.

Download Full Size | PDF

Two timing modes of the ICCD camera are used to optimize signal acquisition in this experiment, as is shown in Fig. 1(c). One is the CCD exposure time, during which a frame of image is formed on the CCD chip, here we set the typical exposure time as 0.01 s. The other is the gate delay and width of the intensifier, here we set the gate width as 50 ns, consistent with the full width at half maximum of the laser pulse. Besides, to improve collection efficiency, the gate can be fired many times during each exposure.

There are two critical factors to achieve the precise triggering delay: the periodic pulse signals and the tunable electronic delay mechanism. The photons per pulse are divided into two parts: some are converted into periodic synchronous electrical signals by the Si-PIN device, which triggers the intensifier and labels the initial time; the rest photons make a round-trip through the high-loss underwater channel, then hit the intensifier. The triggering delay estimation is mainly based on the following considerations: the round-trip time in the free space channel and the underwater channel and the inherent insertion delay time in the camera. By tuning the electronic delay, we can determine the precise delay based on the maximum number of photons captured by the ICCD camera during the same exposure.

Three imaging methods are compared based on the whole system and all of them employ a pulsed laser for illumination. The passive method uses an internal camera pulse generator to randomly activate the image intensifier to capture photons, which is often restricted to a range of only a few meters. The active method employs the optimized time-gates synchronized with the pulsed laser to trigger the image intensifier, which can effectively overcome the underwater backscattered photons. As is shown in Fig. 2(a), the ICCD camera periodically switches between “gate on” and “gate off”. The phosphor screen decay on the image intensifier usually follows a multi-exponential function, where the decay time of the phosphor is about 0.1 ms. Therefore, the laser pulse frequency is 10 kHz, which can effectively avoid the fluorescence crosstalk. If the phosphor decay time and the camera time-gate are matched, the photon events can be detected with the maximum intensity shown in Fig. 2(b). In other cases, such as the “late gate”, “early gate” and “no gate” in Fig. 2(c), the intensity of stochastic backscattered photons is lower than that of signal photons. However, the classical intensity integration based on both methods cannot effectively distinguish between signal and backscattered photons.

 figure: Fig. 2.

Fig. 2. Photon threshold strategy. After optimizing the precise activated time-gate, we can obtain the strongest signal photons, while a fraction of the blink responding to backscattered photons appears for the integration period, and therefore the A/D counts are averaged out with the background. (a) Timing of the image intensifier on the ICCD camera. The time intervals of “gate on” and “gate off” have been scaled for better visibility. (b) The integral results corresponding signal photons. (c) The integral results corresponding backscattered photons. “late gate”: the photons hit the intensifier before the activated time gate. “early gate”: the photons hit the intensifier after the activated time gate. “no gate”: the photons hit the intensifier outside the activated time gate. (d) The photon threshold strategy procedure. (i) A sample image including signal photons and noises. The cluster covering an area of 3 pixels $\times$ 3 pixels can avoid the interference of hot pixels. (ii) A distribution of A/D counts on each pixels of the sample image. (iii)-(v) Images obtained for different threshold values of (iii) 560, (iv) 580 and (v) 600. The influence of noise gradually weakens as the threshold value increases, but setting a threshold too high will cause discarding some real photon events. The optimum threshold value in this experiment is 580.

Download Full Size | PDF

Based on the active system, we implement and optimize the photon threshold strategy by the statistics of signals and noises, as is shown in Fig. 2(d). It consists of the differences in A/D counts from individual pixels, including signal photons (i). As the distance increases, the number of returned photons per pulse becomes sparse, so there are only a few signal photons on each frame image. As mentioned before, the response intensity of signal photons on the camera is higher than that of scattered photons and camera noise, and this difference is truly recorded on the camera. To characterize this difference, we collect the A/D counts per pixel from 1017 frames and analyze the distribution of the outputs on each pixel (ii). This distribution can be determined a threshold, counts over which and covering an area of 3 pixels $\times$ 3 pixels are considered as an amplifying electron cluster corresponding to a real signal photon event. Therefore, each original image has a binary outcome. We combine all binary frames to create the final image, the intensity of which is proportional to the probability of successfully hitting the detector from signal photons. Typical integral images for a different threshold of 560, 580 and 600 are shown in (iii)-(v). A low threshold means that more backscattered photons are recognized as signal photons in the final integration, which makes the acquired images noisier. While a high threshold value removes some valuable signals. Therefore, a reasonable threshold is necessary to avoid discarding real signal photons in the final image.

3. Results

We compare the results obtained by the three imaging schemes at the range of 18 m, 24 m, 30 m, 36 m and 42 m respectively, as is shown in Fig. 3. The intensity of the passive and the active methods represents the electronic A/D counts, while the intensity of the thresholded method represents the number of thresholded photons. As the distance increases, the field of view of the target on the camera shrinks. The image quality becomes worse due to the huge spatial geometric loss and the strong backward scattering. The passive method is strongly affected by the backscattering of the laser source and therefore can only image the target in a limited distance (see Fig. 3(a)). However, in the active imaging method, short gate-width can filter out most scattered light noise and get clear images at relatively long distances (see Fig. 3(b)). Interestingly, Fig. 3(c) shows that the photon threshold strategy gives a better performance at the single-photon level and reaches a longer distance. A concentric-ring board is chosen as the test target in this experiment, as is shown in Fig. 3(d), and a laser beam illuminates only the section marked in orange.

 figure: Fig. 3.

Fig. 3. Experimental results for different field-test distances. Three types of imaging strategies are shown in (a)-(c) for different distance (18 m, 24 m, 30 m, 36 m, 42 m). (a) shows that passive method is strongly affected by the backscattering of laser source and can only image the target in a very limited distance. (b) shows that short gate-width is beneficial to filter the scattered light noise and can get clear images at relative long distances. (c) shows that photon threshold strategy gives a better performance in underwater optical imaging. (d) shows the test target and orange section is the imaging area illuminated by the laser beam.

Download Full Size | PDF

In the same conditions, the classical intensity methods integrate all A/D counts in the final image, which is inevitably affected by the multiplication noise. In contrast, the single-photon threshold strategy can register the existence of every single photon. After a succession of exposures, we accumulate these individual binary images into the final image. This methodology means that the signal intensity is proportional to the probability of photons successfully hitting the detector, which is not affected by the multiplication noise. Thus, it is beneficial to extract and optimize weak signals in the underwater environment, compared to a single integral image with the same overall exposure time.

The peak signal-to-noise ratio (PSNR) is usually used to compare image compression quality. To further compare the performance of the above three imaging schemes, we compare the PSNR shown in Table 1. Active intensity imaging can provide a better image quality than passive one; in the long-range, where the ICCD camera collects only a few single photons, active and passive imaging based on intensity accumulation cannot obtain the target information. In the limited-photon condition, the single-photon threshold strategy with area-array pixels has the fast-tracking ability and good robustness to the noisy environment.

Tables Icon

Table 1. The PSNR for all imaging schemes

Since the image data determined by photon threshold strategy is a collection of discrete photons hitting a detector. The measurement results are strongly affected by shot noise due to the inherent nature of coherent light at a low light level. According to prior knowledge, the number of photons in a pulsed laser beam is subject to Poisson statistical distribution [5254]. Combined with the practical measurement process, it can be written the modified format as follows

$$P\left(I_{(i, j)} \mid n_{(i, j)}\right)=\frac{\left(I_{(i, j)}+\beta\right)^{n_{(i, j)}} e^{-\left(I_{(i, j)}+\beta\right)}}{n_{(i, j)} !},$$
where $n_{(i,j)}$ is the number of photons on the ${(i,j)}$ pixel in measured image, and $I_{(i,j)}$ is the modified intensity of each pixel, $\beta$ is an estimated dark-count rate with the typical order of $10^{-4}$ per pixel.

Poisson statistical distribution can be optimized to reconstruct salient features of targets from these noisy measurements as accurately as possible. It can be rewritten as a negative Poisson log-likelihood form

$$\mathcal{F}\left(I_{(i, j)}\right)=\left(I_{(i, j)}+\beta\right)-\sum_{i=1}^{M} \sum_{j=1}^{N} n_{(i, j)} \ln \left(I_{(i, j)}+\beta\right),$$

The term $\ln \left (n_{(i,j)} !\right )$ can be neglected since it is constant with respect to $I_{(i,j)}$, which does not impact the global optimization. The image data is modified within the bounds set by Poisson statistics of the original data. The optimization process for $I_{(i,j)}$ is based on minimizing a convex function, $\Phi$, as follows [55]

$$\begin{array}{ll} \underset{I \in \mathbb{R}^{n}}{\operatorname{minimize}} & \Phi\left(I_{(i, j)}\right) \equiv \mathcal{F}\left(I_{(i, j)}\right)+\tau\left\|I_{(i, j)}\right\|_{\mathrm{TV}}, \\ \textrm{subject to } & I_{(i, j)} \geq 0. \end{array}$$

The first term is the negative Poisson log-likelihood function of the measured image, which facilitates the subsequent optimization process. The second term is a total variation (TV) seminorm regularization, which often converts the image denoising process into a well-posed problem through ensuring the existence and uniqueness of the optimized image. $\tau$ is a balance factor that satisfies the measured data and the sparsity condition constraints. When $\tau$ is at a low level, the reconstructed image is underfitting. In contrast, when $\tau$ is at a very high level, the reconstructed image corresponds to an overfitting result.

We reconstruct the target images of a submarine model and a concentric ring placed in 36 m and 42 m respectively in Fig. 4(a). The original measurements and images for a low value, optimized value and high value of $\tau$ are shown in Fig. 4(b)-(e). These results highlight the trade-off by changing the balance factor $\tau$ between two modified terms in Eq. (3). When $\tau$ is zero, the reconstructed images correspond to the initial measured data in Fig. 4(b). For small values of $\tau$ in Fig. 4(c), the reconstructed images are still noisy in the spatial frequency domain due to undersmoothing. As $\tau$ increases, the images become smoother since the TV seminorm regularization term dominates the final image data, thus the resolution gets degraded for high values of $\tau$. For high values of $\tau$ in Fig. 4(e), the intensity estimate has overfitting distortions. As is shown in Fig. 4(d), the regularization parameters for the central values of $\tau$ give the best-optimized images.

 figure: Fig. 4.

Fig. 4. Regularized image results. (a) Original test objects. (b) Raw measurements captured by the ICCD camera. Reconstructed images for increasing values of $\tau$ in columns (c)-(e). (c) The reconstructed images weighted towards maximizing the log likelihood. (e) The reconstructed images obtained when the optimization algorithm is overly weighted towards increasing the sparsity. (d) The reconstructed images with $\tau$ adjusted to give subjectively the best images.

Download Full Size | PDF

In addition to underwater imaging, single-photon threshold strategy is also sensitive to detect targets by reflected photon statistics. We choose a 2.56 m $\times$ 2.56 m scene area corresponding to 1024 pixels $\times$ 1024 pixels of the ICCD camera with the target about 0.2 m $\times$ 0.2 m at the center of the scene. The target is 50 meters away from the camera, which means the loss is more than 9 attenuation lengths.

We perform the basic calculation about the number of photons per pulse. The maximum energy of a single laser pulse is 300 $\mu$J/pulse. By measuring the power between the sender section and receiver section, we estimate the underwater channel loss is about 40 dB at a distance of 50 meters. The loss wound be double due to the round-trip distance that the laser travels in water. The spatial geometric loss is about 60 dB in the reflection channel. Besides, in the imaging system, the effective quantum efficiency of the ICCD camera is about 10 $\%$, and the reflectivity of mirrors on the cloud stage is about 90 $\%$. Half of the area of targets can return photons. Considering all the above loss, we can calculate about 0.8 photons per pulse on average. In such a photon-limited case, traditional active intensity imaging has been unable to obtain the location information about the target, as is shown in Fig. 5(a). However, the existence of the target will cause the statistical change of the number of reflected photons, which can help us determine whether there are objects by the single-photon threshold strategy.

 figure: Fig. 5.

Fig. 5. The photon threshold strategy for identifying the targets at an ultimately high loss. (a) The reflected photons captured by the ICCD camera at a distance of 50 m away from the scene. (a) The camera can cover an area of 2.56 m $\times$ 2.56 m, a 0.2 m $\times$ 0.2 m target is placed in a center of the scene. After thousands of frames acquisition, signals and background noise in different areas are analyzed and shown in (c)-(f). The number of reflected photons are slightly higher than the surrounding areas. Based on above result, we can lock the location of targets even in a very high-loss scenario.

Download Full Size | PDF

To illustrate this capability, we divide 400 pixels $\times$ 400 pixels into 25 subareas on average, where the target covers an area of 80 pixels $\times$ 80 pixels in the center, as is shown in Fig. 5(b). We add up the total number of photons in these subareas for statistics, respectively. When the number of photons is sparse, the signal is submerged in the fluctuations of noise. The peak area is random so that it is difficult to determine the location of a target shown in Fig. 5(c)-(d). As frames increase, the statistical results change in each subarea as is shown in Fig. 5(e)-(f). Signal photons exceed the fluctuation of environmental noise, and the peak area is locked in the center and no longer changes. Surprisingly, an ultimately high loss may induce a failure of imaging, it is still possible to identify the existence of the target by photon statistics. The photon threshold strategy ensures that the target information can still be obtained even though the signal photons are incredibly sparse, where the classical methods fail to meet the challenge.

4. Conclusion

In summary, we experimentally demonstrate quantum-enhanced underwater imaging and detection using the single-photon threshold method. This strategy does not detect quantum states directly but filters out low counts representing backscattered photons, which is a unique way to break through the classical limits and exhibit quantum-enhanced properties. We verify the unique advantage by analyzing the enhancement on the PSNR of underwater sensing, imaging and other light-detection applications. The method always improves the PSNR in principle in the high-noise regime and can transform underwater optical imaging to the single-photon regime. Besides, we show that the statistical change of the number of reflected photons can identify the existence of targets. The scheme of thresholded single-photon underwater imaging and detection, together with our field test in the longest-ever distance, suggests a promising route for air-to-sea and deep-sea exploration.

Funding

National Key Research and Development Program of China (2017YFA0303700, 2019YFA0308700, 2019YFA0706302); National Natural Science Foundation of China (11690033, 11761141014, 11904229, 61734005); Shanghai Municipal Science and Technology Major Project (2019SHZDZX01); Shanghai Municipal Education Commission (2017-01-07-00-02-E00049); Science and Technology Commission of Shanghai Municipality (17JC1400403).

Acknowledgments

The authors thank Jian-Wei Pan for helpful discussions. X.-M. J. acknowledges additional support from a Shanghai Talent Program and support from Zhiyuan Innovative Research Center of Shanghai Jiao Tong University.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. J. S. Jaffe, K. D. Moore, J. McLean, and M. P. Strand, “Underwater optical imaging: status and prospects,” Oceanography 14(3), 64–75 (2001). [CrossRef]  

2. J. S. Jaffe, “Underwater optical imaging: the past, the present, and the prospects,” IEEE J. Oceanic Eng. 40(3), 683–700 (2015). [CrossRef]  

3. S. Reed, Y. Petillot, and J. Bell, “An automatic approach to the detection and extraction of mine features in sidescan sonar,” IEEE J. Oceanic Eng. 28(1), 90–105 (2003). [CrossRef]  

4. R. D. Richmond and S. C. Cain, Direct-detection LADAR systems, vol. 1 (SPIE, 2010).

5. B. Cochenour, S. P. O’Connor, and L. J. Mullen, “Suppression of forward-scattered light using high-frequency intensity modulation,” Opt. Eng. 53(5), 051406 (2013). [CrossRef]  

6. D. M. Kocak, F. R. Dalgleish, F. M. Caimi, and Y. Y. Schechner, “A focus on recent developments and trends in underwater imaging,” Mar. Technol. Soc. J. 42(1), 52–67 (2008). [CrossRef]  

7. G. R. Fournier, D. Bonnier, J. L. Forand, and P. W. Pace, “Range-gated underwater laser imaging system,” Opt. Eng. 32(9), 2185–2191 (1993). [CrossRef]  

8. G. Fournier, J. Forand, and P. Mathieu, “Range-gated active underwater imaging: Evolution, performance and perspectives,” in C. Proc. Ocean Sciences Meeting’s Presentation, Conference report, (2008), pp. 1–29.

9. B. A. Swartz, “Laser range gate underwater imaging advances,” in Proceedings of OCEANS’94, vol. 2 (IEEE, 1994), pp. II–722.

10. E. McLean, H. Burris, and M. Strand, “Short-pulse range-gated optical imaging in turbid water,” Appl. Opt. 34(21), 4343–4351 (1995). [CrossRef]  

11. K.-L. Chu, Y.-W. Huang, W.-Q. Jin, F.-M. Cao, and H.-L. Li, “System and experimental study of underwater laser range-gated imaging,” Infrared and Laser Engineering, Papers 39, 213–217 (2010).

12. L. Wu, Y. Zhao, Y. Zhang, C. Jin, and J. Wu, “Multipulse gate-delayed range gating imaging lidar,” Opt. Lett. 36(8), 1365–1367 (2011). [CrossRef]  

13. J. Busck and H. Heiselberg, “Gated viewing and high-accuracy three-dimensional laser radar,” Appl. Opt. 43(24), 4705–4710 (2004). [CrossRef]  

14. J. Busck, “Underwater 3-d optical imaging with a gated viewing laser radar,” Opt. Eng. 44(11), 116001 (2005). [CrossRef]  

15. X.-W. Wang, Y. Zhou, F. Song-Tao, J. He, and Y.-L. Liu, “Four-dimensional flash trajectory imaging using time-delay-modulated range-gated viewing,” Opt. Lett. 36(3), 364–366 (2011). [CrossRef]  

16. A. Maccarone, F. M. Della Rocca, A. McCarthy, R. Henderson, and G. S. Buller, “Three-dimensional imaging of stationary and moving targets in turbid underwater environments using a single-photon detector array,” Opt. Express 27(20), 28437–28456 (2019). [CrossRef]  

17. A. Maccarone, A. McCarthy, X.-M. Ren, R. E. Warburton, A. M. Wallace, J. Moffat, Y. Petillot, and G. S. Buller, “Underwater depth imaging using time-correlated single-photon counting,” Opt. Express 23(26), 33911–33926 (2015). [CrossRef]  

18. A. Maccarone, A. McCarthy, J. Tachella, D. A. Garcia, F. M. Della Rocca, Y. Altmann, S. McLaughlin, R. Henderson, and G. S. Buller, “Underwater depth imaging using time-correlated single-photon counting at video frame rates,” in Electro-Optical Remote Sensing XIII, vol. 11160 (International Society for Optics and Photonics, 2019), pp. 111600I.

19. M. Aßmann, F. Veit, M. Bayer, M. van der Poel, and J. M. Hvam, “Higher-order photon bunching in a semiconductor microcavity,” Science 325(5938), 297–300 (2009). [CrossRef]  

20. J. Wiersig, C. Gies, F. Jahnke, M. Aßmann, T. Berstermann, M. Bayer, C. Kistner, S. Reitzenstein, C. Schneider, S. Höfling, and A. Forchel, “Direct observation of correlations between individual photon emission events of a microcavity laser,” Nature 460(7252), 245–249 (2009). [CrossRef]  

21. M. Aßmann, F. Veit, M. Bayer, C. Gies, F. Jahnke, S. Reitzenstein, S. Höfling, L. Worschech, and A. Forchel, “Ultrafast tracking of second-order photon correlations in the emission of quantum-dot microresonator lasers,” Phys. Rev. B 81(16), 165314 (2010). [CrossRef]  

22. O. Schwartz, J. M. Levitt, R. Tenne, S. Itzhakov, Z. Deutsch, and D. Oron, “Superresolution microscopy with quantum emitters,” Nano Lett. 13(12), 5832–5836 (2013). [CrossRef]  

23. M. P. Edgar, D. S. Tasca, F. Izdebski, R. E. Warburton, J. Leach, M. Agnew, G. S. Buller, R. W. Boyd, and M. J. Padgett, “Imaging high-dimensional spatial entanglement with a camera,” Nat. Commun. 3(1), 984 (2012). [CrossRef]  

24. R. Chrapkiewicz, M. Jachura, K. Banaszek, and W. Wasilewski, “Hologram of a single photon,” Nat. Photonics 10(9), 576–579 (2016). [CrossRef]  

25. C.-S. Yang, Y.-H. Yao, C.-Z. Jin, D.-L. Qi, F.-Y. Cao, Y.-L. He, J.-L. Yao, P.-P. Ding, L. Gao, T.-Q. Jia, J.-Y. Liang, Z.-R. Sun, and S.-A. Zhang, “High-fidelity image reconstruction for compressed ultrafast photography via an augmented-lagrangian and deep-learning hybrid algorithm,” Photonics Res. 9(2), B30–B37 (2021). [CrossRef]  

26. N. Bhusal, S. Lohani, C.-L. You, M.-Y. Hong, J. Fabre, P. Zhao, E. M. Knutson, R. T. Glasser, and O. S. Magaña-Loaiza, “Spatial mode correction of single photons using machine learning,” Adv. Quantum Technol. 4(3), 2000103 (2021). [CrossRef]  

27. Z.-M. Li, S.-B. Wu, J. Gao, H. Zhou, Z.-Q. Yan, R.-J. Ren, S.-Y. Yin, and X.-M. Jin, “Fast correlated-photon imaging enhanced by deep learning,” Optica 8(3), 323–328 (2021). [CrossRef]  

28. P. Wang, J.-Y. Liang, and L.-H. V. Wang, “Single-shot ultrafast imaging attaining 70 trillion frames per second,” Nat. Commun. 11, 2091 (2020). [CrossRef]  

29. M. A. Broome, A. Fedrizzi, S. Rahimi-Keshari, J. Dove, S. Aaronson, T. C. Ralph, and A. G. White, “Photonic boson sampling in a tunable circuit,” Science 339(6121), 794–798 (2013). [CrossRef]  

30. M. Bentivegna, N. Spagnolo, C. Vitelli, F. Flamini, N. Viggianiello, L. Latmiral, P. Mataloni, D. J. Brod, E. F. Galvão, A. Crespi, R. Ramponi, R. Osellame, and F. Sciarrino, “Experimental scattershot boson sampling,” Sci. Adv. 1(3), e1400255 (2015). [CrossRef]  

31. K. Sun, J. Gao, M.-M. Cao, Z.-Q. Jiao, Y. Liu, Z.-M. Li, E. Poem, A. Eckstein, R.-J. Ren, X.-L. Pang, I. A. Walmsley, H. Tang, and X.-M. Jin, “Mapping and measuring large-scale photonic correlation with single-photon imaging,” Optica 6(3), 244–249 (2019). [CrossRef]  

32. A. Peruzzo, M. Lobino, J. C. F. Matthews, N. Matsuda, A. Politi, K. Poulios, X.-Q. Zhou, Y. Lahini, N. Ismail, K. Wörhoff, Y. Bromberg, Y. Silberberg, M. G. Thompson, and J. L. OBrien, “Quantum walks of correlated photons,” Science 329(5998), 1500–1503 (2010). [CrossRef]  

33. K. Poulios, R. Keil, D. Fry, J. D. A. Meinecke, J. C. F. Matthews, A. Politi, M. Lobino, M. Gräfe, M. Heinrich, S. Nolte, A. Szameit, and J. L. O’Brien, “Quantum walks of correlated photon pairs in two-dimensional waveguide arrays,” Phys. Rev. Lett. 112(14), 143604 (2014). [CrossRef]  

34. J. O. Owens, M. A. Broome, D. N. Biggerstaff, M. E. Goggin, A. Fedrizzi, T. Linjordet, M. Ams, G. D. Marshall, J. Twamley, M. J. Withford, and A. G. White, “Two-photon quantum walks in an elliptical direct-write waveguide array,” New J. Phys. 13(7), 075003 (2011). [CrossRef]  

35. J. Gao, L.-F. Qiao, X.-F. Lin, Z.-Q. Jiao, Z. Feng, Z. Zhou, Z.-W. Gao, X.-Y. Xu, Y. Chen, H. Tang, and X.-M. Jin, “Non-classical photon correlation in a two-dimensional photonic lattice,” Opt. Express 24(12), 12607–12616 (2016). [CrossRef]  

36. Z.-Q. Yan, C.-Q. Hu, Z.-M. Li, Z.-Y. Li, H. Zheng, and X.-M. Jin, “Photon-inter-correlation optical communication,” arXiv preprint arXiv:2104.08913 (2021).

37. C.-Z. Peng, T. Yang, X.-H. Bao, J. Zhang, X.-M. Jin, F.-Y. Feng, B. Yang, J. Yang, J. Yin, Q. Zhang, N. Li, B.-L. Tian, and J.-W. Pan, “Experimental free-space distribution of entangled photon pairs over 13 km: towards satellite-based global quantum communication,” Phys. Rev. Lett. 94(15), 150501 (2005). [CrossRef]  

38. X.-M. Jin, J.-G. Ren, B. Yang, Z.-H. Yi, F. Zhou, X.-F. Xu, S.-K. Wang, D. Yang, Y.-F. Hu, S. Jiang, T. Yang, H. Yin, C.-Z. Peng, and J.-W. Pan, “Experimental free-space quantum teleportation,” Nat. Photonics 4(6), 376–381 (2010). [CrossRef]  

39. O. S. Magaña-Loaiza and R. W. Boyd, “Quantum imaging and information,” Rep. Prog. Phys. 82(12), 124401 (2019). [CrossRef]  

40. P.-A. Moreau, E. Toninelli, T. Gregory, and M. J. Padgett, “Ghost imaging using optical correlations,” Laser Photonics Rev. 12(1), 1700143 (2018). [CrossRef]  

41. P.-A. Moreau, E. Toninelli, T. Gregory, and M. J. Padgett, “Imaging with quantum states of light,” Nat. Rev. Phys. 1(6), 367–380 (2019). [CrossRef]  

42. M. Gilaberte Basset, F. Setzpfandt, F. Steinlechner, E. Beckert, T. Pertsch, and M. Gräfe, “Perspectives for applications of quantum imaging,” Laser Photonics Rev. 13(10), 1900097 (2019). [CrossRef]  

43. L. Ji, J. Gao, A.-L. Yang, Z. Feng, X.-F. Lin, Z.-G. Li, and X.-M. Jin, “Towards quantum communications in free-space seawater,” Opt. Express 25(17), 19795–19806 (2017). [CrossRef]  

44. F. Bouchard, A. Sit, F. Hufnagel, A. Abbas, Y. Zhang, K. Heshami, R. Fickler, C. Marquardt, G. Leuchs, R. W. Boyd, and E. Karimi, “Underwater quantum key distribution in outdoor conditions with twisted photons,” arXiv preprint arXiv:1801.10299 (2018).

45. C.-Q. Hu, Z.-Q. Yan, J. Gao, Z.-Q. Jiao, Z.-M. Li, W.-G. Shen, Y. Chen, R.-J. Ren, L.-F. Qiao, A.-L. Yang, H. Tang, and X.-M. Jin, “Transmission of photonic polarization states through 55-m water: towards air-to-sea quantum communication,” Photonics Res. 7(8), A40–A44 (2019). [CrossRef]  

46. C.-Q. Hu, Z.-Q. Yan, J. Gao, Z.-M. Li, H. Zhou, J.-P. Dou, and X.-M. Jin, “Decoy-state quantum key distribution over a long-distance high-loss air-water channel,” Phys. Rev. Appl. 15(2), 024060 (2021). [CrossRef]  

47. Y. Chen, W.-G. Shen, Z.-M. Li, C.-Q. Hu, Z.-Q. Yan, Z.-Q. Jiao, J. Gao, M.-M. Cao, K. Sun, and X.-M. Jin, “Underwater transmission of high-dimensional twisted photons over 55 meters,” PhotoniX 1(1), 5–11 (2020). [CrossRef]  

48. T.-W. Lee, S. D. Huver, H. Lee, L. Kaplan, S. B. McCracken, C. Min, D. B. Uskov, C. F. Wildfeuer, G. Veronis, and J. P. Dowling, “Optimization of quantum interferometric metrological sensors in the presence of photon loss,” Phys. Rev. A 80(6), 063803 (2009). [CrossRef]  

49. U. Dorner, R. Demkowicz-Dobrzanski, B. Smith, J. Lundeen, W. Wasilewski, K. Banaszek, and I. Walmsley, “Optimal quantum phase estimation,” Phys. Rev. Lett. 102(4), 040403 (2009). [CrossRef]  

50. L. Cohen, E. S. Matekole, Y. Sher, D. Istrati, H. S. Eisenberg, and J. P. Dowling, “Thresholded quantum lidar: Exploiting photon-number-resolving detection,” Phys. Rev. Lett. 123(20), 203601 (2019). [CrossRef]  

51. Z. Bao, Y. Liang, Z. Wang, Z. Li, E. Wu, G. Wu, and H. Zeng, “Laser ranging at few-photon level by photon-number-resolving detection,” Appl. Opt. 53(18), 3908–3912 (2014). [CrossRef]  

52. M. Fox, Quantum optics: an introduction, vol. 15 (Oxford University, 2006).

53. H.-A. Bachor and T. C. Ralph, A guide to experiments in quantum optics (John Wiley & Sons, 2019).

54. M. O. Scully and M. S. Zubairy, Quantum Optics (Cambridge University, 1997).

55. Z. T. Harmany, R. F. Marcia, and R. M. Willett, “This is SPIRAL-TAP: Sparse Poisson Intensity Reconstruction ALgorithms—Theory and Practice,” IEEE Trans. on Image Process. 21(3), 1084–1096 (2012). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1.
Fig. 1. Experimental implementation. A 532 nm pulsed laser beam illuminates the target, which is generated by a Q-switched laser operated at a repetition rate of 10 kHz with the full width of half-maximum pulse duration of 50 ns. A part of the laser output is sent to a linear optical detector to trigger the ICCD camera. The reflected light from the target is gathered by a telescope with an optical aperture of 150 mm. After spectral filtering with a bandpass filter at 532 nm (FWHM $\pm$ 10 nm), photons reflected by the target are focused on the ICCD camera. The focus modulator (FM) is utilized to adjust the size of laser spot by two lens with focal lengths of 30 mm and 75 mm, and the collimation modulator (CM) can guide the laser to illuminate the center of targets. The intensifier is triggered by a Si-PIN detector. (a) Real field-test environment. (b) Sketch of the experimental setup. HWP: half wave plate; QWP: quarter wave plate; PBS: polarization beam splitter; BF: band-pass filter ; ICCD: intensified charge coupled device. (c) Timing diagram of synchronization control signals.
Fig. 2.
Fig. 2. Photon threshold strategy. After optimizing the precise activated time-gate, we can obtain the strongest signal photons, while a fraction of the blink responding to backscattered photons appears for the integration period, and therefore the A/D counts are averaged out with the background. (a) Timing of the image intensifier on the ICCD camera. The time intervals of “gate on” and “gate off” have been scaled for better visibility. (b) The integral results corresponding signal photons. (c) The integral results corresponding backscattered photons. “late gate”: the photons hit the intensifier before the activated time gate. “early gate”: the photons hit the intensifier after the activated time gate. “no gate”: the photons hit the intensifier outside the activated time gate. (d) The photon threshold strategy procedure. (i) A sample image including signal photons and noises. The cluster covering an area of 3 pixels $\times$ 3 pixels can avoid the interference of hot pixels. (ii) A distribution of A/D counts on each pixels of the sample image. (iii)-(v) Images obtained for different threshold values of (iii) 560, (iv) 580 and (v) 600. The influence of noise gradually weakens as the threshold value increases, but setting a threshold too high will cause discarding some real photon events. The optimum threshold value in this experiment is 580.
Fig. 3.
Fig. 3. Experimental results for different field-test distances. Three types of imaging strategies are shown in (a)-(c) for different distance (18 m, 24 m, 30 m, 36 m, 42 m). (a) shows that passive method is strongly affected by the backscattering of laser source and can only image the target in a very limited distance. (b) shows that short gate-width is beneficial to filter the scattered light noise and can get clear images at relative long distances. (c) shows that photon threshold strategy gives a better performance in underwater optical imaging. (d) shows the test target and orange section is the imaging area illuminated by the laser beam.
Fig. 4.
Fig. 4. Regularized image results. (a) Original test objects. (b) Raw measurements captured by the ICCD camera. Reconstructed images for increasing values of $\tau$ in columns (c)-(e). (c) The reconstructed images weighted towards maximizing the log likelihood. (e) The reconstructed images obtained when the optimization algorithm is overly weighted towards increasing the sparsity. (d) The reconstructed images with $\tau$ adjusted to give subjectively the best images.
Fig. 5.
Fig. 5. The photon threshold strategy for identifying the targets at an ultimately high loss. (a) The reflected photons captured by the ICCD camera at a distance of 50 m away from the scene. (a) The camera can cover an area of 2.56 m $\times$ 2.56 m, a 0.2 m $\times$ 0.2 m target is placed in a center of the scene. After thousands of frames acquisition, signals and background noise in different areas are analyzed and shown in (c)-(f). The number of reflected photons are slightly higher than the surrounding areas. Based on above result, we can lock the location of targets even in a very high-loss scenario.

Tables (1)

Tables Icon

Table 1. The PSNR for all imaging schemes

Equations (3)

Equations on this page are rendered with MathJax. Learn more.

P ( I ( i , j ) n ( i , j ) ) = ( I ( i , j ) + β ) n ( i , j ) e ( I ( i , j ) + β ) n ( i , j ) ! ,
F ( I ( i , j ) ) = ( I ( i , j ) + β ) i = 1 M j = 1 N n ( i , j ) ln ( I ( i , j ) + β ) ,
minimize I R n Φ ( I ( i , j ) ) F ( I ( i , j ) ) + τ I ( i , j ) T V , subject to  I ( i , j ) 0.
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.