## Frame-filtered ghost imaging with a SPAD array used both as a multiple "bucket" detector and an imaging camera

V. S. Starovoitov, V. N. Chizhevsky, D. Mogilevtsev, A. Smaliakou, and S. Kilin B.I.Stepanov Institute of Physics, NAS of Belarus, Nezavisimosti ave. 68, 220072 Minsk, Belarus

## M. Perenzoni

Fondazione Bruno Kessler FBK, 38122 Trento, Italy and now at Sony Semiconductor Solutions Europe, Trento, Italy

## L. Gasparini

Fondazione Bruno Kessler FBK, 38122 Trento, Italy

## D. B. Horoshko

Univ. Lille, PhLAM - Physique des Lasers Atomes et Molecules, F-59000 Lille, France (Dated: April 2024)

An approach to ghost imaging with a single SPAD array used simultaneously as a several-pixel "bucket" detector and an imaging camera is described. The key points of the approach are filtering data frames used for ghost-image reconstruction by the number of per-frame counts and superposing correlation images obtained for different "bucket" pixels. The imaging is performed in an experiment with a pseudo-thermal light source where the light intensity is so low that the dark counts have a noticeable effect on imaging. We demonstrate that the approach is capable to significantly reduce the destructive effect of dark counts on the ghost image and improve image contrast, spatial resolution, and image similarity to a reference image.

## I. INTRODUCTION

Ghost imaging (GI) is already a mature research field dating back nearly a quarter of century. Since the first demonstration for GI with twin-photons [1] and subsequent GI realization with classically correlated light [24], a large corpus of experimental and theoretical works has appeared addressing particular aspects of GI and its applications for imaging/sensing (see, for example, review works [5-10]). The major practical appeal of GI is a possibility to use for image detection a maximally simple bucket detector unable to produce an image by itself. One has to use high-resolution detection only with the reference beam (a part of the correlated light not interacting with the object). Moreover, even correlations of the bucket-detector photo-current with the reference current are able to provide data sufficient for building a correlation image [11] (this possibility gave rise to a rapidly developing "single-pixel" imaging technique with such possible applications as detection of gas leaks and sensors for situational awareness for vehicles [12]).

ple, in near-infrared fluorescence lifetime imaging [16], super-resolving quantum imaging [17, 18], demonstrating of EPR inequality violation by measurement of photon correlations [19], correlation plenoptic imaging [20], Hong-Ou-Mandel interference microscopy [21], LIDAR applications [22]. In these applications, the image is reconstructed from data acquired in fixed temporal windows, hereafter referred to as frames.

One of the most perspective GI ways to enhance imaging sensitivity and account for non-classical light correlations is to apply recently developed single photon avalanche diode (SPAD) arrays for the reference-beam detection [7]. Currently, the SPAD arrays operated in various detection regimes (such as time-stamping, counting, and gating) have reached a high technological maturity [13], featuring observation rates up to 1 MHz, 200 ps time resolution [14], close to 100 % fill factor with megapixel spatial resolution [15]. SPAD cameras have already found a widespread use in imaging, for exam-

The dark counts (DC) are an inevitable effect that limits the sensitivity of SPAD-based imaging. The DC are caused by intrinsic generation of carriers within the detector in absence of illumination. The DC are not correlated with the photocounts, but they, however, cannot be separated from them. The higher the rate of DC in a pixel, the lower the correlation with the light for the events it detects. In large-size detector arrays, a significant percentage of detector-pixels (such pixels are called 'hot' pixels) could show a DC rate much higher than the median DC rate of the array, sometimes even exceeding the typical detected signal. In correlation-based imaging applications, the destructive action of DC cannot be reduced by data averaging over time. Eventually, the effect of dark counts on measurement may be minimized by disabling noisy pixels [23-25].

In this work we use the SPAD array for a ghost imaging modality that makes a step further toward full-size quantum correlation imaging with aim to improve image quality. Namely, instead of a single-pixel bucket detector we use several pixels of the imaging SPAD camera (as in [4, 26]) and perform averaging over a ghost images obtained for each "bucket" pixel. Also, we propose procedures allowing to deal with intrinsic shortcomings of SPADs in the low-light regime. The procedures are

FIG. 1: The scheme of the ghost imaging experiment: ( LD ) laser diode; ( LDM ) laser diode mount; C LD fiber collimator; ( F BB ) broadband filter; ( R B ) beam reducer; ( D ) rotating ground-glass disk; ( L C , L T ) lenses; ( D I ) iris diaphragm; ( S 1 , S 2 ) beam splitters; ( R Ref , R T ) reflectors; ( T ) target; ( F NB ) narrow-band filter; ( A SPAD ) SPAD array; ( Cam ) camera "SuperEllen"; ( PC ) computer. The inset shows how the SPAD array is illuminated by the object and reference beams. The dashed line in the inset shows boundaries of the 'Flare' area; the "bucket" pixels are there.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01153v1-figure-a56ca2af-812f-47e7-a3d8-343d8d44b278.jpg)

frame-filtering on the number of counts per frame and algorithms for dealing with dark counts in-homogeneously distributed over the SPAD pixels. We demonstrate that for a standard pseudo-thermal source of correlated photons (namely, a rotating ground-glass disk) such an approach can significantly rise the ghost-image quality.

It should be emphasized that key points of our approach, namely, filtering frames with respect to a specific number of counts and averaging over several simultaneously obtained intensity correlation images look to be rather general and applicable not only for different ghost imaging modalities, but also for generic correlation imaging with light on the few-photon level, first of all, for quantum microscopy/sensing applications.

The outline of the paper is as follows. In the second Section we describe the ghost imaging set-up, the object and the SPAD camera used for imaging. In the third Section we outline the data post-processing procedure, give the definition of the ghost image for the case of multiple bucket detectors and define two measures of image quality: the contrast-to-noise ratio and the correlation coefficient to the reference image. In the fourth Section the ghost image inference and quantification of the image contrast and quality are discussed. Here we demonstrate ghost imaging with filtered and unfiltered data, and also compare results averaged over multiple "bucket" pixels and just a single "bucket" pixel. In the Appendix, the influence of dark counts on the ghost image is discussed.

## II. SETUP

## A. Optical layout

