## Abstract

Single-pixel imaging (SPI) is a novel method for capturing high-quality 2D images of scenes using a non-spatially-resolved detector. While implementing conventional SPI, a huge number of illuminated patterns are projected onto the object to reconstruct a sharp image. For a situation in which the object occupies part of the illuminated region, we propose an adaptive regional SPI method (ARSI) to decrease the number of projected patterns. In the ARSI scheme, the object region is adaptively located based on the Fourier slice theorem. Then, the illuminated patterns are projected only onto the object region to facilitate imaging efficiency. Experiments demonstrate that the proposed ARSI method can achieve sharp image reconstruction with a substantial reduction in pattern number, thereby improving imaging efficiency.

© 2017 Optical Society of America

## 1. Introduction

Single-pixel imaging (SPI) is a novel method for capturing 2D images using a single-pixel detector without spatial resolution. In SPI, known illuminated patterns are projected onto an object, and the intensity of the reflected light is measured using a single-pixel detector [1]. Then, the image can be reconstructed by calculating the correlation between projected patterns and recorded signals [2]. SPI has been applied to many fields, including 3D imaging [3–5], terahertz imaging [6,7], and multispectral imaging [8], due to its excellent anti-noise and high signal-to-noise ratio (SNR) characteristics.

Despite these excellent characteristics, SPI is required to project large amounts of patterns to reconstruct a sharp image. We divide SPI into two categories based on the orthogonality of illuminated patterns: SPI using patterns generated from non-orthogonal bases and SPI using patterns generated from orthogonal bases. The first SPI method typically uses random patterns [9], which actually comprise an overcomplete set. In this method, a huge number of patterns should be projected to reconstruct a sharp image because of the non-orthogonality of the patterns [1]. Sun [3] used 1 million random patterns to reconstruct an image with a resolution of 256 × 192. Hence, the numbers of used patterns are considerably higher than the resolution of the reconstructed image.

SPI that uses patterns from complete orthogonal bases includes SPI based on Fourier transform [10], discrete cosine transform [11], and Hadamard transform [12]. The number of patterns from a complete orthogonal basis is less than that from an overcomplete non-orthogonal basis. Thus, the use of patterns generated from a complete orthogonal basis can reduce the projection counts of SPI compared with the former SPI method. Recently, the Fourier SPI (FSI) method [10], which uses Fourier basis patterns for illumination, is proposed to achieve high-quality 2D imaging. However, data are obtained by sequence for all existing SPI methods, whereas millions of data are acquired simultaneously using conventional digital cameras. The time consumption of SPI is proportional to the total number of patterns given a certain rate of pattern projection, thereby more patterns require more time for projection.

A large number of projected patterns create a bottleneck in the application of SPI. The compressive sampling [13] method can be used to decrease the number of projected patterns, whereas the signal-to-noise ratio (SNR) of the reconstructed image will be reduced. Bian *et al.* [14] developed an efficient SPI technique (eSPI) to reduce the number of patterns. This technique captures only the most informative data in spatial frequency bands.

The number of projected patterns is reduced if the resolution of the image is decreased because the number of projected patterns for SPI is proportional to the spatial resolution of the reconstructed image. In fact, the object to be imaged typically takes up a partial region in the reconstructed image [1,3]. If SPI is conducted only on the object region, then the number of projected patterns and the time consumption will be significantly decreased without degrading image quality. This finding motivates us to develop a method for locating the object region in a scene and adaptively conduct SPI on the object region to reduce the number of projected patterns.

In this study, we propose an adaptive regional single-pixel imaging (ARSI) method to facilitate the efficiency of SPI. The proposed method applies the Fourier slice theorem [15] to adaptively locate the object region and then conducts SPI only on the object region. Compared with the conventional SPI method that images the whole scene, the proposed method significantly reduces the number of patterns, and thus, increases imaging efficiency.

The rest of this paper is organized as follows. The principles of the proposed method are explained in Section 2. The measurement system and the experimental results are discussed in Section 3. Lastly, the conclusions drawn from the study are presented in Section 4.

## 2. Principles

The proposed ARSI method is suitable for the situation in which an object occupies a partial area in the scene region. Unlike the conventional SPI method that images the whole scene [Fig. 1(a)], ARSI adaptively conducts SPI on the object region after locating the object region [Fig. 1(b)]. ARSI can reduce the number of illuminated patterns and time consumption via adaptive regional imaging.

The framework for ARSI comprises the following steps.

