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

Foveated panoramic ghost imaging

Open Access Open Access

Abstract

Panoramic ghost imaging (PGI) is a novel method by only using a curved mirror to enlarge the field of view (FOV) of ghost imaging (GI) to 360°, making GI a breakthrough in the applications with a wide FOV. However, high-resolution PGI with high efficiency is a serious challenge because of the large amount of data. Therefore, inspired by the variant-resolution retina structure of human eye, a foveated panoramic ghost imaging (FPGI) is proposed to achieve the coexistence of a wide FOV, high resolution and high efficiency on GI by reducing the resolution redundancy, and further to promote the practical applications of GI with a wide FOV. In FPGI system, a flexible variant-resolution annular pattern structure via log-rectilinear transformation and log-polar mapping is proposed to be used for projection, which can allocate the resolution of the region of interest (ROI) and the other region of non-interest (NROI) by setting related parameters in the radial and poloidal directions independently to meet different imaging requirements. In addition, in order to reasonably reduce the resolution redundancy and avoid the loss of the necessary resolution on NROI, the variant-resolution annular pattern structure with a real fovea is further optimized to keep the ROI at any position in the center of 360° FOV by flexibly changing the initial position of the start-stop boundary on the annular pattern structure. The experimental results of the FPGI with one fovea and multiple foveae demonstrate that, compared to the traditional PGI, the proposed FPGI not only can improve the imaging quality on the ROIs with a high resolution and flexibly remain a lower-resolution imaging on the NROI with different required resolution reduction; but also reduce the reconstruction time to improve the imaging efficiency due to the reduction of the resolution redundancy.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Ghost imaging (GI) is an unconventional imaging technique that obtains object information by the correlation calculation between the modulated light pattern and the backscattered light intensity using a bucket detector [15]. GI has the advantages of high detection efficiency, wide spectrum range and low cost. Thus, it has been considered for using in many fields, especially in areas where the traditional imaging methods are technically unavailable or costly, such as infrared imaging [6], X-ray imaging [7], and terahertz imaging [8]. Recently, the proposal of panoramic ghost imaging (PGI) using only a curved mirror [9,10] extends the field of view (FOV) of GI to 360°, making GI a breakthrough in the related applications of wide FOV, such as security monitoring, visual navigation, panoramic situation awareness, etc. Due to the required multiple samples, most conventional GIs are degraded by low resolution or long sampling times [11]. There are already many works on GI optimization to improve the imaging quality with reduced measurements from reconstruction algorithms (such as compress sensing [1215], deep learning [1618]) to measurement matrices [1921]. However, high-resolution GI and particularly high-resolution PGI with a wide FOV of 360° with reduced measurements and less time remains a serious challenge because of the large amount of data. Therefore, it is urgent to develop methods that achieve the coexistence of a wide FOV, and high resolution and high efficiency on GI used for wide FOV and high-resolution applications.

In practical applications of the wide FOV (such as monitoring [22], reconnaissance [23], VR display [24]), there is no need to acquire high-resolution information of all the scenes in the whole FOV. The region of interest (ROI), which is required to be of high resolution, occupies only a part of the whole FOV, while the other region of non-interest (NROI) does not need such a high resolution to provide the required information. Taking the monitoring of sea ship safety as an example, the sea occupies most of the area on the FOV, but we focus on ship information. Then, the imaging of the sea can maintain lower resolution, while the ships, even the ship's personnel, need for the higher-resolution imaging, which can make an efficient monitoring function. In other words, high-resolution imaging of only the ROIs is sufficient to meet the imaging requirements for most applications of wide FOV, and the resolution redundancy on the NROIs can be reduced to improve the imaging efficiency. One effective way to reduce the resolution redundancy is foveated imaging [25]. Inspired by the variant-resolution retina structure of the human eye, foveated imaging has been proposed with high resolution on the fovea area surrounded by the peripheral area with low resolution and has been widely used in many fields [2629]. Recently, foveated imaging has also been introduced into the research on pattern design to enhance the data gathering capacity of traditional GI [3037].

Based on the above observation and our previous work on PGI [9], we report a foveated panoramic ghost imaging (FPGI) to obtain foveated panoramic ghost images with high imaging efficiency by designing variant-resolution annular patterns flexibly and suitably. Specifically, first of all, to better fit the rectangular scenes in real applications, especially panoramic scenes that we focus on, and the ROIs having different lengths and widths, a flexible variant-resolution allocation structure via log-rectilinear transformation [38,39] is proposed, which can allocate the resolution of the ROI and NROI by setting parameters in the horizontal and vertical directions independently to flexibly meet different resolution requirements of GI in practical application. Later, since the pattern used for the projection in the PGI system is designed with a log-polar annular structure to match the convex geometry of the curved mirror [9], the variant-resolution annular pattern structure should be designed via the log-polar mapping from the proposed variant-resolution allocation structure to project on the curved mirror in the FPGI system, as shown in Fig. 1. In addition, since the FOV of the FPGI is 360°, the annular pattern structure is borderless without a fixed start-stop boundary. If the ROI is closed to the traditionally defined boundary, the traditional variant-resolution allocation structure would make the resolution loss of the NROI around the ROI but out of the boundary. Therefore, the variant-resolution allocation structure is further optimized to keep the ROI at any position in the center of the 360° FOV, called as a real fovea, by flexibly changing the initial position of the start-stop boundary of the annular pattern structure. The real fovea can further develop the central advantage of the fovea, especially in the applications where all the areas close to the ROI need to be of high resolution, reasonably reducing the resolution redundancy to avoid the resolution loss. In summary, the FPGI not only improves the imaging quality on the ROI with a high resolution and remain the lower-resolution imaging on the NROI with different required resolution reductions flexibly, but also reduces the reconstruction time due to the reduction of resolution redundancy, which provides an effective way to achieve the coexistence of a wide FOV, high resolution and high efficiency on GI, and further to promote the practical applications of GI with wide FOV.

 figure: Fig. 1.

Fig. 1. Schematic of the FPGI system. Given the designed variant-resolution annular pattern sequence via log-rectilinear transformation, the light modulated by the DMD is projected on the curved mirror through a projecting lens. After the multi-reflection of light between the curved mirror and the omnidirectional scene, the backscattered light intensity is collected by two single-pixel detectors to collect uniformly more and more light intensity from all directions. And then the unwrapping-free foveated panoramic ghost image can be reconstructed by the correlation calculation between the collected light intensity and the variant-resolution rectangular patterns via the log-polar mapping.

Download Full Size | PDF

2. Methodology

2.1 Variant-resolution annular pattern structure via log-rectilinear transformation

The conventional geometric transformation to obtain a retina-like variant-resolution allocation structure on GI is log-polar transformation [30,32], which reduces the resolution of the periphery along only a radial direction with a circular fovea. However, this conventional variant-resolution structure cannot be fully adapted to the rectangular scenes in real applications, especially panoramic scenes. For one thing, the circular geometry of fovea cannot be particularly suitable for the ROI where the length and width differ significantly. For another, the reduction of resolution only along the radial directional cannot meet the imaging requirements of some applications with different imaging resolution around the ROI. Therefore, we propose a variant-resolution allocation structure via log-rectilinear transformation to obtain more flexibility in the resolution allocation to meet different resolution requirements of GI in practical application. The proposed variant-resolution allocation structure via log-rectilinear transformation has several properties. First, the fovea can be designed as a rectangle according to the size of the required ROI to void the sampling redundancy of the extra area designed by the conventional circular fovea. Second, according to the imaging requirements for the NROI, the reduction of the resolution redundancy can be flexibly designed with different parameters in the horizontal and vertical directions independently based on the rectilinear transformation. Third, the resolution decay of the NROI in each direction is still exponential from the center point of the ROI to the sides based on the log transformation.