In the work, the GI is performed by an experimental way using correlated light beams produced with help of a laser and a rotating ground-glass disk. The initial light beam is generated by a CW-operated laser source. The experimental setup is shown in Fig 1. The source is a fiber-pigtailed single-transverse-mode FabryPerot laser diode (LDS-670-FP-5-U-3-SMP04-FU-CW0.5, LaserCom LLC) operated near a wavelength λ ≈ 671 nm with an output power less than 8 mW. The diode is thermally stabilized inside a laser mount (Thorlabs LDM9T) and operates at constant current in nearthreshold mode (power supply Thorlabs LDC201CU). The laser beam is output from a PM-fiber and a fiber collimator ( C LD , efficient focus length of 7.5 mm) mounted in a Thorlabs PAF2-7B fiber-port. This beam is collimated. Its cross section has a Gaussian profile with a diameter of approximately 1.4 mm. The beam is directed through a broadband filter F BB and a free-space beam reducer ( R B , of about 3 times reduction factor) towards a rotating ground-glass disk (D). The distances | C LD D | and | R D B | are about 100 mm and 20 mm, respectively. The diameter of spot illuminated by this beam on the rotating disk is d d ≈ 0 47 . mm. Disc D has a matte surface with a random grained transparent pattern that scatters the input optical beam over a wide solid angle. The luminous laser spot on the surface of disk D is a pseudo-thermal source of correlated photons that form random structure of light spots (also called speckles) in space. The distance between the optical axis and the axis of disk rotation is 52 mm. The disk rotates at a frequency of 10 Hz. A lens ( L C , focal length F LC = 80 mm) positioned at a distance | DL C | ≈ F LC collects photons scattered by the disk, collimates and directs them through an iris diaphragm ( D I , | DD I | ≈ 20 mm) towards a beam splitter ( S 1 , | DS 1 | ≈ 55 mm).

Splitter S 1 divides the collimated speckle beam into two arms: object and reference. In the reference arm, the beam is directed through a reflector ( R Ref , | S R 1 Ref | ≈ 105 mm) and another beam splitter ( S 2 , | R S T 2 | ≈ 85 mm) towards a SPAD array ( A SPAD ) mounted in a camera (Cam). In the object arm, the beam is directed by a reflector ( R T , | S R 1 T | ≈ 85 mm) to a target (T). A lens ( L T , focus length F LT = 67 mm) collects the light transmitted by target T and focuses it (through beam splitter S 2 ) on array A SPAD ( | L A T SPAD | ≈ F LT ). Target T is positioned as close as possible to splitter S 2 . The distances | TL T | and | L S T 2 | are 13 and 15 mm, respectively. A 20 nm narrowband filter ( F NB ) is located in front of the camera for suppressing influence of stray light. The distance | S A 2 SPAD | ≈ 50 mm. The beam splitters S 1 and S 2 are 50:50 non-polarizing beamsplitter cubes BS011 (Thorlabs). The reflectors R Ref and R T are Thorlabs right-angle prism FS910L-B.

Splitter S 2 directs photons from both the reference and

FIG. 2: Resolution test target (left) and high-resolution contour plot of the object with sizes (right): a ≈ 1 3 . mm; b ≈ 0 94 . mm; c ≈ 0 20 . mm; d ≈ 0 16 . mm; e ≈ 0 116 . mm; f ≈ 0 12 . mm; g ≈ 0 19 . mm.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01153v1-figure-6f127e73-3822-4d1e-aa55-d9c3e4cb5e0f.jpg)

object arms to array A SPAD . The reference beam uniformly fills the entire surface of array A SPAD . The object arm is adjusted so that the light spot focused by lens L T on the surface of array A SPAD falls on the upper right part of A SPAD (see the inset in Fig. 1(a)). The diameter of this spot is d T = d F d LT /F C ≈ 0 39 . mm. The central part of this spot (referred to as 'bucket') is intended to be used as a bucket detector. A larger area (an area limited by dashed lines in the inset in Fig. 1(a)) around the object-beam spot is considered as a flare area, unsuitable for image reconstruction. The remaining part of array A SPAD (called 'Ghost') lighted only by the reference beam is reserved for ghost image reconstruction. A similar arrangement of having a "bucket" and "ghost" images in the region of just a single camera (albeit a CCD one) was used in an earlier work on ghost imaging with quasi-classical field [4]. Recently, a similar set-up was used in works on ghost polarimetry [26-28].

## B. Imaged object and SPAD array

The imaged object is a part of a 1951 USAF resolution test target (negative chrome-on-glass version, see Fig. 2) discussed in detail in [29]. We reproduce the images for the digit template '2' (Element 2 of Group 0) of the test target. In order to illuminate the template '2' selectively from other target patterns, the light beams are restricted by diaphragm D I . The height and width of the selected template are 1.3 and 0.94 mm. The font-line thickness for this object is a value varying from 0.11 to 0.20 mm. The image of the selected digit template is shown in Fig. 2.

The SPAD camera (Cam) is based on the "SuperEllen" sensor, specifically developed by Fondazione Bruno Kessler to target quantum imaging applications [14, 18] and implemented recently also for ghost imaging

[30, 31]. The array consists of 32 × 32 = 1024 pixels with a pitch of δ s = 44 54 . µm and a total sensitive area of 1 4 . × 1 4 . mm 2 manufactured in 150 nm CMOS standard technology. The pixels are addressed by their linear index k ∈ { 1 , . . . , 1024 } such that k = y k + 32( x k -1) where x , y k k ∈ { 1 , . . . , 32 } are the discrete Cartesian coordinates of pixels. For every pixel there is a time-to-digital converter (TDC) which timestamps the first count with up to 204 ps resolution within a frame of exposure time of up to 100 ns. Using on-chip features for reading empty rows or frames, a frequency of exposure repetition (i.e. an observation rate) of 800 kHz is achieved leading to a measurement duty-cycle of 3.6 % . The photon detection efficiency reaches 5% at 400 nm and 0 8% . at 810 nm. We did not specifically measure the efficiency at the used wavelength 671 nm. However, we estimate it as being between 1 5% . and 2% . The median dark count rate per pixel is approximately 0 7 . kHz over the whole pixel population. The timestamped counts are registered in a frame-by-frame way.

In our experiment, each frame is a result of light exposition during the time of 45 ns. As it is natural to expect with our pseudo-thermal light source, there are timecorrelation between frames. We have addressed these correlations in our recent work [32]. However, their presence does not affect current consideration. The data related to an individual frame is converted into a compressed sparse structure

$$\mathcal { D } _ { f } \coloneqq \{ f, n _ { f }, \{ x _ { k } ^ { ( i ) }, y _ { k } ^ { ( i ) }, \tau ^ { ( i ) } \} _ { i = 1, \dots, n _ { f } } \}, \quad ( 1 )$$

where f is the serial number of frame repetition; n f is the total number of events (counts) detected by all pixels in frame f ; x ( ) i k , y ( ) i k and τ ( ) i are the discrete Cartesian coordinates and the detection time for the i -th event registered in frame f . The registered data is collected and transferred to a computer (PC) as a stream of multiple frames F = {D f ∣ ∣ n f ≥ 1 } . The frames {D | f n f = 0 } (frames for which no events were detected) are not included in the stream for accelerating the data transfer process.