- ● Step 2. Conducting SPI only on the located object region. After the object region is located, the image of the object can be reconstructed by conducting any SPI method on the object region. In this study, FSI is used to obtain the object image.
- ● Step 3. Placing the reconstructed image on the object location in the scene. The reconstructed image is placed on the object location in the scene using the calculated object region location to generate a complete image.

#### 2.1 Adaptively locating the object region in the scene based on the Fourier slice theorem

In accordance with the Fourier slice theorem, the 1D Fourier transform of the projection is a slice of the 2D Fourier transform of the image along the same orientation as the projection [15]. As shown in Fig. 2, the projection of the scene image *f* (*x*,*y*) [Fig. 2(a)] is obtained via the 1D inverse Fourier transform of the values extracted from *F*(*u*,*v*) [Fig. 2(d)] along the line oriented at the same angle used to generate the projection. The projection of the image *f* (*x*,*y*) onto the *x*-axis is expressed as *p*(*x*) [Fig. 2(c)], whereas the projection onto the *y*-axis is *p*(*y*) [Fig. 2(b)]. Projections *p*(*x*) and *p*(*y*) can be obtained using

*IFT*represents the 1D inverse Fourier transform,

*F*(

*u*,

*v*) denotes the 2D Fourier coefficients of the image, and

*F*(

*u*,0) and

*F*(0,

*v*) are the slices of the Fourier coefficients. After

*F*(

*u*,0) and

*F*(0,

*v*) are acquired, projections

*p*(

*x*) and

*p*(

*y*) can be calculated based on the Fourier slice theorem. The method for acquiring

*F*(

*u*,0) and

*F*(0,

*v*) is introduced in the next part.

On the projection line of the scene, the mutation phenomenon occurs on the positions where the object edges are located because the grayscale distribution differs between the object and the background. The edges of the object region are located by using the edge detection algorithm [16] on projections *p*(*x*) and *p*(*y*); the edges are shown as points surrounded by a red circle in Figs. 2(b) and 2(c). The object region is located by the red rectangle shown in Fig. 2(a). To guarantee a correct object location, the object should have a high contrast with the background which has a uniform grayscale distribution.

Let the object region be Ω_{a}. Then, we have

*x*and

_{1}*x*represent the edges of the object region along the

_{2}*x*-axis, and

*y*and

_{1}*y*represent the edges of the object region along the

_{2}*y*-axis.

#### 2.2 Acquisition of the slice of the Fourier coefficients

We adopt the FSI method to acquire the required Fourier coefficients. As shown in Fig. 3, vertical and horizontal 2D sinusoidal patterns, i.e., *PV _{φ} and PH_{φ}*, are projected onto the scene to obtain the slice of the Fourier coefficients

*F*(

*u,*0) and

*F*(0

*,v*). The patterns of

*PV*and

_{φ}*PH*are generated using Eqs. (4) and (5) as follows:

_{φ}*x*,

*y*) represents the 2D Cartesian coordinates in the scene, (

*u*,

*v*) represents the spatial frequency,

*a*represents the DC component,

*b*denotes the contrast of the pattern, and

*φ*denotes the initial phase.

The response of the single-pixel detector to the illuminated sinusoidal patterns can be expressed by Eqs. (6) and (7) as follows:

*RV*and

_{φ}*RH*respectively represent the response of the single-pixel detector to the reflected light from the vertical and horizontal 2D sinusoidal patterns,

_{φ}*k*is a coefficient related to the sensitivity and location of the single-pixel detector,

*R*is the response of the single-pixel detector to environmental light, and Ω represents the scene region.

_{n}The light reflected from four patterns, which has the same frequency and π/2 phase shifting, is collected individually by the single-pixel detector. In accordance with the FSI approach, the slice of the Fourier coefficients can be calculated using Eqs. (8) and (9) as follows:

Moreover, *p*(*x*) and *p*(*x*) are calculated using Eqs. (1) and (2). Then, the object region Ω_{a} can be located using the edge detection algorithm on *p*(*x*) and *p*(*y*).

#### 2.3 Image reconstruction

To reconstruct the object image, we conduct SPI on the located object region. Any type of SPI method is suitable for reconstructing images. Similar to the conventional SPI method, the illuminated patterns with the same size as the object region are projected onto the object region. Once the single-pixel detector accomplishes data acquisition, the object image can be reconstructed using the corresponding reconstruction method of SPI. In this study, we conduct FSI on the object region, and the object image is reconstructed by performing 2D inverse Fourier transform on the acquired coefficients. Let the reconstructed object image be *f _{o}*(

*x,y*).