The variant-resolution allocation structure via log-rectilinear transformation is shown in Figs. 2 (b)-(c). Assume the uniform full resolution of PGI is P × Q, and the required variant resolution of FPGI is U × V. In the p(u)-q(v) coordinate system of one-to-one variant-resolution allocation, given that the central point of the ROI is (p0, q0), the radius of the ROI in the p direction is m_0 and that in the q direction is n_0; the resolution decay factor of the NROI in the p direction is α1 and that in the q direction is α2. Then, the model of the proposed variant-resolution allocation via log-rectilinear transformation can be expressed as:

$$\left\{ \begin{array}{l} R_1^u = \left\{ \begin{array}{l} u,\quad \quad \quad \quad \quad \quad \quad \quad \;\;\;\;\;\;\;when\;1 \le u \le m\_0,\;\;\;\;\;\;\;\;\;\;\;\;\\ m\_0 \cdot \alpha_1^{u - m\_0},\quad \quad \quad \quad \quad \;\,\,\,when\;m\_0 \le u \le {{(U - 1)} / 2},\; \end{array} \right.\\ R_2^v = \left\{ \begin{array}{l} v,\quad \quad \quad \quad \quad \quad \quad \quad \quad \,\,\,\,\,\;when\;1 \le v \le n\_0,\;\;\;\;\;\;\;\;\;\;\;\;\\ n\_0 \cdot \alpha_2^{v - n\_0},\quad \quad \quad \quad \quad \quad \,when\;n\_0 \le v \le {{(V - 1)} / 2},\; \end{array} \right.\\ p(u) = \left\{ \begin{array}{l} \max (1,round({p_0} - R_1^u)),\,\,\;\;\;\;when\;1 \le u \le {{(U - 1)} / 2},\;\;\;\,\\ {p_0},\quad \quad \quad \quad \quad \quad \quad \quad \,\quad \;when\;u = {{(U + 1)} / 2},\;\;\;\;\;\;\;\,\\ \min (round({p_0} + R_1^u),P),\,\,\;\;\;\;when\;{{(U + 1)} / {2 < u \le U}}, \end{array} \right.\\ q(v) = \left\{ \begin{array}{l} \max (1,round({q_0} - R_2^v)),\,\;\,\quad when\;1 \le v \le {{(V - 1)} / 2},\quad \\ {q_0},\quad \quad \quad \quad \quad \quad \quad \quad \;\,\quad when\;v = {{(V + 1)} / 2},\\ \min (round({q_0} + R_2^v),Q),\;\;\quad when{{\;\;(V + 1)} / {2 < v \le V}}, \end{array} \right.\\ U \in 2N + 1,V \in 2N + 1 \end{array} \right.$$
where $R_1^u$ represents the radius of the uth layer extending from the central point of the ROI in the p direction and $R_2^v$ represents the radius of the vth layer extending from the central point of the ROI in the q direction. p(u) provides the uth value of variant-resolution allocation in the p direction mapped in the uniform-resolution pattern structure, and q(v) provides the vth value of variant-resolution allocation in the q direction mapped in the uniform-resolution pattern structure. Note that, given a fixed ROI, the larger the value of α1 or α2 is, the smaller the corresponding values of U or V is, and the more resolution redundancy of the NROI is reduced to improve the imaging efficiency of the FPGI.

 figure: Fig. 2.

Fig. 2. Schematic of the variant-resolution annular pattern structure via log-rectilinear transformation. (a) Uniform-resolution rectangular pattern structure in log-polar coordinate with P × Q pixels. (b) Variant-resolution allocation structure with U × V pixels. (c) Detail of the one-to-one variant-resolution allocation structure in the p(u)-q(v) coordinate system. (d) Uniform-resolution annular pattern structure via log-polar mapping. (e) Variant-resolution annular pattern structure via log-polar mapping. (f) Detail of part of the variant-resolution annular structure.

Download Full Size | PDF

Furthermore, to match the convex geometry of the curved mirror used in the FPGI system shown in Fig. 1, the previous pattern used for the projection is designed with a log-polar annular structure [9]. The variant-resolution pattern used for the projection in the FPGI system is also designed with a log-polar annular structure accordingly. Therefore, based on the above variant-resolution allocation structure, the variant-resolution annular pattern structure in the FPGI with U rings and V non-uniform cells in each ring can be obtained via log-polar mapping, as shown in Fig. 2 (d)-(f). Combining Eq. (1), as p(u) represents the uth ring in the variant-resolution annular pattern structure mapped in the uniform-resolution annular pattern structure with P rings and q(v) represents the vth cell in each ring, and the variant-resolution annular pattern structure can be obtained with the following model:

$$\left\{ \begin{array}{l} \varepsilon = \frac{{1\textrm{ + }\sin ({\pi /Q} )}}{{1 - \sin ({\pi /Q} )}}\\ r{c_1} = \frac{{{r_0}}}{{1 - \sin ({\pi /Q} )}}\\ {r_{p(u)}} = {r_0} \cdot {\varepsilon^{p(u)}}\quad \\ r{c_{p(u)}} = r{c_1} \cdot {\varepsilon^{p(u)\textrm{ - 1}}}\\ {\theta_{q(v)}} = q(v) \cdot \frac{{2\pi }}{Q}\textrm{ } \qquad (v = 1,2,3\ldots V)\\ {\xi_{p(u)}} = {\log_\varepsilon }(r{c_{p(u)}}) = {\log_\varepsilon }(r{c_1}) + p(u) - 1\textrm{ } \quad(u = 1,2,3\ldots U)\; \end{array} \right.,$$
where r0 is the radius of the fovea blind hole in the annular pattern structure used to cover the center of the curved mirror with the self-reflection, ɛ is the increasing coefficient of the radius between adjacent rings in the uniform-resolution annular pattern structure, rp(u) is the outer radius of the uth ring in the variant-resolution annular pattern structure, rcp(u) is the center radius of the uth ring, θq(v) is the angle of the vth cell of each ring starting at 0, and ξp(u) is the mapped value of rcp(u) in the log-polar coordinate.

The above variant-resolution annular pattern structure is based on the design idea of the traditional foveated structure, which has a fixed start-stop boundary of FOV. However, the annular 360° FOV of the FPGI is borderless. Given that the designed annular pattern always starts at 0 in the angular direction of θ, there is an imaginary boundary from 0 to 2π. If the angular position of ROI is near 2π, in the log-polar coordinate, the areas of NROI near 0 would be the one with the farthest distances from ROI. With the variant-resolution allocation structure with traditional fovea, the area near 0 is allocated the lowest resolution to obtain the least information, filled in green in Fig. 3 (a). In the Cartesian coordinate system, however, the areas near 0 surrounding the ROI are also close to the boundary 2π, because the cyclic period is 2π, as shown in Fig. 3 (b), which should provide a higher resolution to obtain more information closely around the ROI. Therefore, considering the arbitrariness of the required ROI position, the variant-resolution allocation can be optimized by flexibly changing the initial position in the log-polar coordinate to center the ROI, as shown in Fig. 3 (c). Since the ROI is always kept in the center of the whole FOV flexibly, the optimized fovea is called as a real fovea. Based on the rotation invariance of the log-polar transformation [40], if the abscissa of the log-polar coordinate system shifts horizontally, the position of the Cartesian coordinate system will be rotated (anticlockwise is positive). Given that the central point of the required ROI is (p0, q0), the central point of the whole FOV is (pc, qc), and the corresponding rotated angle φ0 to center the ROI can be calculated as follows:

$$\left\{ \begin{array}{l} {\varphi_0} = ({{q_c} - {q_0}} )\cdot \frac{{2\pi }}{Q}\\ \theta^{\prime} = \theta + {\varphi_0} \end{array} \right.,$$
where θ represents the angle in the initial coordinate system and θ’ represents the angle in the rotated coordinate system where the ROI is in the center of the whole FOV. Then, the variant-resolution allocation can be optimized, as shown in Fig. 3 (e), where the area near 0 is given a higher-resolution allocation than that on the traditional allocation structure in Fig. 3 (a). The variant-resolution annular structure with a real fovea shown in Fig. 3 (f) can develop the central advantage of the fovea, reasonably reducing the resolution redundancy to avoid the resolution loss, especially in the applications where all the areas close to the ROI need to be of a high resolution.

 figure: Fig. 3.

Fig. 3. Schematic of variant-resolution annular pattern structure with real fovea. (a) The variant-resolution rectangular structure with traditional fovea in the log-polar coordinate, where the ROI near 2π is at the edge and the area in green is the farthest from the ROI. (b) The variant-resolution annular structure with traditional fovea in the log-polar coordinate, where the area in green is close to the ROI. (c) Variant-resolution structure with the real fovea via the horizontal shifting of the abscissa of the log-polar coordinate system to center the ROI in the whole FOV. (d) Variant-resolution annular structure with the real fovea with a rotated Cartesian coordinate system. (e) After optimization with the real fovea, the variant-resolution allocation structure in the original log-polar coordinate, where the area in green has a higher-resolution allocation than the one in (a). (f) Variant-resolution annular pattern structures with the real fovea in the original Cartesian coordinate system. Note that ξ'-θ’ is the shifted coordinate to center the ROI in the log-polar coordinate system from φ0 to φ0 + 2π and X'-Y’ is the rotated coordinate axis in the Cartesian coordinate system with the rotated angle φ0.

Download Full Size | PDF

2.2 Unwrapping-free foveated panoramic ghost imaging

Based on our previous work on unwrapping-free PGI [9], Fig. 1 illustrates a schematic of the FPGI system. Given the imaging requirements of ROI and NROI in the omnidirectional scene, a variant-resolution annular pattern can be obtained according to the design method in Section 2.1. With the variant-resolution annular pattern projected onto the curved mirror, the back-scattered light intensity after the multi-reflection between the curved mirror and the omnidirectional scene, denoted by I, can be written as follows:

$$I = \sum\limits_x {\sum\limits_y {{P_{v - annular}}(x,y) \cdot {O_{omnidirectional}}(x,y)} } ,$$
where x and y index the position of the omnidirectional scene in the Cartesian coordinate system, Pv-annular represents a variant-resolution annular pattern with U rings and V non-uniform cells of each ring, and Oomnidirectionalrepresents the omnidirectional scene reflected on the curved mirror.

Although the variant-resolution allocation structure is designed with U × V pixels, the variant-resolution pattern is created by reformatting into the uniform high-resolution grid with P × Q pixels to ensure high-resolution imaging on the ROI. Therefore, given the two transformation matrices in the p direction and the q direction respectively, the variant-resolution rectangular pattern via log-polar mapping can be obtained as follows:

$${P_{v - rectangular}} = {A_p}{P_n}A_q^T,$$
where Pn is the variant-resolution allocation structure matrix with U × V, Pv-rectangular is the variant-resolution rectangular pattern reformatted in the high-resolution grid with P × Q, Ap is a P × U transformation matrix mapping the U-element variant-resolution structure to the P-element high-resolution grid in the p direction, and Aq is a Q × V transformation matrix mapping the V-element variant-resolution structure to the Q-element high-resolution grid in the q direction. Ap and Aq are binary matrices where the locations of the ‘ones’ in Column n denote the high-resolution pixels that belong to the nth variant-resolution structure in the corresponding direction.

According to the previously proposed unwrapping-free reconstruction of PGI, the variant-resolution panoramic ghost image can be directly reconstructed by the correlation between the collected light intensity using the projected variant-resolution annular pattern and the variant-resolution rectangular pattern via log-polar mapping. To obtain better imaging quality with fewer measurements, an efficiency compress sensing reconstruction algorithm of total variation (TV) is used on the unwrapping-free reconstruction of the FPGI.

$$\begin{aligned} \min \quad \textrm{ }&{||c ||_{{l_1}}}\\ s.t.\textrm{ }\quad&H{O_{v - panoramic}} = c\quad \quad \quad ,\\& \textrm{ }{{\boldsymbol P}_{v - rectangular}}{O_{v - panoramic}} = {\boldsymbol I} \end{aligned}$$
where H represents the gradient calculation matrix, c is the corresponding coefficient vector, l1 represents the l1 norm, IRT × 1 denotes the measurement matrix of the projected variant-resolution annular patterns (T is the number of measurements), Pv-rectangularRT × M represents the variant-resolution rectangular pattern sequence matrix via log-polar mapping (M = P × Q is the pixel number of the variant-resolution pattern on the high-resolution grid), and Ov-panoramic represents the required variant-resolution panoramic ghost image of the omnidirectional scene.

3. Results

3.1 Experimental setup

The experimental setup of the FPGI system is the same as that used in our previous work on PGI system [9]. Based on the schematic of FPGI system shown in Fig. 1, the light emitted from a light emitting diode (LED) operating at 500-700nm (@20W) is modulated with the proposed variant-resolution annular pattern, which is designed based on the imaging requirement and generated by a DMD (Texas Instruments DLP4100; refresh rate, 10KHz). Through a projecting lens (focus length, 150 mm), the modulated light illuminates the curved mirror made of a plane-convex lens coated with a strong reflective film (diameter, 50.8 mm; central thickness, 16.3 mm; radius of curvature, 30.91 mm), and then generates the multi-reflection between the curved mirror and the omnidirectional scene. Later, the backscattered light intensity is collected by two photodetectors (Thorlabs PDA36A; active area, 13 mm2) to collect uniformly more and more light intensity from all directions and converted to a digital signal by a data acquisition board (DAQ, PICO6404E). Finally, the unwrapping-free foveated panoramic ghost image is reconstructed on the computer (6226R CPU × 2, 128GB RAM). Note that, in order to observe the omnidirectional scene conveniently, the object is attached to a hollow cylinder (height, 50 mm; base diameter, 75 mm).

Moreover, since the FPGI focuses on improving the imaging quality of the required ROI on the wide FOV by reducing the resolution redundancy of the NROI, the two classic evaluation indexes to analyze the imaging quality, peak signal to noise ratio (PSNR) and structural similarity index measure (SSIM), are used to compare quantitatively the difference on the ROI between the uniform-resolution PGI and variant-resolution FPGI, which are defined as follows:

$$\left\{ \begin{aligned} &{\text{PSNR}} = 10{\log _{10}}\frac{{{{\left( {{2^k} - 1} \right)}^2}}}{{{\text{MSE}}}} \hfill \\ &{\text{MSE}} = \frac{1}{S}\sum\limits_{i,j} {{{\left( {F'(i,j) - F(i,j)} \right)}^2}} \; \hfill \\ &\;{\text{SSIM}}{{\text{M}}_{i,j}} = \frac{{\left( {2\mu \mu ' + {c_1}} \right)\left( {2w + {c_2}} \right)}}{{\left( {{\mu ^2} + u{'^2} + {c_1}} \right)\left( {{\sigma ^2} + \sigma {'^2} + {c_2}} \right)}} \hfill \\ \end{aligned} \right.,$$
where F(i, j) is the ground truth of the ROI, F'(i, j) is the fovea in the reconstructed image, MSE represents the mean square error; S is the total pixels of the fovea; k is the number of bits and set as 8. µ and µ’ are the average value of F(i, j) and F'(i, j); σ and σ’ are the variance of F(i, j) and F'(i, j). w is the covariance between F(i, j) and F'(i, j). c1 = (k1×L)2 and c2 = (k2×L)2 are the constant with k1 = 0.01, k2 = 0.03 and L = 1. Note that, the higher the PSNR, the better the imaging quality of the fovea. In addition, the larger the value of SSIM is from 0 to 1, the reconstructed fovea is closer to the object of ROI.

3.2 FPGI with variant resolution

With the experimental setup of the FPGI system, FPGI experiments with variant-resolution are conducted. First, the simple object is a string of letters “PANORAMIC” attached to a hollow cylinder to make an omnidirectional 360° FOV, as shown in Fig. 4 (a). Accordingly, the uniform-resolution annular pattern structure is designed with 20 rings and 102 cells in each ring, as shown in Fig. 4 (b). To observe the effect of the FPGI with different reductions in the resolution redundancy, the three sets of parameters on the variant-resolution pattern structure can be seen in Table 1. Given the position of the required ROI “ORA”, highlighted by a red box in Fig. 4(a), and the different resolution decay factors of NROI in the p direction and q direction, the variant-resolution annular pattern structures via log-rectilinear transformation can be designed based on Eq. (1) and Eq. (2), and the corresponding pattern samples randomly filled by 0 or 1 are given in Fig. 4 (c)-(e). As seen in Table 1, as the resolution decay factors of the NROI increase, the total pixels of the three variant-resolution pattern structures decrease in the corresponding direction, resulting in a greater reduction in the resolution redundancy on the NROI. For example, for ‘Variant-resolution 2’, the resolution decay factor in the q direction α2 increases from 1.08 to 1.16, and the pixel number of the variant-resolution pattern in the q direction is reduced from 66 to 50.

 figure: Fig. 4.

Fig. 4. FPGI with variant resolution. (a) A sample object to make an omnidirectional 360° FOV, where the ROI is highlighted by the red box. (b) A uniform-resolution annular pattern sample randomly filled by 0 or 1 with 20 × 102 pixels and the corresponding uniform-resolution rectangular pattern via log-polar mapping. (c) With the parameters of ‘Variant-resolution 1’ in Table 1, a random variant-resolution annular pattern sample with 19 × 66 pixels and the log-polar mapped variant-resolution rectangular pattern. (d) With the parameters of ‘Variant-resolution 2’ in Table 1, a random variant-resolution annular pattern sample with 19 × 50 pixels and the log-polar mapped variant-resolution rectangular pattern. (e) With the parameters of ‘Variant-resolution 3’ in Table 1, a random variant-resolution annular pattern sample with 16 × 38 pixels and the log-polar mapped variant-resolution rectangular pattern. (f) With the same measurements of full sampling, the unwrapping-free reconstructed panoramic ghost images using the four annular pattern sequences with different structures of (b)-(e) respectively, and the corresponding magnified views of the ROIs. (g) Quantitative comparison on the PSNR of fovea of the reconstructed image using a uniform-resolution annular pattern structure and the three variant-resolution pattern structures with various parameters shown in (b)-(e) at different sampling rates. (h) Quantitative comparison on the reconstruction time of unwrapping-free PGI using the four different pattern structures at different sampling rates. Note that, the sampling ratio refers to the uniform-resolution panoramic ghost image with 20 × 102 pixels, because the variant-resolution panoramic ghost image is also reformatted into the uniform-resolution grid to ensure high-resolution imaging on the ROI.

Download Full Size | PDF

Tables Icon

Table 1. Parameters of the four pattern structures in Fig. 4

Using the four annular pattern sequences with different parameters given in Table 1, the corresponding unwrapping-free panoramic ghost images were reconstructed. Figure 4 (f) shows the four reconstructed images with the same measurements of full sampling and the corresponding magnified views of the ROIs. According to the observation, both subjectively and quantitatively, the imaging quality of the three foveated variant-resolution panoramic images is better than that of the uniform-resolution panoramic image with noisy background. In addition, the quantitative curves of the PSNR at different sampling rates are given in Fig. 4 (g) to further analyze the imaging quality using the three variant-resolution annular pattern structure. With the greater reduction in the resolution redundancy of NROI from ‘variant-resolution 1’ to ‘variant-resolution 3’, the PSNR of the fovea is mainly higher at the same sampling ratio, because the obtained redundant information on the NROI is less; however, the resolution on the NROI is too low to be distinguished, especially in ‘variant-resolution 3’. Therefore, it is necessary to balance the PSNR of the ROI and the required information on the NROI by setting the parameters of the variant-resolution pattern structure flexibly in accordance with the various imaging requirements in different situations. Moreover, due to the reduction in the total resolution, the reconstruction time on the FPGI is always less than that on the traditional PGI using the same unwrapping-free reconstruction algorithm, as shown in Fig. 4 (h). In particular, at a sampling rate of 1, all the reconstruction time of the three variant-resolution foveated panoramic ghost images are closed to 0.4 s, while the reconstruction time of the traditional uniform-resolution panoramic ghost image is larger than 1.2 s, which is 3 times more than that of the image used on the FPGI.

A complex scene is also used to conduct the experiments to observe the effect of the FPGI once more, as shown in Fig. 5. Here, the uniform-resolution annular pattern structure is designed with 40 rings and 204 cells in each ring to obtain higher-resolution images. The required ROI is highlighted by the red box in Fig. 5 (a), and the random pattern samples with different pattern structures designed via the parameters in Table 2 are shown in Figs. 5 (b)-(e) respectively.

 figure: Fig. 5.

Fig. 5. FPGI with variant resolution. (a) A complex object to make an omnidirectional 360° FOV, where the ROI is highlighted by the red box. (b) A uniform-resolution annular pattern sample randomly filled by 0 or 1 with 40 × 204 pixels and the log-polar mapped uniform-resolution rectangular pattern. (c) With the parameters of ‘Variant-resolution 1’ in Table 2, a random variant-resolution annular pattern sample with 29 × 98 pixels and the log-polar mapped variant-resolution rectangular pattern. (d) With the parameters of ‘Variant-resolution 2’ in Table 2, a random variant-resolution annular pattern sample with 20 × 98 pixels and the log-polar mapped variant-resolution rectangular pattern. (e) With the parameters of ‘Variant-resolution 3’ in Table 2, a random variant-resolution annular pattern sample with 20 × 50 pixels and the log-polar mapped variant-resolution rectangular pattern. (f) With the same measurements of full sampling, the unwrapping-free reconstructed panoramic ghost images with the four annular pattern sequences with different structures of (b)-(e) respectively, and the corresponding magnified views of the ROIs. (g) Quantitative comparison on the SSIM of fovea of the reconstructed image using a uniform-resolution annular pattern structure and the three variant-resolution pattern structures with different parameters shown in (b)-(e) at different sampling rates. (h) Quantitative comparison on the reconstruction time of PGI with the four different pattern structures at different sampling rates. Note that, the sampling ratio refers to the uniform-resolution panoramic ghost image with 40 × 204 pixels.

Download Full Size | PDF

Tables Icon

Table 2. Parameters of the four pattern structures in Fig. 5

Figure 5 (f) shows the four unwrapping-free reconstructed panoramic ghost images with the same measurements at the full sampling rate, which use the uniform-resolution annular pattern sequence and the three variant-resolution annular pattern sequences with different parameters given in Table 2 respectively. It is clear that the uniform-resolution panoramic ghost image is much noisier, and the letter ‘B’ on the ROI is flat. In contrast, the backgrounds of the three variant-resolution panoramic ghost images are clean, and the letter ‘B’ on the ROI is structurally fuller. In addition, the three variant-resolution foveated panoramic ghost images also show the flexible acquisition on the NROI according to the imaging requirements in the horizontal and vertical directions independently by the variant-resolution allocation via log-rectilinear transformation. If the NROI around the ROI in the vertical direction is required to obtain major information, such as ‘Variant-resolution 1’, the resolution decay factor in the p direction α1 is set with a small value of 1.16. Conversely, such as ‘Variant-resolution 2’, α1 is set with a larger value of 1.69 to reduce the resolution redundancy in the vertical direction that is not needed, while α2 is set with a small value of 1.06 to independently obtain the required information on the NROI around the ROI in the horizontal direction.

Moreover, to quantitatively analyze the structural information of the ROI, the SSIM of fovea at different sampling rates are given in Fig. 5 (g). First, the SSIM of fovea of the uniform-resolution PGI is always lower than that of the variant-resolution FPGI, which proves that the imaging quality of the ROI in our proposed FPGI is always better than that of the traditional PGI to achieve GI with high imaging quality and wide FOV. In addition, from the trends of the three curves on the FPGI with different variant-resolution pattern structures, the SSIM of ‘variant-resolution 3’ increases most rapidly and then remains steady with increasing sampling rate, where the reduction in the resolution redundancy of NROI on ‘variant-resolution 3’ is the greatest with the fewest total pixels shown in Table 2. For example, for the SSIM of fovea value reaching 0.7, the SSIM of fovea on ‘variant-resolution 1’ reaches the value at the sampling rate of 0.6, the SSIM of fovea on ‘variant-resolution 2’ reaches the value at the sampling rate of 0.5, while the SSIM on ‘variant-resolution 3’ can reach the value at the sampling rate of 0.3. The lower sampling rate means the fewer measurements, therefore, if fewer measurements are required to obtain higher-quality imaging on the ROI, the obtained information on the NROI should be less with a greater reduction in the resolution redundancy. With the same conclusion as the one on the first experiments using a sample object, it is necessary to balance the imaging quality of the ROI and the required information on the NROI by setting the parameters of the variant-resolution pattern structure flexibly in accordance with the various imaging requirements in different situations. Moreover, the obvious difference in the reconstruction time in Fig. 5 (h) shows that the FPGI can reduce the imaging time by reducing the resolution redundancy to improve the imaging efficiency of PGI. Overall, the above experiments show that the proposed FPGI can achieve the coexistence of a wide FOV, high quality and high efficiency on GI to promote the practical applications of GI with a wide FOV.

3.3 FPGI with real fovea

The ROIs in the above experiments are chosen at the center of the whole 360° FOV in the original start-stop boundary of [0, 2π]. However, in the practical applications, the ROI can be in any position of the whole FOV according to the imaging requirements. If the ROI is in the position close to the original start-stop boundary, such as the little engine framed in red in Fig. 6 (a), with the traditional-fovea variant-resolution annular structure, the resolution of the areas circled in green in Fig. 6 (c) would be very low to lose the information, because the areas are supposed to be the farthest from the ROI in the original start-stop boundary from 0 to 2π. However, the areas are close to the ROI in the actual position because it is an annular structure with a period 2π. Therefore, in order to keep the ROI in the center of the whole FOV, the variant-resolution annular pattern structure with a real fovea is optimized to avoid the resolution loss of the areas close to the ROI by flexibly changing the initial position of the start-stop boundary. Given that the central point of the ROI is (19,26), the central point of the whole FOV is (21,103); according to the Eq. (3), the optimized variant-resolution annular pattern structure of this ROI can be obtained, and one random annular pattern sample is shown in Fig. 6 (c).

 figure: Fig. 6.

Fig. 6. FPGI with a real fovea. (a) The object with an ROI located in the position close to the original start boundary 0. (b) A uniform-resolution annular pattern sample randomly filled by 0 or 1 with 40 × 204 pixels. (c) A traditional-fovea variant-resolution annular pattern sample with 39 × 86 pixels, where the resolution of the areas circled in green is lost. (d) A real-fovea variant-resolution pattern sample with 39 × 101 pixels, which avoids the resolution loss of the areas circles in green. (e) The unwrapping-free uniform-resolution panoramic ghost image with full sampling in the original coordinate and the corresponding image where the ROI is at the center. (f) The traditional-fovea variant-resolution panoramic ghost image with full sampling in the original coordinate and the corresponding image with centering the ROI, where the areas marked by the green box are not as clear as the areas marked by the blue box (Quantitatively, the SSIM of the green box is 0.35 while the SSIM of the blue box is 0.61). (g) The real-fovea variant-resolution panoramic ghost image with full sampling and the corresponding image with centering the ROI, where the areas marked by the green box are as clear as the areas marked by the blue box (Quantitatively, the SSIM of the green box is improved to 0.60, which is close to the SSIM of the blue box of 0.61).

Download Full Size | PDF

Using the annular pattern sequences shown in Figs. 6 (b)-(d), the reconstructed panoramic ghost images in the original log-polar coordinate from 0 to 2π and the corresponding image in the changed log-polar coordinate centering the ROI are given in Figs. 6 (e)-(g). First, in the original log-polar coordinate, the areas marked by the blue box are visually very far from the ROI; however, the distance to the ROI is actually the same as the distance between the areas marked by the blue box with the ROI in the omnidirectional FOV. With the traditional-fovea variant-resolution annular patterns, as shown in Fig. 6 (f), the areas marked by the green box are not as clear as the areas marked by the blue box which makes the resolution loss close around the ROI, especially in applications that require simultaneous attention to the information closely surrounding the ROI. With the real-fovea variant-resolution annular patterns, as shown in Fig. 6 (g), the areas marked by the green box are as clear as those marked by the blue box. The proposed variant-resolution annular pattern structure with a real fovea can be proven to reasonably reduce the resolution redundancy to avoid the resolution loss.

3.4 FPGI with multiple foveae

In practical applications, there may be more than one required ROI in the wide FOV of 360°, and the variant-resolution annular pattern structure of FPGI needs to be designed with multiple foveae to meet the imaging needs. Given the required multiple ROIs, the corresponding variant-resolution allocation structures with individual fovea can be designed independently according to the imaging requirements. Then, the whole NROI is divided into parts for surrounding each ROI based on the position relation with the required ROIs. Finally, the variant-resolution annular pattern structure with multiple foveae can be obtained by merging the corresponding parts on each variant-resolution allocation structure with the individual fovea.

Here, we present the two experimental examples to show the FPGI with multiple foveae: one with dual foveae and the other with treble foveae, as shown in Fig. 7. As seen from the comparison on the imaging quality and the reconstruction time between the reconstructed uniform-resolution panoramic ghost images and the variant-resolution panoramic ghost images with multiple foveae, although there are multiple ROIs, the FPGI is still of higher quality on each ROI than the one on the traditional PGI and retains the lower-resolution imaging on the surrounding NROIs; and the reconstruction time of the FPGI with multiple foveae is still less than that of the traditional uniform-resolution PGI, especially with the increase of sampling rate, the difference of the reconstruction time is more obvious. Certainly, the reduction of resolution redundancy on the NROIs for different requirements can be achieved by flexibly designing different variant-resolution annular pattern structures with multiple foveae.

 figure: Fig. 7.

Fig. 7. FPGI with multiple foveae. (a) The object with the two ROIs marked in red and green respectively. (b) A random variant-resolution annular pattern sample with double foveae. (c) Reconstructed uniform-resolution panoramic ghost image with full sampling and the magnified views of the two ROIs following the quantitative evaluation of the PSNR and SSIM. (d) Reconstructed variant-resolution panoramic ghost image using the patterns shown in (b) with full sampling and the magnified views of the two ROIs following the quantitative evaluation of the PSNR and SSIM. (e) The object with the three ROIs marked in red, green and blue respectively. (f) A random variant-resolution annular pattern sample with treble foveae. (g) Reconstructed uniform-resolution panoramic ghost image with full sampling and magnified views of the three ROIs following the quantitative evaluation of the PSNR and SSIM. (h) The reconstructed variant-resolution panoramic ghost image using the patterns shown in (f) with full sampling and magnified views of the three ROIs following the quantitative evaluation of the PSNR and SSIM. (i) Comparison on the reconstruction time between the traditional uniform-resolution PGI and the FPGI with the two ROIs at different sampling rates. (j) Comparison on the reconstruction time between the traditional uniform-resolution PGI and the FPGI with the three ROIs at different sampling rates. Note that, the sampling ratio refers to the uniform-resolution panoramic ghost image with 40 × 204 pixels.

Download Full Size | PDF

4. Discussion and conclusion

Here, we report an efficient panoramic ghost imaging method called the FPGI, which provides a way to achieve the coexistence of a wide FOV, high resolution and high efficiency on GI. In the FPGI, the high-resolution imaging is only on the ROIs of the whole FOV, while the NROIs are imaged with low resolution to reduce resolution redundancy to improve the imaging efficiency in the 360° wide FOV. To achieve the FPGI, we perform some work in this paper. First, we propose a variant-resolution allocation structure via log-rectilinear transformation to flexibly meet different resolution requirements of GI, which can allocate the resolution of ROIs and NROIs by setting parameters in the horizontal and vertical directions independently. Then, the variant-resolution annular pattern structure used for projection on FPGI is designed via log-polar mapping from the proposed variant-resolution allocation structure. Next, in order to reasonably reduce the resolution redundancy and avoid the resolution loss, the real-fovea variant-resolution allocation structure is further optimized to keep the ROI at any position in the center of the 360° FOV by flexibly changing the initial position of the start-stop boundary of the annular pattern structure. In addition, we also describe the design of an FPGI with multiple foveae to meet imaging needs with multiple ROIs in practical applications with a wide FOV. The experimental results demonstrate that, compared to the traditional PGI, the proposed FPGI can improve the imaging quality on the ROIs with a high resolution and flexibly remain the lower-resolution imaging on the NROI with different required resolution reductions in the horizonal and vertical direction. Simultaneously, the FPGI can reduce the reconstruction time to improve the imaging efficiency due to the reduction of the resolution redundancy.

Furthermore, the proposed FPGI can potentially be further enhanced by combining existing interesting ideas and mature technologies. For example, with deep learning, first, the pattern optimization with a priori knowledge of the ROIs can be achieved by training with neural networks, such as a deep Q-learning network [17]. Second, the intelligent selection of the ROIs can be achieved by using an object detection algorithm with deep learning, such as generative adversarial nets [34]. Third, the end-to-end deep-learning approach [16] can be introduced to directly correspond to each pixel of the variable-resolution structure to reduce measurements and reduce the imaging time. Besides, we have demonstrated the proposed FPGI in static scenes, an efficient and high-quality FPGI in a dynamic scene can be obtained in the near future by further optimizing patterns and reconstruction algorithms to enhance the functions of FPGI system to open up possibilities for more applications with a wide FOV, such as the motion tracking with an adaptive fovea [30] in a wide FOV. Moreover, the proposed variant-resolution allocation structure and the designed variant-resolution annular structure are not only for achieving the FPGI, but also provide new ideas and references for the practical process of GI and the applications on special scenes, such as circular FOV.

Funding

Beijing Municipal Natural Science Foundation (4222017); Ensan Foundation (2022011); National Natural Science Foundation of China (62275017, 62275022); Funding of foundation enhancement program under Grant (2019-JCJQ-JJ-273).

Acknowledgments

The authors thank the editor and the anonymous reviewers for their valuable suggestions.

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. H. Shapiro, “Computational ghost imaging,” Phys. Rev. A 78(6), 061802 (2008). [CrossRef]  

2. M. P. Edgar, G. M. Gibson, and M. J. Padgett, “Principles and prospects for single-pixel imaging,” Nat. Photonics 13(1), 13–20 (2019). [CrossRef]  

3. R. S. Bennink, S. J. Bentley, and R. W. Boyd, ““Two-Photon” Coincidence Imaging with a Classical Source,” Phys. Rev. Lett. 89(11), 113601 (2002). [CrossRef]  

4. B. I. Erkmen and J. H. Shapiro, “Ghost imaging: from quantum to classical to computational,” Adv. Opt. Photon. 2(4), 405 (2010). [CrossRef]  

5. 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(5), R3429–R3432 (1995). [CrossRef]  

6. R. S. Aspden, N. R. Gemmell, P. A. Morris, D. S. Tasca, L. Mertens, M. G. Tanner, R. A. Kirkwood, A. Ruggeri, A. Tosi, R. W. Boyd, G. S. Buller, R. H. Hadfield, and M. J. Padgett, “Photon-sparse microscopy: visible light imaging using infrared illumination,” Optica 2(12), 1049 (2015). [CrossRef]  

7. Y. Klein, A. Schori, I. P. Dolbnya, K. Sawhney, and S. Shwartz, “X-ray computational ghost imaging with single-pixel detector,” Opt. Express 27(3), 3284 (2019). [CrossRef]  

8. L. Olivieri, J. S. T. Gongora, L. Peters, V. Cecconi, A. Cutrona, J. Tunesi, R. Tucker, A. Pasquazi, and M. Peccianti, “Hyperspectral terahertz microscopy via nonlinear ghost imaging,” Optica 7(2), 186 (2020). [CrossRef]  

9. H. Cui, J. Cao, Q. Hao, D. Zhou, M. Tang, K. Zhang, and Y. Zhang, “Omnidirectional ghost imaging system and unwrapping-free panoramic ghost imaging,” Opt. Lett. 46(22), 5611–5614 (2021). [CrossRef]  

10. Z. Ye, H. B. Wang, J. Xiong, and K. Wang, “Ghost panorama using a convex mirror,” Opt. Lett. 46(21), 5389–5392 (2021). [CrossRef]  

11. D. Li, Z. Gao, and L. Bian, “Efficient large-scale single-pixel imaging,” Opt. Lett. 47(21), 5461 (2022). [CrossRef]  

12. O. Katz, Y. Bromberg, and Y. Silberberg, “Compressive ghost imaging,” Appl. Phys. Lett. 95(13), 131110 (2009). [CrossRef]  

13. C. Zhou, X. Liu, Y. Feng, X. Li, G. Wang, H. Sun, H. Huang, and L. Song, “Real-time physical compression computational ghost imaging based on array spatial light field modulation and deep learning,” Opt. Laser Eng. 156, 107101 (2022). [CrossRef]  

14. M. F. Duarte, M. A. Davenport, D. Takhar, J. N. Laska, T. Sun, K. F. Kelly, and R. G. Baraniuk, “Single-pixel imaging via compressive sampling,” IEEE Signal Proc. Mag. 25(2), 83–91 (2008). [CrossRef]  

15. W. Gong and S. Han, “High-resolution far-field ghost imaging via sparsity constraint,” Sci. Rep. 5(1), 9280 (2015). [CrossRef]  

16. F. Wang, H. Wang, H. Wang, G. Li, and G. Situ, “Learning from simulation: An end-to-end deep-learning approach for computational ghost imaging,” Opt. Express 27(18), 25560 (2019). [CrossRef]  

17. Z. Wang, W. Zhao, A. Zhai, P. He, and D. Wang, “DQN based single-pixel imaging,” Opt. Express 29(10), 15463–15477 (2021). [CrossRef]  

18. H. Wu, R. Wang, G. Zhao, H. Xiao, D. Wang, J. Liang, X. Tian, L. Cheng, and X. Zhang, “Sub-Nyquist computational ghost imaging with deep learning,” Opt. Express 28(3), 3846 (2020). [CrossRef]  

19. X. Xu, E. Li, X. Shen, and S. Han, “Optimization of speckle patterns in ghost imaging via sparse constraints by mutual coherence minimization,” Chin. Opt. Lett. 13, 57–60 (2015). [CrossRef]  

20. L. Wang and S. Zhao, “Fast reconstructed and high-quality ghost imaging with fast Walsh-Hadamard transform,” Photonics Res. 4(6), 240–244 (2016). [CrossRef]  

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

22. H. Cui, Q. Hao, J. Cao, Z. Wang, H. Zhang, and Y. Cheng, “Curved retina-like camera array imaging system with adjustable super-resolution fovea,” Appl. Opt. 60(6), 1535 (2021). [CrossRef]  

23. L. Peng, L. Bian, T. Liu, and J. Zhang, “Agile wide-field imaging with selective high resolution,” Opt. Express 29(22), 35602–35612 (2021). [CrossRef]  

24. J. Hua, E. Hua, F. Zhou, J. Shi, C. Wang, H. Duan, Y. Hu, W. Qiao, and L. Chen, “Foveated glasses-free 3D display with ultrawide field of view via a large-scale 2D-metagrating complex,” Light: Sci. Appl. 10(1), 213 (2021). [CrossRef]  

25. Q. Hao, Y. Tao, J. Cao, M. Tang, Y. Cheng, D. Zhou, Y. Ning, C. Bao, and H. Cui, “Retina-like Imaging and Its Applications: A Brief Review,” Appl. Sci. 11(15), 7058 (2021). [CrossRef]  

26. M. G. Bickell, L. Zhou, and J. Nuyts, “Spatially Variant Resolution Modelling for Iterative List-Mode PET Reconstruction,” IEEE Trans. Med. Imaging 35(7), 1707–1718 (2016). [CrossRef]  

27. E. Akbas and M. P. Eckstein, “Object detection through search with a foveated visual system,” PLoS Comput. Biol. 13(10), e1005743 (2017). [CrossRef]  

28. H. Kim, J. Yang, M. Choi, J. Lee, S. Yoon, Y. Kim, and W. Park, “Eye tracking based foveated rendering for 360 VR tiled video,” in Proceedings of the 9th ACM Multimedia Systems Conference (2018), pp. 484–486.

29. F. Huang, H. Ren, X. Wu, and P. Wang, “Flexible foveated imaging using a single Risley-prism imaging system,” Opt. Express 29(24), 40072 (2021). [CrossRef]  

30. D. B. Phillips, M. Sun, J. M. Taylor, M. P. Edgar, S. M. Barnett, G. M. Gibson, and M. J. Padgett, “Adaptive foveated single-pixel imaging with dynamic supersampling,” Sci. Adv. 3(4), e1601782 (2017). [CrossRef]  

31. C. Zhou, T. Tian, C. Gao, W. Gong, and L. Song, “Multi-resolution progressive computational ghost imaging,” J. Opt. (Bristol, U. K.) 21(5), 055702 (2019). [CrossRef]  

32. K. Zhang, J. Cao, Q. Hao, F. Zhang, Y. Feng, and Y. Cheng, “Modeling and Simulations of Retina-Like Three-Dimensional Computational Ghost Imaging,” IEEE Photonics J. 11(1), 1–13 (2019). [CrossRef]  

33. J. Cao, D. Zhou, Y. Zhang, H. Cui, F. Zhang, K. Zhang, and Q. Hao, “Optimization of retina-like illumination patterns in ghost imaging,” Opt. Express 29(22), 36813 (2021). [CrossRef]  

34. X. Zhai, Z. Cheng, Y. Hu, Y. Chen, Z. Liang, and Y. Wei, “Foveated ghost imaging based on deep learning,” Opt. Commun. 448, 69–75 (2019). [CrossRef]  

35. Z. Shin, T. Chai, C. H. Pua, X. Wang, and S. Y. Chua, “Efficient Spatially-Variant Single-Pixel Imaging Using Block-Based Compressed Sensing,” J. Sign. Process Syst. 93(11), 1323–1337 (2021). [CrossRef]  

36. Z. Shin, H. S. Lin, T. Chai, X. Wang, and S. Y. Chua, “Programmable spatially variant single-pixel imaging based on compressive sensing,” J. Electron. Imag. 30(02), 21004 (2021). [CrossRef]  

37. G. Ma, S. Zhao, X. Wang, Y. Qian, J. Zhang, P. Zhang, H. Gao, R. Liu, and F. Li, “Foveated imaging through scattering medium with LG-basis transmission matrix,” Opt. Laser Eng. 159, 107199 (2022). [CrossRef]  

38. D. Li, R. Du, A. Babu, C. D. Brumar, and A. Varshney, “A Log-Rectilinear Transformation for Foveated 360-degree Video Streaming,” IEEE Trans. Visual. Comput. Graphics 27(5), 2638–2647 (2021). [CrossRef]  

39. J. Ye, A. Xie, S. Jabbireddy, Y. Li, X. Yang, and X. Meng, “Rectangular Mapping-based Foveated Rendering,” in IEEE Conference on Virtual Reality and 3D User Interfaces (VR) (2022), pp. 756–764.

40. Q. Hao, M. Tang, J. Cao, S. Rizvi, H. Cui, and Z. Wu, “Study on Rotation and Scaling Invariance of Retina-Like Imaging,” IEEE Photonics J. 12(5), 1–10 (2020). [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 (7)

Fig. 1.
Fig. 1. Schematic of the FPGI system. Given the designed variant-resolution annular pattern sequence via log-rectilinear transformation, the light modulated by the DMD is projected on the curved mirror through a projecting lens. After the multi-reflection of light between the curved mirror and the omnidirectional scene, the backscattered light intensity is collected by two single-pixel detectors to collect uniformly more and more light intensity from all directions. And then the unwrapping-free foveated panoramic ghost image can be reconstructed by the correlation calculation between the collected light intensity and the variant-resolution rectangular patterns via the log-polar mapping.
Fig. 2.
Fig. 2. Schematic of the variant-resolution annular pattern structure via log-rectilinear transformation. (a) Uniform-resolution rectangular pattern structure in log-polar coordinate with P × Q pixels. (b) Variant-resolution allocation structure with U × V pixels. (c) Detail of the one-to-one variant-resolution allocation structure in the p(u)-q(v) coordinate system. (d) Uniform-resolution annular pattern structure via log-polar mapping. (e) Variant-resolution annular pattern structure via log-polar mapping. (f) Detail of part of the variant-resolution annular structure.
Fig. 3.
Fig. 3. Schematic of variant-resolution annular pattern structure with real fovea. (a) The variant-resolution rectangular structure with traditional fovea in the log-polar coordinate, where the ROI near 2π is at the edge and the area in green is the farthest from the ROI. (b) The variant-resolution annular structure with traditional fovea in the log-polar coordinate, where the area in green is close to the ROI. (c) Variant-resolution structure with the real fovea via the horizontal shifting of the abscissa of the log-polar coordinate system to center the ROI in the whole FOV. (d) Variant-resolution annular structure with the real fovea with a rotated Cartesian coordinate system. (e) After optimization with the real fovea, the variant-resolution allocation structure in the original log-polar coordinate, where the area in green has a higher-resolution allocation than the one in (a). (f) Variant-resolution annular pattern structures with the real fovea in the original Cartesian coordinate system. Note that ξ'-θ’ is the shifted coordinate to center the ROI in the log-polar coordinate system from φ0 to φ0 + 2π and X'-Y’ is the rotated coordinate axis in the Cartesian coordinate system with the rotated angle φ0.
Fig. 4.
Fig. 4. FPGI with variant resolution. (a) A sample object to make an omnidirectional 360° FOV, where the ROI is highlighted by the red box. (b) A uniform-resolution annular pattern sample randomly filled by 0 or 1 with 20 × 102 pixels and the corresponding uniform-resolution rectangular pattern via log-polar mapping. (c) With the parameters of ‘Variant-resolution 1’ in Table 1, a random variant-resolution annular pattern sample with 19 × 66 pixels and the log-polar mapped variant-resolution rectangular pattern. (d) With the parameters of ‘Variant-resolution 2’ in Table 1, a random variant-resolution annular pattern sample with 19 × 50 pixels and the log-polar mapped variant-resolution rectangular pattern. (e) With the parameters of ‘Variant-resolution 3’ in Table 1, a random variant-resolution annular pattern sample with 16 × 38 pixels and the log-polar mapped variant-resolution rectangular pattern. (f) With the same measurements of full sampling, the unwrapping-free reconstructed panoramic ghost images using the four annular pattern sequences with different structures of (b)-(e) respectively, and the corresponding magnified views of the ROIs. (g) Quantitative comparison on the PSNR of fovea of the reconstructed image using a uniform-resolution annular pattern structure and the three variant-resolution pattern structures with various parameters shown in (b)-(e) at different sampling rates. (h) Quantitative comparison on the reconstruction time of unwrapping-free PGI using the four different pattern structures at different sampling rates. Note that, the sampling ratio refers to the uniform-resolution panoramic ghost image with 20 × 102 pixels, because the variant-resolution panoramic ghost image is also reformatted into the uniform-resolution grid to ensure high-resolution imaging on the ROI.
Fig. 5.
Fig. 5. FPGI with variant resolution. (a) A complex object to make an omnidirectional 360° FOV, where the ROI is highlighted by the red box. (b) A uniform-resolution annular pattern sample randomly filled by 0 or 1 with 40 × 204 pixels and the log-polar mapped uniform-resolution rectangular pattern. (c) With the parameters of ‘Variant-resolution 1’ in Table 2, a random variant-resolution annular pattern sample with 29 × 98 pixels and the log-polar mapped variant-resolution rectangular pattern. (d) With the parameters of ‘Variant-resolution 2’ in Table 2, a random variant-resolution annular pattern sample with 20 × 98 pixels and the log-polar mapped variant-resolution rectangular pattern. (e) With the parameters of ‘Variant-resolution 3’ in Table 2, a random variant-resolution annular pattern sample with 20 × 50 pixels and the log-polar mapped variant-resolution rectangular pattern. (f) With the same measurements of full sampling, the unwrapping-free reconstructed panoramic ghost images with the four annular pattern sequences with different structures of (b)-(e) respectively, and the corresponding magnified views of the ROIs. (g) Quantitative comparison on the SSIM of fovea of the reconstructed image using a uniform-resolution annular pattern structure and the three variant-resolution pattern structures with different parameters shown in (b)-(e) at different sampling rates. (h) Quantitative comparison on the reconstruction time of PGI with the four different pattern structures at different sampling rates. Note that, the sampling ratio refers to the uniform-resolution panoramic ghost image with 40 × 204 pixels.
Fig. 6.
Fig. 6. FPGI with a real fovea. (a) The object with an ROI located in the position close to the original start boundary 0. (b) A uniform-resolution annular pattern sample randomly filled by 0 or 1 with 40 × 204 pixels. (c) A traditional-fovea variant-resolution annular pattern sample with 39 × 86 pixels, where the resolution of the areas circled in green is lost. (d) A real-fovea variant-resolution pattern sample with 39 × 101 pixels, which avoids the resolution loss of the areas circles in green. (e) The unwrapping-free uniform-resolution panoramic ghost image with full sampling in the original coordinate and the corresponding image where the ROI is at the center. (f) The traditional-fovea variant-resolution panoramic ghost image with full sampling in the original coordinate and the corresponding image with centering the ROI, where the areas marked by the green box are not as clear as the areas marked by the blue box (Quantitatively, the SSIM of the green box is 0.35 while the SSIM of the blue box is 0.61). (g) The real-fovea variant-resolution panoramic ghost image with full sampling and the corresponding image with centering the ROI, where the areas marked by the green box are as clear as the areas marked by the blue box (Quantitatively, the SSIM of the green box is improved to 0.60, which is close to the SSIM of the blue box of 0.61).
Fig. 7.
Fig. 7. FPGI with multiple foveae. (a) The object with the two ROIs marked in red and green respectively. (b) A random variant-resolution annular pattern sample with double foveae. (c) Reconstructed uniform-resolution panoramic ghost image with full sampling and the magnified views of the two ROIs following the quantitative evaluation of the PSNR and SSIM. (d) Reconstructed variant-resolution panoramic ghost image using the patterns shown in (b) with full sampling and the magnified views of the two ROIs following the quantitative evaluation of the PSNR and SSIM. (e) The object with the three ROIs marked in red, green and blue respectively. (f) A random variant-resolution annular pattern sample with treble foveae. (g) Reconstructed uniform-resolution panoramic ghost image with full sampling and magnified views of the three ROIs following the quantitative evaluation of the PSNR and SSIM. (h) The reconstructed variant-resolution panoramic ghost image using the patterns shown in (f) with full sampling and magnified views of the three ROIs following the quantitative evaluation of the PSNR and SSIM. (i) Comparison on the reconstruction time between the traditional uniform-resolution PGI and the FPGI with the two ROIs at different sampling rates. (j) Comparison on the reconstruction time between the traditional uniform-resolution PGI and the FPGI with the three ROIs at different sampling rates. Note that, the sampling ratio refers to the uniform-resolution panoramic ghost image with 40 × 204 pixels.

Tables (2)

Tables Icon

Table 1. Parameters of the four pattern structures in Fig. 4

Tables Icon

Table 2. Parameters of the four pattern structures in Fig. 5

Equations (7)

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

{ R 1 u = { u , w h e n 1 u m _ 0 , m _ 0 α 1 u m _ 0 , w h e n m _ 0 u ( U 1 ) / 2 , R 2 v = { v , w h e n 1 v n _ 0 , n _ 0 α 2 v n _ 0 , w h e n n _ 0 v ( V 1 ) / 2 , p ( u ) = { max ( 1 , r o u n d ( p 0 R 1 u ) ) , w h e n 1 u ( U 1 ) / 2 , p 0 , w h e n u = ( U + 1 ) / 2 , min ( r o u n d ( p 0 + R 1 u ) , P ) , w h e n ( U + 1 ) / 2 < u U , q ( v ) = { max ( 1 , r o u n d ( q 0 R 2 v ) ) , w h e n 1 v ( V 1 ) / 2 , q 0 , w h e n v = ( V + 1 ) / 2 , min ( r o u n d ( q 0 + R 2 v ) , Q ) , w h e n ( V + 1 ) / 2 < v V , U 2 N + 1 , V 2 N + 1
{ ε = 1  +  sin ( π / Q ) 1 sin ( π / Q ) r c 1 = r 0 1 sin ( π / Q ) r p ( u ) = r 0 ε p ( u ) r c p ( u ) = r c 1 ε p ( u )  - 1 θ q ( v ) = q ( v ) 2 π Q   ( v = 1 , 2 , 3 V ) ξ p ( u ) = log ε ( r c p ( u ) ) = log ε ( r c 1 ) + p ( u ) 1   ( u = 1 , 2 , 3 U ) ,
{ φ 0 = ( q c q 0 ) 2 π Q θ = θ + φ 0 ,
I = x y P v a n n u l a r ( x , y ) O o m n i d i r e c t i o n a l ( x , y ) ,
P v r e c t a n g u l a r = A p P n A q T ,
min   | | c | | l 1 s . t .   H O v p a n o r a m i c = c ,   P v r e c t a n g u l a r O v p a n o r a m i c = I
{ PSNR = 10 log 10 ( 2 k 1 ) 2 MSE MSE = 1 S i , j ( F ( i , j ) F ( i , j ) ) 2 SSIM M i , j = ( 2 μ μ + c 1 ) ( 2 w + c 2 ) ( μ 2 + u 2 + c 1 ) ( σ 2 + σ 2 + c 2 ) ,
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.