## III. POST-PROCESSING

## A. Datasets and correlation functions

For ghost imaging both unfiltered and filtered datasets will be used. The unfiltered dataset F 0 = { D f ∣ ∣ n f ≥ 0 } , contains all frames registered in the experiment (including the frames of original stream F and the frames {D | f n f = 0 } omitted at data transfer. The total number of frames in this dataset is N 0 .

The frame-filtered datasets

$$F ( n ^ { * } ) = F _ { 0 } | _ { n _ { f } = n ^ { * } } = \left \{ \mathcal { D } _ { f } | n _ { f } = n ^ { * } \right \} \quad \ \ ( 2 )$$

are non-overlapping sub-sets of F 0 with fixed number of

FIG. 3: Functions g (2) ( m ,k ∗ ) depending on the distance ∆ y m k ∗ for the pixel m ∗ = 496 ( x m ∗ = y m ∗ = 16 ). The functions are acquired from a frame-filtered dataset F (8) (blue squares) and unfiltered data F 0 (red stars). The blue solid ( F (8) ) and red dashed ( F 0 ) lines indicate approximation (6) to the acquired functions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01153v1-figure-db054490-6c17-42dc-bdad-5027e1fa6ef3.jpg)

counts n ∗ in frame. The ratio

$$p ( n _ { f } = n ^ { * } ) = N ( n ^ { * } ) / N _ { 0 } \quad \quad ( 3 ) \quad \text{pro}$$

is normalised distribution of filtered sub-sets number N n ( ∗ ) .

The prepared dataset is used to calculate both a per-frame distribution of counts and their pairwise correlations[32]:

$$G ^ { ( 1 ) } ( k ) = \langle P _ { k } \rangle _ { \{ f \} }, \ \ G ^ { ( 2 ) } ( m, k ) = \langle P _ { m } P _ { k } \rangle _ { \{ f \} }, \ \ ( 4 ) \quad \text{the p}$$

where k and m are linear pixel indices, and ⟨ . . . ⟩ { f } denotes averaging over the dataset frames. The quantities P k in (4) are random binaries taking the value 1 for a count registration at the pixel k or zero in its absence. In the following consideration we use the normalized secondorder correlation function

$$g ^ { ( 2 ) } ( m, k ) = G ^ { ( 2 ) } ( m, k ) / G ^ { ( 1 ) } ( m ) G ^ { ( 1 ) } ( k ), \quad ( 5 ) \quad \text{where}$$

conditioned by the choice of dataset.

The ghost imaging implies knowledge of parameters (spatial resolution, correlation strength and background) that specify the correlation over SPAD pixels for the applied light beams. The correlation g (2) ( m ,k ∗ ) between a 'probing' low-noise pixel m = m ∗ , selected inside the "ghost" area (Fig. 1), and all other pixels of the SPAD array is well approximated by the Gaussian

$$g ^ { ( 2 ) } ( m ^ { * }, k ) \approx A _ { m ^ { * } } \exp \left \{ - \frac { | r _ { m ^ { * } k } | ^ { 2 } } { 2 B _ { m ^ { * } } ^ { 2 } } \right \} + C _ { m ^ { * } }, \quad ( 6 ) \quad \sin ^ { \infty c } \text{sim} \text{ }$$

with the distance | r m k ∗ | = √ ∆ x 2 m k ∗ +∆ y 2 m k ∗ between k -th and m ∗ -th pixels. Typical functions g (2) ( m ,k ∗ ) depending on the distance ∆ y m k ∗ acquired from different datasets for a '"probing" pixel m ∗ (the pixel is located in the middle of the SPAD array) are demonstrated in Fig. (3).

## B. Ghost image and its quality

In the work, we define the ghost image as a weighted sum of several correlation functions

$$\epsilon \Big | \quad \mathcal { I } _ { k } ^ { \prime } = \sum _ { m \in \{ B \} } G ^ { ( 1 ) } ( m ) g ^ { ( 2 ) } ( m, k ) / \sum _ { m \in \{ B \} } G ^ { ( 1 ) } ( m ). \quad ( 7 )$$

The pixels labeled by m are "bucket" pixels belonging to the set { B } . These pixels are chosen as located near the center of the object-beam spot on the SPAD area. The image pixels labeled by k are in the "ghost" area.

The scaling factor and background of the image I ′ k are dataset dependent. In order to minimize such an uncertainty and, hence, to facilitate analyzing the images produced from different datasets the array (7) is converted into a background-subtracted normalized image structurally identical to I ′ k :

$$\mathcal { I } ( x _ { k }, y _ { k } ) \equiv \mathcal { I } _ { k } = ( \mathcal { I } _ { k } ^ { \prime } - \bar { C } ) / \bar { A }. \text{ \quad \ \ } ( 8 )$$

In our experiment, quality assessment for the reconstructed ghost images is performed by analyzing the properties of I k . In accordance with (6), the spatial resolution of GI is limited by a quantity of ∼ 2 ¯ Bδ s .

We apply two measures specifying the image quality. The contrast-to-noise ratio ( CNR ) is a measure traditionally applied to characterize the quality of spatially restricted image parts. In imaging applications, these parts usually refer to transmitting ('In') or non-transmitting ('Out') areas of the imaged object. In our experiment, the pixels of these parts are addressed by the indices k belonging to the sets { In } or { Out } , respectively. The CNR defined as a ratio of the image contrast K to the total error σ is found in the form:

$$C N R = \frac { K } { \sigma } = \frac { \langle \mathcal { I } _ { k } \rangle _ { \{ I n \} } - \langle \mathcal { I } _ { k } \rangle _ { \{ O u t \} } } { \sqrt { \sigma _ { I n } ^ { 2 } + \sigma _ { O u t } ^ { 2 } } }, \quad \ ( 9 )$$

where the averaging ⟨ . . . ⟩ { In,Out } is carried out over the sets { In } or { Out } . The quantities σ 2 In,Out are corresponding variances of I k . At evaluating CNR , we cannot properly assess the influence of strong-noise pixels on the ghost image since they are sparsely dispersed over the SPAD array and the area 'In' is essentially small. Therefore, for definiteness, the sets { In } and { Out } includes only indices of low-noise pixels.

The quality of image can also be assessed with help of a correlation coefficient which pixel-by-pixel measures the strength of linear relationship (that is, a structural similarity [33]) between the analyzed image and a higherquality reference image of the same object. Unlike CNR ,