To acquire a complete image, a blank image with the same resolution as the scene image is generated as follows:

Then, the complete image can be obtained by*M*,

*N*) is the resolution of the object region, (

*x*,

_{1}*y*) is the location of the object region (Fig. 4).

_{1}#### 2.4 Number of projected patterns

When conducting SPI on the whole illuminated region, the count of the SPI projected patterns, *C*, is proportional to the resolution of the reconstructed image and can be expressed as

*k*and

_{1}*k*are the scale factors that are equal to or greater than one. The total number of projected patterns of ARSI,

_{2}*C*, iswhere

_{a}*C*represents the number of illuminated patterns while locating the object region,

_{L}*C*is the number of illuminated patterns while conducting SPI on the object region.

_{S}*C*and

_{L}*C*can be calculated using Eqs. (15) and (16) as follows:

_{S}In most cases, *C _{L}* is sufficiently small compared with

*C*, therefore the time consumption of ARSI is mainly determined by the number of illuminated patterns

_{S}*C*. Through the adaptive projection, the number of illuminated patterns

_{S}*C*can be greatly reduced and the efficiency of SPI can be significantly improved.

_{S}## 3. Experiments

The structure of our experiment setup is shown in Fig. 5. To produce high-contrast structured patterns, a digital projector is selected as the SLM. The used digital projector contains a light source with a wavelength of 455 nm. An illuminated pattern is projected onto the scene by the digital projector every 0.1 s. The reflected light is collected by a lens, and detected by the single-pixel detector, which is actually a photodiode. The 455 nm bandpass filter is fixed in front of the single-pixel detector. The collected analog signal is converted into a digital signal by an analog-to-digital converter (ADC) and processed using a computer. The object is placed 0.5 m away from the lens of the digital projector.

We choose the FSI approach as the SPI method to reconstruct the object image given its anti-noise and high-quality imaging characteristics. To acquire a coefficient, four patterns with a π/2 phase that shifts in the same frequency are projected onto the scene by the digital projector and the reflected light is detected sequentially by the single-pixel detector. The use of a four-step phase-shifting technology [10] can eliminate the effect of environment illumination and enhance measurement accuracy.

#### 3.1 Experiment on object location via ARSI

To verify the environmental adaptability of the proposed location method, the object location experiments with dark background, colored background and non-background were conducted respectively.

Firstly, the white toy bear with dark background, as shown in Fig. 6(a), was located. The resolution of scene image is 320 × 200. After projecting vertical patterns and horizontal patterns, the projection of the scene image in both directions was obtained; the result is shown in Figs. 6(b) and 6(c). For this certain resolution of the reconstructed image, the number of projected patterns for object location is invariable. In the projection line [blue lines shown in Figs. 6(b) and 6(c)], the positions where the mutation phenomenon occurs are the edges of the object region. The edges of the object are located by the red circle in Figs. 6(b) and 6(c) by using the edge detection algorithm. The object region is obtained as follows:

The located object region, which is illuminated by the digital projector, is shown in Fig. 6(d). The illuminated region is larger than the real object region to guarantee that the whole object region can be reconstructed.

The second experiment result with colored background is shown in Fig. 7. The resolution of scene image is 320 × 200. By applying proposed location method, object location is obtained as Eq. (18). In this experiment, the background of the object has a uniform grayscale distribution. As shown in Figs. 7(b) and 7(c), the projection lines on the non-object region have a gentle fluctuation. Whereas the object still can be located correctly by the proposed method.

The third experimental result with non-background is shown in Fig. 8, and the located object region is represented by Eq. (19). As Fig. 8(a) shows, only the light reflected from object can be collected by the detector. The projection lines [Figs. 8(b) and 8(c)] on the non-object region is near 0, since few light was reflected from non-object region. And the ambient light was greatly eliminated by bandpass filter. After location, a white board was posed behind the object to better reveal the located object region.

The proposed ARSI method has an excellent environmental adaptivity, which has been proved by the above experiments. The proposed method can be widely used in the situation where the object has a background with uniform grayscale distribution or non-background, since the object edges are judged by the mutation on the projection lines of scene image.

#### 3.2 Experiment on image reconstruction

After locating the object region, FSI was conducted on the located object region to acquire the coefficients of the image in the Fourier domain. The object image was obtained by conducting inverse 2D Fourier transform on the coefficients.