FIG. 4: A high-resolution conventional image (a) and a 32 × 32 -pixel reference image I ( r ) ( x , y k k ) (b) of the object.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01153v1-figure-95ecee1d-5c85-4591-be65-36f2c2f33ff6.jpg)

this measure is applicable to the entire image. We apply the correlation coefficient to evaluate the quality for the entire "ghost" area including highly noised pixels. To produce a reference image of suitable quality, we project a conventional image of the object (back-lighted by a laser beam) with help of a converging lens on the photosensitive matrix of a high-resolution digital camera (UI3240CP-NIR-GL Rev.2, DS Imaging Development Systems GmbH, pixel size of 5 3 . µm ) at a lens magnification of ∼ 1 5 . . A high-resolution image (spatial resolution of ∼ 3 5 . µm /pixel, this image is illustrated in Fig. 4(a)) stored by the digital camera is converted into a 32 × 32 -pixel image I ( r ) ( x , y k k ) ≡ I ( r ) k (resolution of ∼ 45 µm /pixel). The obtained high-contrast and lownoise image I ( r ) k (the image is shown in Fig. 4(b)) is then exploited as a reference in order to find the correlation coefficient:

$$R = \frac { \left | \sum _ { k \in \{ g h o s t \} } ( I _ { k } ^ { ( r ) } - \langle I _ { k } ^ { ( r ) } \rangle ) ( \mathcal { I } _ { k } - \langle \mathcal { I } _ { k } \rangle ) \right | } { \sqrt { \sum _ { l, k \in \{ g h o s t \} } ( I _ { l } ^ { ( r ) } - \langle I _ { l } ^ { ( r ) } \rangle ) ^ { 2 } ( \mathcal { I } _ { k } - \langle \mathcal { I } _ { k } \rangle ) ^ { 2 } } }, \quad ( 1 0 ) \quad \text{Th} _ { \substack { F _ { 0 }, \text{ is} } }.$$

where the summation is done over the "ghost" part of the SPAD area. In equation (10) brackets ⟨ . . . ⟩ denote averaging over the "ghost" area pixels. To achieve an accurate correspondence between the ghost and reference images, (i.e., to ensure the same position and the same scaling factor for the both images on the 32 × 32 pixel grid) the correlation coefficient R is maximized at image converting. The error of R evaluation is less than 0 01 . .

Overall the following procedure of ghost-image reconstruction is implemented. First of all, a dataset is prepared. The prepared dataset is applied to calculate G (1) ( m ) and G (2) ( m,k ) and find normalized correlation function g (2) ( m,k ) as described in equations (4) and (5). After that the averaged parameters ¯ A B C , ¯ , ¯ and A/C of correlation function (6) are evaluated. Then, a ghost image I k is computed in accordance with equations (7) and (8). Finally, quality of the obtained image is examined visually and estimated using CNR (9) and correlation

## IV. RESULTS AND DISCUSSION

In our GI experiment, the light beams are registered under conditions where ⟨ n f ⟩ ≈ 8 . The registered data stream F is collected over 186 minutes, stored on the computer and later post-processed. At post-processing, this dataset is used to prepare an unfiltered dataset F 0 and then to produce a number of frame-filtered sets F n ( ∗ ) .

The maximum count rate for the object-beam spot on the SPAD array is realized near the point (20 29) , (here and further we denote the pixels coordinates x and y just by a pair of indices ( x, y ) ). The "flare" area is taken as a rectangle with vertices at the pixels (12,32), (12,24), (26,24), (26,32). The rectangle contains 15 × 9 = 135 pixels. So, the remaining "ghost" area includes J max = 1024 -135 = 889 pixels. Twenty four low-noise pixels uniformly and randomly distributed across the "ghost" area are selected as the "probing" pixels m ∈ { m ∗ } required to determine the averaged parameters ¯ A B C , ¯ , ¯ and A/C .

In the study, we apply two sets { B } of "bucket" pixels. Eight low-noise pixels around the pixel (20 29) , and this pixel are used as a "bucket" set. On the other hand, to highlight advantages offered by our approach, we also demonstrate a traditional GI technique. For that we just choose as a "bucket" detector a single pixel in the "bucket" region. This pixel is the point (20 29) , . The "flare" area, the exploited "probing" and "bucket" pixels are indicated in Fig 5(a).

## A. GI with unfiltered data

The total number of frames in the unfiltered dataset F 0 is N 0 ∼ 1 9 . × 10 9 . The per-frame probabilities G (1) ( x , y k k ) ≡ G (1) ( k ) obtained for all pixels of the SPAD array from this dataset are shown in Fig. 5(a). The averaged values of shape parameters and the corresponding standard deviations are as follows: ¯ A ≈ 0 780 . , σ A ≈ 0 047 . ; ¯ B ≈ 1 72 . , σ B ≈ 0 085 . ; and ¯ C ≈ 1 085 . , σ C ≈ 0 0046 . . According to the obtained data, the spatial resolution of the image computed from the data F 0 is of ∼ 2 ¯ Bδ s ≈ 3 44 . δ s ≈ 0 15 . mm. The average perframe probability of counts detected by the "bucket" pixels ⟨ G (1) ( k ) ⟩ k ∈{ B } ≈ 0 057 . (this quantity is obtained by averaging of ⟨ G (1) ( k ) ⟩ over the nine 'bucket' pixels) is ∼ 300 times higher than the corresponding averaged probability of dark counts ⟨ G (1) DC ( k ) ⟩ k ∈{ B } ≈ 1 8 . × 10 -4 .

The computed ghost image I ( x , y k k ) is represented in Fig. 5(b). The visual examination of this figure shows that the imaged object is recognizable in the ghost image. However, the presented figure clearly demonstrates a serious drawback. The number of noisy pixels is so high

FIG. 5: GI based on the unfiltered dataset F 0 . The white x's show the pixels belonging to the "bucket' region. The white squares and black circles show, respectively, the pixels belonging to the sets { In } and { Out } . Panel (a) : distribution of log 2 [ G (1) ( x , y k k )] ; the white cross shows the point (20,29) corresponding to the maximal intensity of the object-beam spot; the white dashed line indicates a rectangle containing the 'flare' region; the white circles show the 'probing' pixels. Panel (b) : ghost image I ( x , y k k ) ; white crosses indicate positions of 9 pixels chosen as a "bucket". Panel (c) : ghost image reconstructed with the single-pixel "bucket" marked by x in the "flare" area.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01153v1-figure-755e4df8-1669-4820-9169-42a0827097d8.jpg)

that they have a strong destructive effect on the quality of the entire image. According to Appendix, the efficient number of pixels distinguishable (due to the dark-count effect) in the ghost image reconstructed from the dataset F 0 is of ∼ 100 . In addition, the spatial resolution is clearly suffering. The resolution is not high enough to resolve all parts of the object. The "tail" of the imaged digit "2" (a ∼ 0 1 . mm wide strip) is not visible in the reconstructed image.

The quality of the obtained image I k is estimated numerically in terms of CNR (9) and structure similarity with the reference image (10). In the CNR evaluation we exploit 6 low-noise pixels for the set { In } and 43 low-noise pixels for the set { Out } (these pixels are indicated in Fig. 5(b)). The { In } pixels are located near the maxima of the observed ghost image. For these pixels sets we have ⟨I k ⟩ { In } ≈ 0 1170 . , σ In ≈ 0 0046 . and ⟨I k ⟩ { Out } ≈ 0 0163 . , σ Out ≈ 0 0073 . . According to (9), the image CNR ≈ 11 7 . is achieved at the contrast of K ≈ 0 101 . and the total error σ ≈ 0 0086 . . Along our estimation (10), the coefficient of correlation R between the image I k and the reference image I ( r ) k is of ∼ 0 62 . .

The unfiltered data F 0 is a composition of frames with different values of n f . Fig. 6 shows the probability distribution that expresses a per-frame probability (3). The shape of this distribution is rather close to the shape of the Poisson distribution with the same average number of counts per frame. The obtained p n ( f = n ′ ) distribution has a maximum near n ′ = 7 and a half-maximum width of ∼ 9 .

## B. GI with filtered data

In this subsection, we exploit frame-filtered datasets F n ( ∗ ) that are produced by filtering the data F 0 consid-

FIG. 6: Distribution of per-frame probabilities p n ( f = n ′ ) over realizations of n f (blue squares) for the unfiltered dataset F 0 . The Poisson probability distribution of the same ⟨ n f ⟩ are shown by the red solid line. The Poisson distribution that represents the dark counts typical for the applied SPAD array is depicted by the black dashed line

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01153v1-figure-2dfc92d3-e743-4f23-8290-f8f8228ea984.jpg)

ered in the previous subsection. For the exploited data, the parameter n ∗ ranges from 2 to 24 (see Fig.6).

The datasets F n ( ∗ ) are used to clarify how the dataset composition may affect the correlation function (6). Fig.7 shows the dependence of the half-width ¯ B and the ratio A/C on n ∗ (the scatters show the corresponding rootmean-squared deviations). As can be seen from this figure, the parametric ratio A/C , within the error of its determination, does not depend on n ∗ . On the other hand, the half-width ¯ B monotonically increases with n ∗ from 1.4 to 2 pixels. As a result, the spatial resolution of GI will worsen from ∼ 0 12 . mm up to ∼ 0 18 . mm with the parameter n ∗ . The GI based on a frame-filtered dataset with n ∗ &lt; 12 must exhibit better spatial resolution than the GI reconstructed from the unfiltered data F 0 . The observed dispersion-like phenomenon for the half-width

FIG. 7: Dependence on n ∗ for averaged ratio A/C ( left Y axis: blue circles) and half-width ¯ B ( right Y axis: red stars) obtained from frame-filtered datasets F n ( ∗ ) . The relevant parameters acquired from the unfiltered data F 0 are shown by a blue dash-dotted line with a blue square (left Y axis: A/C ) and a red dashed line with a red square (right Y axis: ¯ B ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01153v1-figure-506db5be-98ec-40b0-89e6-e7ce1e2ef73c.jpg)

FIG. 8: Ghost images reconstructed from the frame-filtered datasets F (4) (a), F (8) (b) and F (16) (c). Ghost image computed from the frame-filtered dataset F (8) with the single-pixel "bucket" (d). The white squares, circles and x's show, respectively, the pixels belonging to the sets { In } , { Out } and { B } .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01153v1-figure-555f5610-5b73-4703-869a-01504cd5a5df.jpg)

¯ B may be associated with specific optical properties of the rotating disk D. The datasets F n ( ∗ ) are used to reconstruct ghost images. Fig. 8 demonstrates images I k computed from frame-filtered datasets F n ( ∗ ) at n ∗ = 4 , 8 and 16 . In terms of contrast, spatial resolution, and number of distinguishable noised pixels, the presented images look significantly better compared to the ghost image obtained from unfiltered data. Image qualities (such as spa-

FIG. 9: CNR (left Y axis: blue stars) and correlation coefficient R (right Y axis: red squares) as functions of n ∗ for ghost images reconstructed from frame-filtered datasets F n ( ∗ ) . The blue dashed and red solid lines indicates CNR (left Y axis) and R (right Y axis) for the unfiltered data F 0 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01153v1-figure-998a69f6-2402-4fc7-b0ec-55c6cf967dc4.jpg)

tial resolution and contrast) depend on the parameter n ∗ . Based on visual inspection of the presented images one can conclude that ghost images of the highest quality are realized for datasets F n ( ∗ ) with n ∗ ∼ 8 . As can be seen from Fig 8 (b), the spatial resolution of the image I k based on the dataset F(8) is high enough to resolve all parts of the imaged object. In accordance with Appendix, only a few pixels of the image reconstructed from F(8) can be identified as pixels affected by dark counts.

This conclusion is corroborated with the results of numeric evaluation of CNR (9) and the correlation coefficient (10) represented in Fig. 9. In the CNR evaluation, the pixels k ∈ { In } and k ∈ { Out } are the same pixels as in the previous subsection. Generally, for the ghost image based on the frame-filtered data, the maximum CNR and correlation coefficient R are achieved near those values of n ∗ that correspond to the maximum probability p n ( f = n ′ ) of the distribution analyzed in the previous subsection and represented in Fig.6. As shown in Fig. 9 (left Y axis), compared to the "unfiltered" image, the "frame-filtered" image I k has a higher CNR when the parameter n ∗ ranges from 5 to 10. At n ∗ = 8 , the filtered GI reaches its maximal CNR ≈ 17 0 . ( ⟨I k ⟩ { In } ≈ 0 0744 . , σ In ≈ 0 0033 . , ⟨I k ⟩ { Out } ≈ 0 0009 . , σ Out ≈ 0 0028 . , K ≈ 0 0735 . , σ ≈ 0 0043 . ) which is nearly 45% higher as compared to the CNR of unfiltered GI.

Unlike CNR , the measure of structural similarity (10) takes into account the effect of strong-noise pixels on the image. Therefore, the range 3 ≤ n ∗ ≤ 19 , where the "low-noise" filtered GI leads to a higher coefficient R (as compared to the "noisy" unfiltered GI), covers almost all of the considered values of n ∗ . Like CNR , the structural similarity between the reference and ghost images attains its maximum at n ∗ = 8 leading to R ≈ 0 8 . (right Y axis of Fig. 9).