We first acquired the whole coefficients in Fourier domain to reconstruct the image. Figure 9(a) shows the final reconstructed image with a resolution of 320 × 200 via ARSI, whereas the resolution of the located object region is 112 × 90. During ARSI, 20160 patterns were projected to obtain the whole Fourier coefficients and reconstruct the object image. For comparison, we conducted normal FSI on the whole scene. The reconstructed result is shown in Fig. 9(b). A total of 128000 patterns were used to obtain the whole coefficients in the Fourier domain.

As shown in Figs. 9(a) and 9(b), the object image obtained via ARSI is nearly the same as that obtained via conventional FSI, whereas the projected number of patterns for normal FSI is 6.3 times higher than that for ARSI. The time consumption of SPI is proportional to the projected number of patterns; hence, ARSI used 1/6.3 of the projection and the acquisition time of normal FSI.

Furthermore, we compared the reconstructed image obtained via ARSI and FSI using the same number of illuminated patterns.

As shown in Fig. 10, the experimental results demonstrate that ARSI produces images with higher qualities than those produced via FSI when both methods acquired the same number of illuminated patterns.

We conducted another experiment to better compare the reconstructed images obtained via ARSI and normal FSI. The object region, which has a resolution of 98 × 46, was located via ARSI:

Figure 11(a) shows the reconstructed image obtained via ARSI using 9016 patterns. Figure 11(b) presents the reconstructed image obtained via FSI using 128000 patterns. The two experimental results produce nearly the same image quality, although the ARSI method projects considerably less patterns than the normal FSI. When the same number of patterns is used, ARSI produces clearer and sharper images than the normal FSI method, as shown in Fig. 12.

To quantitatively compare the imaging results of ARSI and normal FSI, we use a fragment of white paper with three black bars as the target object. The height and width of the black bars are 3.5 cm and 0.7 cm, respectively. Image quality can be represented by contrast-to-noise ratio (*CNR*) [17], which is defined as

*σ*is the standard deviation of pixel intensity.

The resolution of the object region located via ARSI is 82 × 86, and the object region is located as follows:

The ARSI method projects 14104 patterns to acquire all the Fourier coefficients of the object, whereas the normal FSI method uses 128000 patterns. The experimental result is presented in Figs. 13(a) and 13(b). The *CNR* values of the two images are nearly close. In Fig. 13(c), the grayscale distributions of the two results indicate that image contrast remains stable; ARSI used considerably fewer patterns.

The comparison of the experimental results using the two methods is presented in Fig. 14. In the case where the same number of patterns is projected onto the object, *CNR* and the contrast of the image reconstructed via ARSI are higher than those of the image reconstructed via normal FSI.

The aforementioned experiments illustrate that the number of illuminated patterns is significantly reduced in ARSI to reconstruct an image sharper than that derived from the conventional SPI method. The proposed ARSI can be applied to any situation in which the object region is smaller than the whole illuminated area. The efficiency of SPI can be significantly improved by adaptively locating the object region.

## 4. Discussions and conclusions

We adaptively conduct SPI on the located object region to facilitate imaging efficiency. The proposed ARSI method exhibits the advantages of less projected patterns and higher-quality reconstructed image compared with the conventional SPI method.

The core of the proposed ARSI is locating the object region via the Fourier slice theorem. In practice, only low-frequency vertical and horizontal Fourier patterns are required to project onto the scene, which can smooth the acquired projection line of the scene image and reduce measurement numbers for region location.

Any SPI method is suitable to perform on the located region. Nature images are sparse in the Fourier domain; hence, the FSI method can reconstruct high-quality images with a part of low-frequency coefficients via advanced sampling methods [14].

The positions of detector and digital projector are not overlapped in our experimental set-up, which might cause visible shadows in the reconstructed image, since the object can block the reflected light from some special region. The shadows are small to have little effect on the obtained projection lines, and thus won’t affect the location correctness.

Our proposed method also has limitations. Firstly, the object is required to have a background with uniform grayscale distribution or non-background to guarantee a correct location. Secondly, the proposed ARSI is not applicable for multiple objects in the scene for now, because the location algorithms for one object and separate objects are different. In our future work, the location method for multiple objects will be investigated.

ARSI is experimentally capable of producing high-quality images of an object with less number of measurements than the conventional SPI method. The proposed method may be widely used in areas where the object only takes up a partial region of the scene to increase SPI efficiency. Moreover, if the detector used in the SPI system only observes part of the illuminated region, the proposed method could be applied to reduce the projected pattern number, because the unobserved region can be deemed as the dark background.