## C. GI with a single-pixel "bucket"

The per-frame probability of counts detected by the single "bucket" pixel G (1) (20 29) , ≈ 0 069 . obtained from the unfiltered data F 0 is substantially high compared to the corresponding probability of dark counts G (1) DC (20 29) , ≈ 4 5 . × 10 -5 . The ghost images computed with the single-pixel "bucket" from the unfiltered data F 0 and the filtered dataset F (8) are shown in Fig 5(c) and Fig 8(d), respectively. The resulting ghost image obtained with the data F 0 is the image that one gets with traditional ghost imaging technique for the case. It seems close to the GI of Fig 5(b) obtained with our approach for the data F 0 . Indeed, the correlation coefficient is only slightly lower, R ≈ 0 6 . in comparison to the correlation for the multiple-pixel "bucket" image ( R ≈ 0 62 . ). However, CNR is nearly twice lower, CNR ≈ 6 73 . .

Using the filtered dataset F (8) improves the correlation giving R ≈ 0 724 . . However, simultaneously it leads to degradation of the contrast, CNR ≈ 5 46 . . Thus, a single-pixel "bucket" imaging modality for the case cannot ensure the quality of image given by averaging over several GI in our approach.

## V. CONCLUSIONS

We have presented an approach to ghost imaging and analyzed the approach capabilities in an experiment where the light intensity is so low that the dark counts have a noticeable effect on the imaging. We have demonstrated that filtering the acquired data with respect to the number of counts per frame and obtaining the image as a weighted average over the images for different "bucket" pixels, one can significantly reduce the destructive effect of dark counts on the ghost image, improve the image contrast, spatial resolution and image similarity to the reference image. The proposed approach allowed us inferring ghost images for dark counts reaching 20% of the the total "bucket" counts and with the photon fluxes of just few photons per time-frame. In our opinion, using multiple bucket detectors instead of just one is compatible with most applications of ghost imaging and will allow one to greatly enhance imaging quality in the extremely-low-light condition. These results can be extended to computational [5, 11-13] and temporal [34-36] ghost imaging. Our findings are important for developing correlation imaging schemes with light on a-few-photon level, first of all, for quantum sensing applications.

We feel that as a concluding remark it would be useful to emphasize our motivation. As it is shown in Appendix, the proposed frame-filtering approach improves considerably the reconstructed ghost images and, simultaneously, has no effect on the quality of G (1) . Thus, we do not expect any improvement compared to the conventional imaging with the same light source and detection system [37]. Our motivation was to improve the quality of the ghost imaging with very low photon fluxes using very realistic noisy SPADs with inhomogeneous quality of pixels. We see this task as quite actual and practical for improving ghost imaging in its application niche.

## FUNDING

V.S., S.K., V.C. and A. S. acknowledge financial support from the State Research Program "Convergence" subprogram 3.01.2. D.M. acknowledges financial support from the BRFFR project F23UZB-064. This work was supported by the European Commission through the SUPERTWIN Project under Grant 686731.

## 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.

## Appendix A: Dark counts of the SPAD array and their effect on reconstructed ghost images

To analyze the effect of dark-counts (DC) for the SPAD array we, first of all, measure the probability of counts for each SPAD pixel under conditions where the detection events are DC. To ensure this condition, the SPAD array is blocked by a light-tight screen from any light source (including ambient light). The registered data are stored in the computer as a frame-unfiltered dataset F DC including ∼ 2 8 . × 10 9 frames. The per-frame probabilities of DC G (1) DC ( k ) for all SPAD pixels are found from the dataset F DC , as described in Section III. The distribution of measured G (1) DC ( x , y k k ) over the pixels of the SPAD array is shown in Fig 10. The bright pixels of this distribution are strongly noised pixels of the SPAD array. The brightest pixel is the pixel k = 212 with coordinates x k = 7 , y k = 20 . One can see in Fig 10 that distribution of DC rates over the SPAD pixels is quite inhomogeneous and randomized. The measured probabilities G (1) DC ( k ) are in a wide range of values from ∼ 10 -5 to ∼ 10 -2 . The average DC probability per pixel and per frame reaches a value of ⟨ G (1) DC ⟩ ≈ 4 9 . × 10 -4 (or approximately 0.5 per frame). The dependence of per-frame probability p n ( f = n ′ ) on n ′ for this data is represented by the Poisson distribution with ⟨ n f ⟩ ≈ 0 5 . and shown in Fig.6.

In general, despite the substantial inhomogeneity and randomness, the influence of DC on the measurement can be evaluated if the evaluation operates with specially

FIG. 10: Per-frame dark counts distribution G (1) over the SPAD pixels. The colorbar shows the corresponding normalized values.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01153v1-figure-0d151083-4992-4c53-85f9-5a5416dbe96a.jpg)

FIG. 11: 'DC-ordered' sequences G (1) DC ( S j ) (a, black dashed line), G (1) ( S j ) (a) and ∆ ( I S j ) (b) as functions of the serial number j . The red circles show the sequences calculated from the unfiltered data F 0 . The blue solid line and dots indicate the sequences found from the frame-filtered dataset F (8) . The inset of panel (a) demonstrates the same data (as in the main panel) at a larger Y-axis scale. The red dashed line in this inset indicates the average probability ⟨ G (1) DC ⟩ . The inset of panel (b) shows the same data (as in the main panel) in a log scale on the X-axis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01153v1-figure-b976b151-2fd5-408c-a452-302fa83623ea.jpg)

organized sequences of pixels. In these sequences, the pixels are arranged so that their dark count rates G (1) DC ( k ) are ordered in descending order ('DC-ordering'):

$$G _ { D C } ^ { ( 1 ) } ( k ) \xrightarrow { D C - o r d e r i n g } G _ { D C } ^ { ( 1 ) } ( S _ { j } ), \quad ( \text{A} 1 ) \quad \overset { \psi \nu } { \text{is} \nu_{\ell} }$$

where j is the serial number for each term of the se- quence G (1) DC ( S j ) . The sort index S j specifies how the data G (1) DC ( k ) should be rearranged to obtain the ordered sequence G (1) DC ( S j ) . The first terms of the DC-ordered sequence belong to strongly noisy pixels, the last to lownoise pixels.

The DC effect on a measured quantity Q k ( ) is evaluated by analyzing a 'DC-ordered' sequence Q S ( j ) . For the pixels with a weak influence of DC on Q (e.g., for low-noise pixels with j ≫ 1 ) the measured Q S ( j ) fluctuates around a certain value Q 0 , independent of j . On the other hand, for a pixel k ′ where the influence of DC on Q is significant, the quantity Q k ( ′ ) is a function of G (1) DC ( k ′ ) . For several pixels with strong DC effect, the DC-ordered sequence Q S ( j ) explicitly reproduces the sequence G (1) DC ( S j ) within a j index range ν Q . In a sense, the deviation χ Q = max ( | Q S ( j ) -Q 0 | ) realized near j = 1 and the width of the range ν Q can be accepted as measures of the DC effect. The width ν Q can be associated, for instance, with an efficient number of noised pixels that have a destructive effect on the measured quantity Q . The deviation χ Q characterizes the effect strength for the "hot" pixels.

To evaluate the DC effect on a reconstructed ghost image I ( x , y k k ) we compute a pixel-by-pixel image difference ∆ ( ) I k ≡ ∆ ( I x , y k k ) = ˜ ( I x , y k k ) - I ( x , y k k ) . Here, ˜ ( I x , y k k ) is an image obtained from I ( x , y k k ) by a 2D smoothing in which the convolution kernel is a 3 × 3 pixel square. The quantity ∆ ( ) I k is insensitive to pattern variation for the imaged object and approximately equal to fluctuation in I ( x , y k k ) between the neighboring pixels. Owing to these properties, the difference ∆ I is well adapted to register one-pixel-ranged correlation drops typical for strongly noised pixels sparsely dispersed across the SPAD array. We analyze the dependence of the 'DC-ordered' difference ∆ ( I S j ) on the index j . In a similar way, the DC influence on the measured probabilities G (1) can be evaluated by considering the sequence G (1) ( S j ) as a function of j . The 'DCordered' sequences G (1) ( S j ) and ∆ ( I S j ) calculated for the pixels k ∈ { ghost } ( j = 1 , . . . , J max ) from framefiltered and unfiltered datasets (they are datasets F (8) and F 0 applied for ghost image computation in Section IV) are represented in Fig. 11. For illustrative purposes (to demonstrate the effect manifestation in comparison with the image contrast K ), the difference ∆ ( I S j ) in the figure is normalized by K .

Regardless of the type of filtering (with or without frame filtering), the DC have an equal influence on the probabilities G (1) obtained from data with the same average numbers ⟨ n f ⟩ . As shown in Fig. 11 (a), the dependencies of G (1) ( S j ) on the j index for filtered and unfiltered data ( F (8) and F 0 with ⟨ n f ⟩ ≈ 8 ) are very similar. The efficient number of pixels distinguishable due to the DC effect in the imaged distribution G (1) ( x , y k k ) is ν G 1 ∼ 100 for the both datasets. Contrary to G (1) , the manifestation of DC effect in ghost images (see Fig. 11 (b)) depends significantly on what data (frame-filtered or

not) is used to compute the images. The efficient number of pixels distinguishable (due to the DC effect) in the ghost image computed from the unfiltered data F 0 is ν ∆ I ∼ ν G 1 ∼ 100 . For the frame-filtered dataset F (8) , the manifestation of the DC effect is essentially small. In fact, only a few pixels can be identified that are affected

- [1] T. B. Pittman, Y. H. Shih, D. V. Strekalov, and A. V. Sergienko, Optical imaging by means of two-photon quantum entanglement, Phys. Rev. A 52 , R3429 (1995).
- [2] R. S. Bennink, S. J. Bentley, and R. W. Boyd, Twophoton coincidence imaging with a classical source, Phys. Rev. Lett. 89 , 113601 (2002).
- [3] A. Valencia, G. Scarcelli, M. D'Angelo, and Y. Shih, Twophoton imaging with thermal light, Phys. Rev. Lett. 94 , 063601 (2005).
- [4] A. Gatti, M. Bache, D. Magatti, E. Brambilla, F. Ferri, and L. A. Lugiato, Coherent imaging with pseudothermal incoherent light, Journal of Modern Optics 53 , 739 (2006).
- [5] B. I. Erkmen and J. H. Shapiro, Ghost imaging: from quantum to classical to computational, Adv. Opt. Photon. 2 , 405 (2010).
- [6] J. H. Shapiro and R. W. Boyd, The physics of ghost imaging, Quantum Information Processing 11 , 949 (2012).
- [7] M. J. Padgett and R. W. Boyd, An introduction to ghost imaging: quantum and classical, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences 375 , 20160233 (2017).
- [8] D. S. Simon, G. Jaeger, and A. V. Sergienko, Ghost imaging and related topics, in Quantum Metrology, Imaging, and Communication (Springer International Publishing, Cham, 2017) pp. 131-158.
- [9] G. M. Gibson, S. D. Johnson, and M. J. Padgett, Singlepixel imaging 12 years on: a review, Opt. Express 28 , 28190 (2020).
- [10] M. Padgett, Quantum imaging overview, in Quantum Sensing, Imaging, and Precision Metrology , Vol. 12447, edited by J. Scheuer and S. M. Shahriar, International Society for Optics and Photonics (SPIE, 2023) p. 1244702.
- [11] J. H. Shapiro, Computational ghost imaging, Phys. Rev. A 78 , 061802 (2008).
- [12] M. P. Edgar, G. M. Gibson, and M. J. Padgett, Principles and prospects for single-pixel imaging, Nature Photonics 13 , 13 (2019).
- [13] C. Bruschini, H. Homulle, I. M. Antolovic, S. Burri, and E. Charbon, Single-photon avalanche diode imagers in biophotonics: review and outlook, Light: Science and Applications 8 , 1 (2019).
- [14] M. Zarghami, L. Gasparini, L. Parmesan, M. MorenoGarcia, A. Stefanov, B. Bessire, M. Unternährer, and M. Perenzoni, A 32 x 32-pixel CMOS imager for quantum optics with Per-SPAD TDC, 19.48 % fill-factor in a 44.64µ m pitch reaching 1-MHz observation rate, IEEE Journal of Solid-State Circuits 55 , 2819 (2020).
- [15] K. Morimoto, J. Iwata, M. Shinohara, H. Sekine, A. Abdelghafar, H. Tsuchiya, Y. Kuroda, K. Tojima, W. Endo, Y. Maehashi, Y. Ota, T. Sasago, S. Maekawa, S. Hikosaka, T. Kanou, A. Kato, T. Tezuka, S. Yoshizaki, T. Ogawa, K. Uehira, A. Ehara, F. Inui, Y. Matsuno,

by DC. The deviation χ ∆ I realized for the "hottest" pixel k = 212 is more than three times smaller than for unfiltered data. Such a dramatic reducing (compared to F 0 ) in the DC effect is observed for all frame-filtered datasets used in Section IV at GI reconstructing. All pixels with the index j &gt; 300 are accepted to be low-noise.