## Funding

National Natural Science Foundation of China (NSFC) (61475013); “High-end CNC Machine Tools and Basic Manufacturing Equipment” Special Science and Technology Major (2013ZX04001-011); Program for Changjiang Scholars and Innovative Research Teams in University (IRT1203); Science and Technology Project of Shenzhen Under Grant (JCYJ20140424173418463).

## References and links

**1. **S. M. M. Khamoushi, Y. Nosrati, and S. H. Tavassoli, “Sinusoidal ghost imaging,” Opt. Lett. **40**(15), 3452–3455 (2015). [CrossRef] [PubMed]

**2. **J. H. Shapiro, “Computational ghost imaging,” Phys. Rev. A **78**(6), 061802 (2008). [CrossRef]

**3. **B. Sun, M. P. Edgar, R. Bowman, L. E. Vittert, S. Welsh, A. Bowman, and M. J. Padgett, “3D computational imaging with single-pixel detectors,” Science **340**(6134), 844–847 (2013). [CrossRef] [PubMed]

**4. **Z. Zhang and J. Zhong, “Three-dimensional single-pixel imaging with far fewer measurements than effective image pixels,” Opt. Lett. **41**(11), 2497–2500 (2016). [CrossRef] [PubMed]

**5. **M. J. Sun, M. P. Edgar, G. M. Gibson, B. Sun, N. Radwell, R. Lamb, and M. J. Padgett, “Single-pixel three-dimensional imaging with time-based depth resolution,” Nat. Commun. **7**, 12010 (2016). [CrossRef] [PubMed]

**6. **W. Chan, K. Charan, D. Takhar, K. Kelly, R. Baraniuk, and D. Mittleman, “A single-pixel terahertz imaging system based on compressed sensing,” Appl. Phys. Lett. **93**(12), 121105 (2008). [CrossRef]

**7. **C. Watts, D. Shrekenhamer, J. Montoya, G. Lipworth, J. Hunt, T. Sleasman, S. Krishna, D. Smith, and W. Padilla, “Terahertz compressive imaging with metamaterial spatial light modulators,” Nat. Photonics **8**(8), 605–609 (2014). [CrossRef]

**8. **L. Bian, J. Suo, G. Situ, Z. Li, J. Fan, F. Chen, and Q. Dai, “Multispectral imaging using a single bucket detector,” Sci. Rep. **6**(1), 24752 (2016). [CrossRef] [PubMed]

**9. **S. S. Welsh, M. P. Edgar, R. Bowman, P. Jonathan, B. Sun, and M. J. Padgett, “Fast full-color computational imaging with single-pixel detectors,” Opt. Express **21**(20), 23068–23074 (2013). [CrossRef] [PubMed]

**10. **Z. Zhang, X. Ma, and J. Zhong, “Single-pixel imaging by means of Fourier spectrum acquisition,” Nat. Commun. **6**, 6225 (2015). [CrossRef] [PubMed]

**11. **B.-L. Liu, Z.-H. Yang, X. Liu, and L.-A. Wu, “Coloured computational imaging with single-pixel detectors based on a 2D discrete cosine transform,” J. Mod. Opt. **64**(3), 259–264 (2017). [CrossRef]

**12. **M. P. Edgar, G. M. Gibson, R. W. Bowman, B. Sun, N. Radwell, K. J. Mitchell, S. S. Welsh, and M. J. Padgett, “Simultaneous real-time visible and infrared video with single-pixel detectors,” Sci. Rep. **5**(1), 10669 (2015). [CrossRef] [PubMed]

**13. **E. J. Candes and M. B. Wakin, “An Introduction To Compressive Sampling,” IEEE Signal Process. Mag. **25**(2), 21–30 (2008). [CrossRef]

**14. **L. Bian, J. Suo, X. Hu, F. Chen, and Q. Dai, “Efficient single pixel imaging in Fourier space,” J. Opt. **18**(8), 085704 (2016). [CrossRef]

**15. **R. C. Gonzalez and R. E. Woods, *Digital Image Processing* (Pearson Education, 2011), Chap. 5.

**16. **R. C. Gonzalez and R. E. Woods, *Digital Image Processing* (Pearson Education, 2011), Chap. 10.

**17. **B. Redding, M. A. Choma, and H. Cao, “Speckle-free laser imaging using random laser illumination,” Nat. Photonics **6**(6), 355–359 (2012). [CrossRef] [PubMed]