- K. Sakurai, and T. Ichikawa, 3.2 megapixel 3D-stacked charge focusing SPAD for low-light imaging and depth sensing, 2021 IEEE International Electron Devices Meeting (IEDM) , 20.2.1 (2021).
- [16] J. T. Smith, A. Rudkouskaya, S. Gao, J. M. Gupta, A. Ulku, C. Bruschini, E. Charbon, S. Weiss, M. Barroso, X. Intes, and X. Michalet, In vitro and in vivo NIR fluorescence lifetime imaging with a time-gated SPAD camera, Optica 9 , 532 (2022).
- [17] H. Defienne, P. Cameron, B. Ndagano, A. Lyons, M. Reichert, J. Zhao, A. Harvey, E. Charbon, J. Fleischer, and D. Faccio, Pixel super-resolution with spatially entangled photons, Nature Communications 13 (2022).
- [18] M. Unternährer, B. Bessire, L. Gasparini, M. Perenzoni, and A. Stefanov, Super-resolution quantum imaging at the heisenberg limit, Optica 5 , 1150 (2018).
- [19] B. Eckmann, B. Bessire, M. Unternährer, L. Gasparini, M. Perenzoni, and A. Stefanov, Characterization of space-momentum entangled photons with a time resolving CMOS SPAD array, Opt. Express 28 , 31553 (2020).
- [20] C. Abbattista, L. Amoruso, S. Burri, E. Charbon, F. Di Lena, A. Garuccio, D. Giannella, Z. Hradil, M. Iacobellis, G. Massaro, P. Mos, L. Motka, M. Paúr, F. V. Pepe, M. Peterek, I. Petrelli, J. Řeháček, F. Santoro, F. Scattarella, A. Ulku, S. Vasiukov, M. Wayne, C. Bruschini, M. D'Angelo, M. Ieronymaki, and B. Stoklasa, Towards quantum 3D imaging devices, Applied Sciences 11 , 10.3390/app11146414 (2021).
- [21] B. Ndagano, H. Defienne, D. Branford, Y. Shah, A. Lyons, N. Westerberg, E. Gauger, and D. Faccio, Quantum microscopy based on Hong-Ou-Mandel interference, Nature Photonics 16 , 384 (2022).
- [22] J. Zhao, A. Lyons, A. C. Ulku, H. Defienne, D. Faccio, and E. Charbon, Light detection and ranging with entangled photons, Opt. Express 30 , 3675 (2022).
- [23] L. H. C. Braga, L. Gasparini, L. Grant, R. K. Henderson, N. Massari, M. Perenzoni, D. Stoppa, and R. Walker, A fully digital 8 × 16 SiPM array for PET applications with per-pixel TDCs and real-time energy output, IEEE Journal of Solid-State Circuits 49 , 301 (2014).
- [24] A. Carimatto, A. Ulku, S. Lindner, E. Gros-Daillon, B. Rae, S. Pellegrini, and E. Charbon, Multipurpose, fully integrated 128 × 128 event-driven MD-SiPM with 512 16-bit TDCs with 45-ps LSB and 20-ns gating in 40nm CMOS technology, IEEE Solid-State Circuits Letters 1 , 241 (2018).
- [25] H. Xu, L. Pancheri, G.-F. D. Betta, and D. Stoppa, Design and characterization of a p + /n-well SPAD array in 150nm CMOS process, Opt. Express 25 , 12765 (2017).
- [26] S. Magnitskiy, D. Agapov, and A. Chirkin, Quantum ghost polarimetry with entangled photons, Opt. Lett. 47 , 754 (2022).
- [27] S. Magnitskiy, D. Agapov, and A. Chirkin, Ghost po-

larimetry with unpolarized pseudo-thermal light, Opt. Lett. 45 , 3641 (2020).

- [28] S. A. Magnitskiy, D. P. Agapov, I. A. Belovolov, P. P. Gostev, D. N. Frolovtsev, and A. S. Chirkin, Ghost Polarimetry in Classical and Quantum Light, Moscow University Physics Bulletin 76 , 424 (2021).
- [29] https://www.savazzi.net/photography/1951usaf. html .
- [30] C. Pitsch, D. Walter, L. Gasparini, H. Bürsing, and M. Eichhorn, 3D quantum ghost imaging, Appl. Opt. 62 , 6275 (2023).
- [31] V. F. Gili, D. Dupish, A. Vega, M. Gandola, E. Manuzzato, M. Perenzoni, L. Gasparini, T. Pertsch, and F. Setzpfandt, Quantum ghost imaging based on a 'looking back' 2D SPAD array, Appl. Opt. 62 , 3093 (2023).
- [32] V. S. Starovoitov, V. N. Chizhevsky, D. B. Horoshko, and S. Y. Kilin, Spatio-Temporal Correlations of Pho-
- tons from a Pseudo-Thermal Source, Journal of Applied Spectroscopy 90 , 377 (2023).
- [33] Z. Wang, A. C. Bovik, H. R. Sheikh, and E. P. Simoncelli, Image quality assessment: from error visibility to structural similarity, IEEE transactions on image processing 13 , 600 (2004).
- [34] P. Ryczkowski, M. Barbier, A. T. Friberg, J. M. Dudley, and G. Genty, Ghost imaging in the time domain, Nature Photonics 10 , 167 (2016).
- [35] F. Devaux, P.-A. Moreau, S. Denis, and E. Lantz, Computational temporal ghost imaging, Optica 3 , 698 (2016).
- [36] D. B. Horoshko, Time-to-space ghost imaging, Opt. Lett. 48 , 3247 (2023).
- [37] P.-A. Moreau, E. Toninelli, P. A. Morris, R. S. Aspden, T. Gregory, G. Spalding, R. W. Boyd, and M. J. Padgett, Resolution limits of quantum ghost imaging, Opt. Express 26 , 7528 (2018).
V. S. Starovoitov V. N. Chizhevsky D. Mogilevtsev A. Smaliakou M. Perenzoni L. Gasparini D. B. Horoshko S. Kilin
Frame-filtered ghost imaging with a SPAD array used both as a multiple "bucket" detector and an imaging camera
An approach to ghost imaging with a single SPAD array used simultaneously as
a several-pixel "bucket" detector and an imaging camera is described. The key
points of the approach are filtering data frames used for ghost-image
reconstruction by the number of per-frame counts and superposing correlation
images obtained for different "bucket" pixels. The imaging is performed in an
experiment with a pseudo-thermal light source where the light intensity is so
low that the dark counts have a noticeable effect on imaging. We demonstrate
that the approach is capable to significantly reduce the destructive effect of
dark counts on the ghost image and improve image contrast, spatial resolution,
and image similarity to a reference image.