Open AccessArticle

Compressive Sensing Based Three-Dimensional Imaging Method with Electro-Optic Modulation for Nonscanning Laser Radar

School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China

Science and Technology on Space Microwave Laboratory, China Academy of Space Technology, Beijing 100081, China

Author to whom correspondence should be addressed.

Symmetry 2020, 12(5), 748; https://doi.org/10.3390/sym12050748

Submission received: 5 March 2020 / Revised: 31 March 2020 / Accepted: 7 April 2020 / Published: 6 May 2020

Browse Figures

Versions Notes

Abstract

Low-cost Laser Detection and Ranging (LiDAR) is crucial to three-dimensional (3D) imaging in applications such as remote sensing, target detection, and machine vision. In conventional nonscanning time-of-flight (TOF) LiDAR, the intensity map is obtained by a detector array and the depth map is measured in the time domain which requires costly sensors and short laser pulses. To overcome such limitations, this paper presents a nonscanning 3D laser imaging method that combines compressive sensing (CS) techniques and electro-optic modulation. In this novel scheme, electro-optic modulation is applied to map the range information into the intensity of echo pulses symmetrically and the measurements of pattern projection with symmetrical structure are received by the low bandwidth detector. The 3D imaging can be extracted from two gain modulated images that are recovered by solving underdetermined inverse problems. An integrated regularization model is proposed for the recovery problems and the minimization functional model is solved by a proposed algorithm applying the alternating direction method of multiplier (ADMM) technique. The simulation results on various subrates for 3D imaging indicate that our proposed method is feasible and achieves performance improvement over conventional methods in systems with hardware limitations. This novel method will be highly valuable for practical applications with advantages of low cost and flexible structure at wavelengths beyond visible spectrum.

Keywords:

compressive sensing; image reconstruction; infrared imaging; laser radar; reconstruction algorithm

1. Introduction

Laser Detection and Ranging (LiDAR) is a remote sensing technology that obtains intensity and range of objects by transmitting laser pulses and receiving return pulses [1]. As it displays advantages of high-reliability, longitudinal resolution, and anti-electric magnetic interference ability, LiDAR is widely used in multiple areas such as target detection and recognition, remote sensing, and three-dimensional (3D) scene imaging [2,3], wherein nonscanning LiDAR system is undergoing major development in recent years because of features such as low requirement of light source, ability to detect moving objects, and image obscuring objects [4,5,6].

At present, an important experimental-research branch of 3D LiDAR system is achieving single-pixel nonscanning imaging and most of this research are based on the new sensing modality called compressive sensing (CS) [7,8,9]. This new mathematical theory was proposed by Donoho and Candès [10,11] initially and has expanded greatly since then. The CS framework for the signal acquisition was a major breakthrough in signal processing community after the famous Shannon sampling theorem. It tried to sense enough information of a sparse signal by projecting the data in a space incoherent with the data’s structure, so that the original signal can be reconstructed perfectly [12,13]. Thanks to the effort of Duarte and Davenport [14] who presented a single-pixel CS camera architecture, CS theory has expanded greatly from methodology to practice in both active and passive imaging technologies. Duarte and Davenport’s work paved the way for various succeeding 3D LiDAR applications such as Howland experimentally demonstrating a photon-counting single-pixel time-of-flight (CS-TOF) laser radar camera [15] and Sun showing a modified TOF 3D imaging system (M-CS-TOF) with an accuracy of ∼3 mm [16].

TOF imaging is performed by correlating the detection time of the back-scattered light with the time of illumination pulses [17]. As a direct ranging method, TOF imaging system depends on narrow bandwidth of illumination laser pulses and the short rise time of detectors to achieve a high range resolution. To overcome such limitations, indirect methods including time slicing (TS), intensity correlation and gain modulation in range-gated imaging systems have become increasingly popular [18,19,20]. The TS technique obtains a full depth map by recording multiple frames whose number is proportional to the required depth resolution [21]. Therefore, tons of 2D images are required to construct one 3D image with high resolution. As to intensity correlation, the range information can be restored by analyzing the intensity variation in the overlapping areas of two range-gates with a small number of distinct images [19]. However, rectangular-shaped or triangular-shaped pulses with high quality are required [22]. Furthermore, the gate widths should be carefully designed in order to adapt to the laser pulse widths [23]. On the contrary, gain modulation method is a pulse-shape-free method. It only needs two gain-modulated images to construct time-resolved imaging, making it possible to achieve accurate imaging reconstruction with less data and low-cost system [24].

In this paper, a CS-based electro-optic modulation scheme in nonscanning LiDAR and a 3D reconstruction method are proposed. This scheme combines the CS technology with gain modulation to achieve the 3D imaging of the target. Laser pulses are emitted toward the target and then the scattered light from the target returns to the system. Peak values are recorded out of total echo pulses which is the sum of echo pulses after spatial modulation and gain modulation are applied in a symmetrical procedure. Next, two modulated images can be recovered by the peak values in the CS framework corresponding to different gain modulation functions respectively. Furthermore, an intensity map and a depth map of the target scene are able to be extracted from the two recovered images which are obtained by solving underdetermined inverse problems. We propose an integrated regularization model for the recovery problems and this minimization functional model is solved by a proposed algorithm applying the alternating direction method of multiplier (ADMM) technique. As an extension of our previous related work [25], the novel resolution to the 3D imaging problem lies in the proposed method that can overcome hardware limitations, by modulating the range information into the reflective light intensity before compressive sampling. Numerical simulations for 3D imaging with various subrates are presented to demonstrate the improved performance of the proposed method.

This paper is organized as follows. Section 2 presents a prototype of the CS-based RGI system applied the electro-optic modulation method for 3D imaging and its mathematical model. Section 3 presents the 3D measurement and reconstruction using CS-based electro-optic modulation method after we review the CS theory. Section 4 shows simulation experiments based on typical discrete and continuous targets to validate the proposed method. Finally, Section 5 presents the conclusion.

2. System Description

In this section, we propose a 3D LiDAR system to implement the CS-based electro-optic modulation method. In addition, the working mechanism of the system is presented by introducing the proposed system setup and its laser pulse model.

2.1. Proposed System Setup

As is shown in Figure 1, the proposed CS range-gated imaging (RGI) system mainly consists of a transmitter, a receiver, a narrowband filter (NBF), an electro-optic modulator (EOM), two polarizers, a digital micromirror device (DMD), and a one-pixel avalanche photodiode (APD) detector. There is a red light path with white arrows and blue data path in dotted line in Figure 1.

First, a pulsed laser is emitted toward the target scene. Second, reflected light returns to the RGI system. With closing the receiver’s gate while the pulsed laser is traveling, the system’s Signal-to-Noise Ratio (SNR) can be improved. Third, an EOM is used to conduct gain modulation and a DMD is applied to conduct spatial modulation with the gate open [26,27,28]. Forth, a NBF is used to remove background radiation from sunlight and other sources. Finally, the reflected light arrives at a one-pixel APD detector.

Figure 2 contains a sequence chart for the RGI system and gain modulation function. In Figure 2a,

V_{i n} (t)

and

V_{d e} (t)

are the applied voltages of the EOM to generate different gain modulations. The index “in” is short for increasing and the index “de” is short for decreasing. The intensities of the passing light are modulated with

g_{i n} (t)

and

g_{d e} (t)

, respectively. Here, we define

R_{0}

is the range according to the optical gate open beginning time; D with subscript is the range between the object and

R_{0}

in Figure 2b,c. R is the range between the RGI system and the object; and L is the gate opening range. The laser pulse is emitted at the time of

T_{0}

. The gate opens at the time of

T_{1}

and the opening time of RGI system is determined by the duration time of gate opening,

T_{G}

. According to the round-trip travel range, L is given by

L = \frac{1}{2} c T_{G},

(1)

where c is the speed of light. Then, we can get

R = R_{0} + D .

(2)

The pulse received by the APD is a sum of overlapped pulses with particular modulation from different distances. Figure 2b,c explain how the components of the pulse of APD are modulated and work together. The pulses in black dash line are echo pulses from objects at different distances and they overlap with each other that results in the pulse of APD,

p_{t} (t)

, in black full lines. The red dash lines are gain modulation functions

g_{i n} (t)

and

g_{d e} (t)

when the gate opens, respectively. Then, the range D and R are able to be resolved from two intensity-modulated CS-based recovered images which will be introduced in detail later. Ultimately, time-resolved imaging can be achieved.

2.2. Proposed System Model

To be more specific, the system model is presented as follows. The time propagation of the laser pulse from the transmitter

p_{t}

is of the form:

p_{t} (t) = \{\begin{matrix} {(t / τ)}^{2} e^{(- \frac{t}{τ})} & t \geq 0 \\ 0 & t < 0 \end{matrix}

(3)

where

τ = T_{1 / 2} / 3.5

and

T_{1 / 2}

is the full width at half maximum (FWHM) of the pulse [3]. Assuming that the target is a Lambert diffuse point-like reflector, the energy of the echo signal

E_{r}

can be represented as the following formula:

E_{r} = \frac{ρ D_{r}^{2} T_{a}^{2} η}{4 R^{4}} \int_{T_{G}} p_{t} (t - 2 R / c) d t,

(4)

where

ρ

is the reflectance of the target,

D_{r}

is the diameter of the receiver,

T_{a}

is the single-pass atmospheric transmittance,

η

is the efficiencies of the optical transmitting and receiving systems and R is the range between the target and the receiver [29]. Since the APD can only obtain the intensity of the pulses rather than their flying time, gain modulation by an EOM based on the electro-optic effect of crystals is applied to realize time-resolved imaging.

Figure 3 presents a typical structure of longitudinal electro-optic modulation using Pockels effect. It mainly consists of the electro-optic crystals and the polarizers,

P_{1}

and

P_{2}

, in a universal coordinate system with x-axis, y-axis and z-axis. The polarizers are fixed and perpendicular to each other. The polarized light whose vibration is parallel to x axis enters the EOM and it is resolved into two components parallel to

x^{'}

axis and

y^{'}

axis respectively. When a voltage is applied to the EOM in the direction of propagation (z axis), phase retardation of the light will take place between the two components and the intensity of the light is modulated. The phase shift

θ

proportional to the applied voltage

V (t)

can be given by

θ = \frac{2 π}{λ} n_{0}^{3} γ_{63} V (t),

(5)

where

λ

is the wavelength of light,

n_{0}

is the ordinary refractive indices of crystal and

γ_{63}

is the crystal electro-optic coefficient [24]. The applied voltages are given in the monotonously increasing and decreasing form of:

\begin{matrix} V_{i n} (t) = \frac{V_{π}}{T_{G}} t, (0 \leq t \leq T_{G}) \\ V_{d e} (t) = V_{π} - \frac{V_{π}}{T_{G}} t, (0 \leq t \leq T_{G}) \end{matrix},

(6)

where the half-wave voltage

V_{π} (θ = π)

is given by

V_{π} = \frac{λ}{2 n_{0}^{3} γ_{63}} .

(7)

Then we can have the phase shift

θ

as a function of time t:

θ = π \cdot \frac{t}{T_{G}}, (0 \leq t \leq T_{G}) .

(8)

Applying Equation (1) to Equation (8), the function of the phase shift

θ

and the depth D are obtained as:

θ = π \cdot \frac{D}{L}, (0 \leq D \leq L) .

(9)

During the gated opening time, the applied voltages are

V_{i n} (t)

and

V_{d e} (t)

. The intensities of the light,

E_{i n}

and

E_{d e}

, are modulated with gain functions

g_{i n} (t)

and

g_{d e} (t)

respectively:

\begin{matrix} E_{i n} = E_{r} {sin}^{2} (θ / 2) \\ E_{d e} = E_{r} {cos}^{2} (θ / 2) \end{matrix} .

(10)

3. CS-Based Electro-Optic Modulation Method for 3D Imaging

3.1. Compressive Sensing Theory

Compressive sensing is different from the traditional approach to digital data acquisition that samples an analog signal uniformly at or above the Nyquist rate. CS is a measurement technique to employ optimization to detect a sparse n-dimensional signal with

m < n

samples [30].

The detection and reconstruction process can be modeled by vectors and matrices as follows. For convenience sake, we use n-dimensional vector x represents the gray value of an original image. In addition, we then take advantage of transform coding to represent the image in terms of the coefficients

\{α_{i}\}

of an orthonormal basis expansion:

x = \sum_{i = 1}^{n} α_{i} ψ_{i},

(11)

where

{\{ψ_{i}\}}_{i = 1}^{n}

are the

n \times 1

basis vectors. We can concisely write the original image as

x = Ψ α,

(12)

by forming the coefficient vector

α

and the

n \times n

transforming matrix

Ψ

which is the stack of vectors

{\{ψ_{i}\}}_{i = 1}^{n}

as columns. For natural images, x is mostly compactly represented by considering only significant elements in the sparse basis, such as discrete cosine transform (DCT) and wavelet transform (DWT) [31,32].

In addition, we then construct a

m ≪ n

dimensional measurement vector y by multiplying x by an

m \times n

matrix

Φ

describing the measurement matrix:

y = Φ x + ε = A α + ε,

(13)

where vector

ε

is the additional noise. As

m ≪ n

, the dimension of the signal has been reduced. Matrix A is called sensing matrix and determined by

A = Φ Ψ

. The subsampling rate is defined by

s u b r a t e = m / n

and conventional recover methods are incapable since

s u b r a t e < 1

Previous research has shown that the sensing matrix obeys the so-called restricted isometry property (RIP) when

Φ

is drawn randomly from a suitable distribution [33]. In this way, the inverse problem of constructing the signal x from the measurements y can be achieved by solving a 1-norm (

p = 1

) or 2-norm (

p = 2

) minimization problem with relaxed constraints:

min_{α \in R^{n}} {| | α | |}_{p}, s . t . | | A α - {y | |}_{2} \leq ε,

(14)

where

ε

bounds the amount of noise in the sampling [34]. Equation (14) is a convex problem and can be solved efficiently by various construction algorithms such as the Total Variation (TV) minimization based on augmented Lagrangian and alternating direction algorithm (TVAL3) which are able to handle different boundary conditions for

α

and widely used in CS systems [35].

3.2. Peak Value Obtained for 3D Reconstruction

Assuming the power of the echo signal is

P_{r}

, the range-discretized form of the total power of echo pulses by sparse sampling in the view of field can be described as a sum of the pulses from each pixel with different amplitudes and delay times:

p_{r, j} (t) = \sum_{i = 1}^{n} Φ_{(j, i)} g (t) P_{r} (R_{i}) p (t - t_{i}),

(15)

where

Φ

is a

m \times n

measurement matrix,

t_{i} = 2 R_{i} / c

and

g (t)

is the gain modulation function applied by the EOM.

For 3D compressive imaging, a 2D extended mathematical model of measuring process is presented as follows according to Equations (13) and (15) in a discretized matrix form:

Y = Φ X F_{g a i n} + ε_{M},

(16)

where Y is a

m \times l

matrix as total echo waveform, X is a

n \times l

matrix as echo waveform from n pixels of the target scene,

F_{g a i n}

is a

l \times l

diagonal matrix as the gain modulation function and

ε_{M}

is the additive noise matrix. Only peak values of the total echo waveform are recorded as the compressive measurements rather than recording the whole waveform. In this instance, we define vector

y = {(y_{1}, y_{2}, \dots y_{m})}^{T}

as the CS measurements in the form of:

y = Φ x_{g} + ε = Φ G x + ε,

(17)

where

x_{g} = {(x_{g, 1}, x_{g, 2}, \dots, x_{g, n})}^{T}

is the

n \times 1

vector used to present the modulated intensity map and

ε

is the additive noise vector.

x_{g}

can be presented as the product of the gain value matrix

G = d i a g (g_{1}, g_{2}, \dots, g_{n})

where

d i a g (\cdot)

represents a diagonal matrix and the

n \times 1

vector x which presents the original echo energy from the target. Applying different voltages,

V_{i n} (t)

and

V_{d e} (t)

, to the EOM, the original intensity map are modulated to

x_{i n}

and

x_{d e}

. Then two sets of measurements

y_{i n}

and

y_{d e}

are achieved in the CS framework shown in Equation (17), respectively.

3.3. 3D Image Reconstruction

In this section, we proposed a recovery method to achieve the goal of 3D imaging according to the CS measurements

y_{i n}

and

y_{d e}

. A set of inverse problem formulations within the CS framework are presented and the alternating direction method of multipliers (ADMM) combined with convex optimization algorithms is employed. In this way, the modulated intensity images are obtained to determine a range value and a gray value for each pixel in the target scene, namely

\hat{x} = {(\hat{x_{1}}, \hat{x_{2}}, \dots, \hat{x_{n}})}^{T}

and

\hat{r} = {(\hat{r_{1}}, \hat{r_{2}}, \dots, \hat{r_{n}})}^{T}

, respectively.

For clarity, we present the electro-optical modulated CS measurement vectors in the form of Equation (17):

\begin{matrix} y_{i n} = Φ x_{i n} + ε_{i n} = Φ G_{i n} x + ε_{i n} \end{matrix}

(18)

\begin{matrix} y_{d e} = Φ x_{d e} + ε_{d e} = Φ G_{d e} x + ε_{i n} . \end{matrix}

(19)

There is a useful property of the gain functions for one particular target according to the proposed system model in Section 2.2, especially Equation (10). The property shows that the sum of the two gain value matrices is not only certain but also known:

G_{i n} + G_{d e} = I,

(20)

where I is the unit matrix.

We enforce the sparsity of the signal to be constructed in some transform domain as a prior and construct the inverse problem generally. However, reconstructions based on primary regularizations may lead to the serious distortion for the modulated images with complex texture which contain the intensity and depth information at the same time. Therefore the shearlet transform and a generalization of TV (TGV) are adopted as the regularization [36]. Furthermore, the optimal solutions of the modulated image and the intensity map are supposed to be achieved simultaneously in the consideration of the requirement of subsequent computation. In addition to the regularization of the modulated image

x_{i n}

x_{d e}

, the intensity map x is regularized by TV norm as a natural signal. Thus, we formulate the proposed 3D reconstruction problem with the regularization terms of both

x_{i n}

and x as follows:

min_{x_{i n}, x} λ_{S H} \sum_{j = 1}^{N} | | S H_{j} (x_{i n}) | | + T G V_{α}^{2} (x_{i n}) + \sum_{i} | | D_{i} x | |_{2}

(21)

where

S H_{j} (\cdot)

is the

j th

subband of the shearlet transform, N is the total number of subbands which is related to the number of scales,

λ_{S H} > 0

is a balancing factor relying on the gradients and the sparsity of the images under the shearlet transform [36] and

D_{i} x \in R^{2}

is the discrete gradient of x at pixel i, namely the TV norm.

It should be pointed out that the optimization problem in Equation (21) is constrained by Equations (18) and (19). Incorporating Equation (20) into them leads to the constrained term with respect to

x_{i n}

and x and the integrated optimization problem with observation constraint is as follows:

min_{x_{i n}, x} λ_{S H} \sum_{j = 1}^{N} | | S H_{j} (x_{i n}) | | + T G V_{α}^{2} (x_{i n}) + \sum_{i} | | D_{i} x | |_{2} s . t . 2 y_{i n} + y_{d e} = Φ x_{i n} + Φ x .

(22)

There is no off-the-shelf algorithm to solve this problem, hence we use the alternating direction method of multipliers (ADMM) to split the problem into multiple subproblems that can be solved easily [37].

It is straightforward to convert the equality constrained minimization problem into an unconstrained regulated Lagrangian optimization problem using the augmented Lagrangian method. For Equation (22), its augmented Lagrange function is

\begin{matrix} L_{A} (x_{i n}, x, λ) = & λ_{S H} \sum_{j = 1}^{N} | | S H_{j} (x_{i n}) | | + T G V_{α}^{2} (x_{i n}) + \sum_{i} | | D_{i} x | |_{2} \\ + λ^{T} (Φ x_{i n} + Φ x - 2 y_{i n} + y_{d e}) + \frac{ρ}{2} | | Φ x_{i n} + Φ x - 2 y_{i n} + y_{d e} {| |}_{2}^{2}, \end{matrix}

(23)

where

λ

is the Lagrangian multiplier and

ρ

is the positive regularization parameter associated with the penalty term. By doing this, the solution of Equation (21) is obtained by seeking a saddle point of

L_{A} (x_{i n}, x)

To solve the problem of minimizing the augmented Lagrangian function efficiently, the ADMM technique is embedded here to decompose it into two subproblems as follows. The

x_{i n} - subproblem

is separable with respect to

x_{i n}

. Given

x^{k}

, the optimization problem associated with

x_{i n}

reduces to

\begin{matrix} min_{x_{i n}} L_{A} (x_{i n}, x^{k}) = & λ_{S H} \sum_{j = 1}^{N} | | S H_{j} (x_{i n}) | | + T G V_{α}^{2} (x_{i n}) + \sum_{i} | | D_{i} x^{k} | |_{2} \\ + λ^{T} (Φ x_{i n} + Φ x^{k} - 2 y_{i n} + y_{d e}) + \frac{ρ}{2} | | Φ x_{i n} + Φ x^{k} - 2 y_{i n} + y_{d e} {| |}_{2}^{2}, \end{matrix}

(24)

where k is the iteration number. The

x_{i n} - subproblem

can be reformulated as a model of several nondifferentiable

l_{1}

terms by deriving another form of

T G V_{α}^{2}

in terms of

l_{1}

minimization. By introducing auxiliary variables and quadratic penalty terms, the subproblem is able to be solved using the detail-preserving regularization scheme-based algorithm in [36].

With

x_{i n}

fixed, the

x - subproblem

can be written as

\begin{matrix} min_{x} L_{A} (x_{i n}^{k}, x) = & λ_{S H} \sum_{j = 1}^{N} | | S H_{j} (x_{i n}^{k}) | | + T G V_{α}^{2} (x_{i n}^{k}) + \sum_{i} | | D_{i} x | |_{2} \\ + λ^{T} (Φ x_{i n}^{k} + Φ x - 2 y_{i n} + y_{d e}) + \frac{ρ}{2} | | Φ x_{i n}^{k} + Φ x - 2 y_{i n} + y_{d e} {| |}_{2}^{2}, \end{matrix}

(25)

which is equivalent to a quadratic function. It is difficult to solve this function directly because of the non-differentiability and non-linearity of the regularization terms. The alternating minimization scheme and multipliers updating derived algorithm may be used to tackle the subproblem efficiently and accurately which is called TVAL3 presented in [38]. In special, it has been proven to accelerate the TVAL3 scheme that the structured measurement matrices derived from Hadamard transform are adopted in the CS framework.

According to the ADMM scheme, the multiplier should be updated at each iteration with the update formula of

λ

as follows:

λ^{k + 1} = λ^{k} + ρ (Φ x_{i n}^{k + 1} + Φ x^{k + 1} - 2 y_{i n} - y_{d e}) .

(26)

All derivations above complete the discussion of the optimization to minimize the augmented Lagrangian function (23) and the reconstruction of the modulated intensity map corresponding to

x_{i n}

{\hat{x}}_{i n}

and the intensity map

\hat{x}

is achieved from the CS measurements. The other modulated intensity map, denoted as

{\hat{x}}_{d e}

, can be achieved in the same method. The computational process can be accelerated due to the knowledge of

\hat{x}

which has been obtained in the previous stage.

We can estimate the distance value of each pixel

{\hat{r}}_{i}

in the target scene from the two modulated maps according to Equations (9), (10) and (17):

{\hat{r}}_{i} = R_{0} + \frac{2 L}{π} arctan \sqrt{\frac{{\hat{x}}_{i n} (i)}{{\hat{x}}_{d e} (i)}} .

(27)

In total, the constructed 3D image of the target scene can be achieved as a combination of the intensity map

\hat{x}

and the depth map

\hat{r}

. A summary of the procedures of the presented 3D reconstruction method is provided as follows in Algorithm 1.

Algorithm 1: The proposed 3D reconstruction algorithm

4. Simulation Results and Discussion

To validate the proposed method for 3D LiDAR imaging, a simulation system is established and the simulation results are presented in this section. The numerical simulations are done with CPU of Intel Core i3-3220 and 4.00 GB RAM by matlab 2016a (64 bits). The pulsed laser works as an illuminator at low repetition rate, high pulse energy and short pulse duration. A series of 2D Hadamard derived patterns are used to modulate the echo pulse for CS measurements [39]. Meanwhile, the EOM is used to modulate the intensity of echo pulses with gain functions

g_{i n} (t)

and

g_{d e} (t)

, respectively. The main parameters of the 3D imaging LiDAR system are give in Table 1. Generally, the range resolution of TOF compressive LiDAR with similar parameters is about

c / 2

times FWHM in seconds which mainly depends on FWHM of the laser pulse. The proposed method can provide range super-resolution ability if we select it as the standard.

4.1. Reconstruction Performance of the Discrete Target

As shown in Figure 4a, a 3D scene is designed as the space discrete target which consists of two square cardboards cut-outs of the letters “U” and “R”. This is a typical target scene in CS-based 3D imaging [4,15,17,40]. The boards are distinguished and located parallelly at different positions in the scene. The distance between the RGI system and the cardboards ranges from 20.3 m to 20.9 m and therefore the gate opening range of 1 m is designated. The 3D scene is modeled by a combination of two matrices which are showed as the original intensity map, Figure 4b and the original depth map, Figure 4c. The original intensity map shows the reflectivity distribution of the target and the depth map shows the distance between each pixel of the target and the imaging system. In Figure 4c, colors correspond to the distance and the bar on the right side is added as an annotation.

The visual comparison of the imaging results of the “U & R” target scene by various methods is given in Figure 5. The measurements are in the presence of white Gaussian noise with SNR = 30 dB and all of the subrates equal to 10%. Colors correspond to the distance from 20.0 m to 21.0 m and the bar on the right side is added as an annotation. Figure 5a,c,e is the reconstructed intensity maps of the target applying CS-TOF, M-CS-TOF and the proposed method, respectively. These figures show that due to the full use of the CS framework, high-quality intensity maps are obtained which not only describe the shape feature of the scene but also have the edge detail of the target. Figure 5b is the reconstructed depth map from the measurements applying CS-TOF. The letters, “U” and “R”, stand out from the background in the figure but most pixels in them have similar depth values. It means the CS-TOF approach fails to successfully recover the scene in this system. Figure 5d is the reconstructed depth map of multiple planar objects applying M-CS-TOF. Two letters at separate positions are successfully recovered, the depth of each pixel in the target is determined accurately except for some acceptable noise. For a specific subrate, two intensity images containing range information are reconstructed in CS framework by two sets of measurements modulated with

g_{i n} (t)

and

g_{d e} (t)

, respectively. A depth map of the target scene can be obtained from them by the proposed method. Figure 5f is the reconstructed map of the objects at separate positions with unknown depth by applying the proposed method. The object planes and edges are recovered accurately to demonstrate the depth and transverse resolution which indicates the feasibility of the proposed method.

Given subjective observation of the reconstructed maps, Peak Signal to Noise Ratio (PSNR) criterion, Normalized Mean Squared Error (NMSE) metric and a measure of Structural SIMilarity (SSIM) are adopted to assess the reconstruction quality as shown in Figure 6. A standard definition of the PSNR for the reconstructed intensity map

\hat{x}

and the original intensity map x is of the form:

P S N R (\hat{x}, x) (d B) = 10 {log}_{10} \frac{N \cdot max {(x)}^{2}}{\sum_{i = 1}^{N} {(\hat{x} (i) - x (i))}^{2}},

(28)

where N is the total number of pixels in x,

\hat{x} (i)

and

x (i)

are the intensity of the ith pixel corresponding to the reconstructed intensity maps and the original intensity maps, respectively [41]. NMSE is adopted to compare the reconstructed depth maps with the original depth map:

N M S E (\hat{d}, d) = \frac{\sum_{i = 1}^{N} {(d (i) - \hat{d} (i))}^{2}}{\sum_{i = 1}^{N} {(d (i))}^{2}},

(29)

where

d (i)

and

\hat{d} (i)

are the depth of the ith pixel corresponding to the original depth maps and the estimated depth maps, respectively [42].

Besides the criterion above which have clear physical meanings in context of optimization, another kind of criterion is adopted that is well matched to perceived visual quality for extracting structural information. SSIM is widely used by comparing local patterns of pixels after normalized for illuminance and contrast. The index is in a specific form as follows:

SSIM (x, \hat{x}) = \frac{(2 μ_{x} μ_{\hat{x}} + C_{1}) (2 σ_{x \hat{x}} + C_{2})}{(μ_{x}^{2} + μ_{\hat{x}}^{2} + C_{1}) (σ_{x}^{2} + σ_{\hat{x}}^{2} + C_{2})},

(30)

where

μ_{*}

is the mean of the image,

σ_{*}

is the standard deviation,

C_{1}

and

C_{2}

are constant to avoid instability under certain circumstances [43]. The index ranges from 0 to 1 and it equals to 1 when two images are exactly the same.

Figure 6 provide the objective assessment of reconstruction qualities of the “U & R” target scene in the presence of white Gaussian noise with SNR = 30 dB. Figure 6a presents the plots of the PSNR values for reconstructed intensity maps on various subrates (from 2% to 14%) using different approaches. It shows that PSNRs of the proposed method are always higher than those of the TOF-based methods as the subrates increase. The plots corresponding to M-CS-TOF and the proposed method stay stubbornly high and the PSNRs are more than 40 dB when the subrates reach 6% so that high-quality intensity maps are achieved. However, the PSNRs of the CS-TOF reconstructing intensity maps are much lower which results from the error of measurements caused by the echo pulses with different delays overlapping each other. Figure 6b presents the plots of the SSIM values for reconstructed intensity maps on various subrates (from 2% to 14%) using different approaches. The plots corresponding to M-CS-TOF and the proposed method stay stubbornly high and the SSIMs are close to 1 when the subrates reach 6% so that the structures of the target and the reconstructed image are extraordinary similar. On the contrary, the SSIMs of the CS-TOF reconstructing intensity maps are much lower. The SSIMs show similar pattern to the PSNR results, further verify the effectiveness and advantage of the proposed method. Figure 6c presents the plots of the NMSE values for reconstructed depth maps on various subrates (from 2% to 14%) using different approaches. It shows that NMSEs of M-CS-TOF and the proposed method are extremely low especially when the subrates reach 6%. These two methods are feasible and high-quality depth maps are achieved. There are significant errors in the depth maps reconstructed by CS-TOF because it fails to conduct multiple peak detection in overlapped echo pulses. Figure 6d presents the plots of the SSIM values for reconstructed depth maps on various subrates (from 2% to 14%) using different approaches. All of the SSIMs are close to 1 because the structures of the depth distribution are simple and the targets are located at only two unique distances. Still, the performance of the proposed method is as good as M-CS-TOF as the subrate increases which surpass CS-TOF obviously. These quantitative comparisons clearly indicate that the proposed method consistently surpasses CS-TOF and achieves the same effect with M-CS-TOF.

From both subjective observation and objective assessment in the above-mentioned, we can assert that the proposed 3D recover method is effective for discrete 3D targets imaging.

4.2. Reconstruction Performance of the Continuous Target

As is shown in Figure 7a, a 3D scene is designed as the space continuous target which is a 3D model of a tank T80. Observing from the left anterior of the tank, its entire length is more than 10 m long. The distance between the tank and RGI system is about 100 m and therefore the gate opening range of 15 m is designed. The 3D scene is modeled by a combination of two matrices which are showed as the original intensity map, Figure 7b and the original depth map, Figure 7c. In Figure 7c, colors correspond to the distance and the bar on the right side is added as an annotation.

The visual comparison of the imaging results of the T80 target scene by various methods is given in Figure 8. The measurements are in the presence of white Gaussian noise with SNR = 30 dB and all of the subrates are equal to 10%. Colors correspond to the distance from 100.0 m to 115.0 m and the bar on the right side is added as an annotation.

Figure 8a shows the reconstructed intensity maps of the target applying CS-TOF. There is a big error in terms of the pixel values and the structures of the target are incomplete between the result and the original intensity map. Figure 8c,e are the reconstructed intensity maps of the target applying M-CS-TOF and the proposed method, respectively. These figures show that due to the full use of the CS framework, high-quality intensity maps are obtained which not only describe the shape feature of the scene but also have the edge detail of the target. Figure 8b is the reconstructed depth map from the measurements applying CS-TOF. Unfortunately, the depth map is not only able to describe the integral outline of the target but also lacks the time-resolved ability. Figure 8d,f are the reconstructed depth maps of the continuous object applying M-CS-TOF and the proposed method. The object planes and edges are recovered well to demonstrate the depth and transverse resolution by determining the depth of each pixel accurately except for some acceptable noise. The errors in the edge of the object in Figure 8f result from some inaccurate pixels obtained in the optimization computation procedure, which will be eliminated as the subrate increases. The result indicates the feasibility of the proposed method.

Figure 9 provides the objective assessment of reconstruction qualities of the T80 target scene in the presence of white Gaussian noise with SNR = 30 dB. Figure 9a presents the plots of the PSNR values for reconstructed intensity maps on various subrates (from 2% to 14%) using different approaches.

It shows that PSNRs of the proposed method are always higher than those of the TOF-based methods as the subrates reach 10%. The plots corresponding to M-CS-TOF and the proposed method stay stubbornly high and the PSNRs are more than 30 dB so that high-quality intensity maps are achieved. However, the PSNRs of the CS-TOF reconstructing intensity maps are much lower which results from the error of the measurements used for reconstruction. The phenomenon that the echo pulses with different delays overlap each other always leads to significant reconstruction errors applying TOF-based methods. Figure 9b presents the plots of the SSIM values for reconstructed intensity maps on various subrates (from 2% to 14%) using different approaches. The plots corresponding to M-CS-TOF and the proposed method stay stubbornly high and the SSIMs are close to 1 when the subrates reach 10% so that the structures of the target and the reconstructed image are extraordinary similar. On the contrary, the SSIMs of the CS-TOF reconstructing intensity maps are much lower. With similar pattern to the PSNR results, the SSIMs further verify the effectiveness and advantage of the proposed method. Figure 9c presents the plots of the NMSE values for reconstructed depth maps on various subrates (from 2% to 14%) using different approaches. It shows that NMSEs of M-CS-TOF and the proposed method are extremely low especially when the subrates reach 10%. These two methods are feasible and high-quality depth maps are achieved. Figure 9d presents the plots of the SSIM values for reconstructed depth maps on various subrates (from 2% to 14%) using different approaches. The SSIMs are close to 1 when the subrates reach 10% except for CS-TOF. Both the pixel-wise and the structure-wise assessment demonstrates the reconstruction quality. Moreover, These quantitative comparisons clearly indicate that the proposed method consistently surpasses CS-TOF and achieves the comparative effect with M-CS-TOF.

From both subjective observation and objective assessment in the above-mentioned, we can assert that the proposed 3D recover method is effective for continuous 3D targets imaging.

Figure 10 provide the 3D reconstruction results of the “U & R” target with varying reflectance. The measurements are in the presence of white Gaussian noise with SNR = 30 dB and all of the subrates equal to 10%. Colors correspond to the distance from 20.0 m to 21.0 m and the bar on the right side is added as an annotation. The objects are with varying reflectance in the intensity map in Figure 10a,b is the depth map where two letters are at separated locations. The results demonstrate that the proposed method is effective for 3D imaging of targets with varying reflectance.

Table 2 is provided for comparison with the performance of several CS-based 3D imaging methods in terms of FWHM, sampling rate and the number of frames. To achieve a particular range resolution as is shown in the simulation above, the required parameters are given when each method is adopted respectively. The required width of the transmitting laser pulses is free and the required sampling frequency of the detector is lower because they are not directly relevant to the range resolution. It means that the proposed method can achieve range super-resolution to reconstruct the 3D scene when conventional methods are ineffective against the limitation of narrow-bandwidth pulsed lasers and high-frequency sampling. Moreover, a frame refers to an intensity image recovered from the measurements which are recorded from the total echo waveforms. Compared with conventional methods, the proposed method leads to smaller data sets and shorter measuring time that it can obtain specific range resolution with much fewer frames. Reconstructing only two frames to achieve 3D imaging reduces the computation expense resulting in a faster procedure.

5. Conclusions

In this paper, we presented a compressive sensing based electro-optic modulation method for nonscanning 3D laser imaging to overcome the limitations of traditional methods. Generally, the flying time values of laser pulses from multiple pixels do not possess linear properties which means they are not compressible. The proposed method transforms the flying time values into intensity values of the echo pulses through electro-optic modulation and peak values of the total echo pulses are recorded in CS framework. As a result, both the time and reflectivity information are compressible and can be reconstructed from much smaller sets of their linear projection. Two sets of CS measurements are obtained to recover the modulated images of the target when a monotone increasing and a monotone decreasing voltage functions are applied respectively. We formulate a minimization model integrated with a detail-preserving regularization and TV regularization for the inverse problem of image retrieval. Then we propose an ADMM imbedded algorithm to solve the unconstrainted regulated optimization problem efficiently. Finally, the 3D form of the target scene is completed as a combination of the intensity map and the depth map extracted from the two reconstructed images with gain modulation. The simulation results demonstrate that the proposed method is able to achieve high-quality 3D imaging for both the discrete target and the continuous target. More importantly, it has the advantage of range super-resolution with fewer required frames compared with CS-TOF method and M-CS-TOF method when echo pulses from targets with different ranges are indistinguishable. With fewer reconstructed frames in need, the proposed method requires less data for further processing, thus reducing the computation expense in reconstruction to achieve a faster performance. The research has important scientific value and is beneficial for the extensive application of 3D laser imaging such as low-cost LiDAR, high-speed image acquisition and remote sensing. A natural progression of this work is to build a realizable system for 3D imaging of various targets..

Author Contributions

Conceptualization, Y.A. and H.G.; methodology, Y.A.; software, Y.A.; validation, Y.A., H.G. and J.W.; formal analysis, Y.A.; investigation, Y.A.; resources, J.W.; data curation, Y.A.; writing—original draft preparation, Y.A.; writing—review and editing, Y.A.; visualization, Y.A.; supervision, Y.Z.; project administration, Y.Z.; funding acquisition, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank the handing Associate Editor and the anonymous reviewers for their valuable comments and suggestions for this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:

LiDAR	Laser Detection and Ranging
3D	Three-dimensional
CS	compressive sensing
TOF	Time-Of-Flight
ADMM	alternating direction method of multiplier
RGI	Range-Gated Imaging
NBF	NarrowBand filter
EOM	Electro-Optic Modulator
DMD	Digital Micromirror Device
APD	Avalanche PhotoDiode
SNR	Signal-to-Noise Ratio
FWHM	Full Width at Half Maximum
DCT	Discrete Cosine Transform
DWT	Discrete Wavelet Transfrom
RIP	Restricted Isometry Property
TV	Total Variation
TVAL3	Total Variation minimization based on Augmented Lagrangian and ALternating direction Algorithm
TGV	Generalization of TV
PSNR	Peak Signal to Noise Ratio
NMSE	Normalized Mean Squared Error
SSIM	Structural SIMilarity

References

Molebny, V.; Mcmanamon, P.; Steinvall, O.; Kobayashi, T.; Chen, W. Laser radar: Historical prospective—From the East to the West. Opt. Eng. 2016, 56, 031220. [Google Scholar] [CrossRef]
Li, L.; Wu, L.; Wang, X.; Dang, E. Gated viewing laser imaging with compressive sensing. Appl. Opt. 2012, 51, 2706–2712. [Google Scholar] [CrossRef] [PubMed]
Gao, H.; Zhang, Y.; Guo, H. Multihypothesis-Based Compressive Sensing Algorithm for Nonscanning Three-Dimensional Laser Imaging. IEEE J. STARS 2018, 11, 311–321. [Google Scholar] [CrossRef]
Colaco, A.; Kirmani, A.; Howland, G.A.; Howell, J.C.; Goyal, V.K. Compressive Depth Map Acquisition Using a Single Photon-Counting Detector: Parametric Signal Processing Meets Sparsity. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Providence, RI, USA, 16–21 June 2012; pp. 96–102. [Google Scholar] [CrossRef] [Green Version]
Gibson, G.M.; Sun, B.; Edgar, M.P.; Phillips, D.B.; Hempler, N.; Maker, G.T.; Malcolm, G.P.A.; Padgett, M.J. Real-time imaging of methane gas leaks using a single-pixel camera. Opt. Express 2017, 25, 2998–3005. [Google Scholar] [CrossRef] [PubMed]
Edgar, M.; Johnson, S.; Phillips, D.; Padgett, M. Real-time computational photon-counting LiDAR. Opt. Eng. 2018, 57, 031304. [Google Scholar] [CrossRef] [Green Version]
Sun, M.J.; Zhang, J.M. Single-Pixel Imaging and Its Application in Three-Dimensional Reconstruction: A Brief Review. Sensors 2019, 19, 732. [Google Scholar] [CrossRef] [Green Version]
Zhang, X.D.; Li, C.L.; Meng, Q.P.; Liu, S.J.; Zhang, Y.; Wang, J.Y. Infrared Image Super Resolution by Combining Compressive Sensing and Deep Learning. Sensors 2018, 18, 2587. [Google Scholar] [CrossRef] [Green Version]
Zhang, T.; Gao, K. MAP-MRF-based super-resolution reconstruction approach for coded aperture compressive temporal imaging. Appl. Sci. 2018, 8, 338. [Google Scholar] [CrossRef] [Green Version]
Donoho, D.L. Compressed sensing. IEEE T. Inform. Theory 2006, 52, 1289–1306. [Google Scholar] [CrossRef]
Candès, E.J. Compressive sampling. In Proceedings of the International Congress of Mathematicians, Madrid, Spain, 22–30 August 2006; Volume 3, pp. 1433–1452. [Google Scholar]
Rani, M.; Dhok, S.B.; Deshmukh, R.B. A Systematic Review of Compressive Sensing: Concepts, Implementations and Applications. IEEE Access 2018, 6, 4875–4894. [Google Scholar] [CrossRef]
Gan, H.P.; Xiao, S.; Zhang, T.; Zhang, Z.M.; Li, J.; Gao, Y. Chaotic Pattern Array for Single-Pixel Imaging. Electronics 2019, 8, 536. [Google Scholar] [CrossRef] [Green Version]
Duarte, M.F.; Davenport, M.A.; Takhar, D.; Laska, J.N.; Sun, T.; Kelly, K.F.; Baraniuk, R.G. Single-pixel imaging via compressive sampling. IEEE Signal Proc. Mag. 2008, 25, 83–91. [Google Scholar] [CrossRef] [Green Version]
Howland, G.A.; Dixon, P.B.; Howell, J.C. Photon-counting compressive sensing laser radar for 3D imaging. Appl. Opt. 2011, 50, 5917–5920. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Sun, M.J.; Edgar, M.P.; Gibson, G.M.; Sun, B.; Radwell, N.; Lamb, R.; Padgett, M.J. Single-pixel three-dimensional imaging with time-based depth resolution. Nat. Commun. 2016, 7. [Google Scholar] [CrossRef] [PubMed]
Howland, G.A.; Lum, D.J.; Ware, M.R.; Howell, J.C. Photon counting compressive depth mapping. Opt. Express 2013, 21, 23822–23837. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Busck, J.; Heiselberg, H. Gated viewing and high-accuracy three-dimensional laser radar. Appl. Opt. 2004, 43, 4705–4710. [Google Scholar] [CrossRef] [PubMed]
Laurenzis, M.; Christnacher, F.; Monnin, D. Long-range three-dimensional active imaging with superresolution depth mapping. Opt. Lett. 2007, 32, 3146–3148. [Google Scholar] [CrossRef]
Jin, C.; Sun, X.; Zhao, Y.; Zhang, Y.; Liu, L. Gain-modulated three-dimensional active imaging with depth-independent depth accuracy. Opt. Lett. 2009, 34, 3550–3552. [Google Scholar] [CrossRef]
Tsagkatakis, G.; Woiselle, A.; Tzagkarakis, G.; Bousquet, M.; Starck, J.L.; Tsakalides, P. Multireturn compressed gated range imaging. Opt. Eng. 2015, 54, 031106. [Google Scholar] [CrossRef] [Green Version]
Wang, X.; Li, Y.; Zhou, Y. Multi-pulse time delay integration method for flexible 3D super-resolution range-gated imaging. Opt. Express 2015, 23, 7820–7831. [Google Scholar] [CrossRef]
Laurenzis, M.; Becher, E. Three-dimensional laser-gated viewing with error-free coding. Opt. Eng. 2018, 57, 7. [Google Scholar] [CrossRef]
Chen, Z.; Liu, B.; Liu, E.; Peng, Z. Electro-optic modulation methods in range-gated active imaging. Appl. Opt. 2016, 55, A184–A190. [Google Scholar] [CrossRef]
An, Y.; Zhang, Y.; Guo, H.; Wang, J. Compressive Sensing-Based Three-Dimensional Laser Imaging With Dual Illumination. IEEE Access 2019, 7, 25708–25717. [Google Scholar] [CrossRef]
Fade, J.; Perrotin, E.; Bobin, J. Polarizer-free two-pixel polarimetric camera by compressive sensing. Appl. Opt. 2018, 57, B102–B113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Zhang, Y.L.; Suo, J.L.; Wang, Y.W.; Dai, Q.H. Doubling the pixel count limitation of single-pixel imaging via sinusoidal amplitude modulation. Opt. Express 2018, 26, 6929–6942. [Google Scholar] [CrossRef] [PubMed]
Chen, Z.; Liu, B.; Liu, E.; Peng, Z. Adaptive Polarization-Modulated Method for High-Resolution 3D Imaging. IEEE Photonic Technol. Lett. 2016, 28, 295–298. [Google Scholar] [CrossRef]
Gao, H.; Zhang, Y.-M.; Guo, H.-C. A compressive sensing algorithm using truncated SVD for three-dimensional laser imaging of space-continuous targets. J. Mod. Optic. 2016, 63, 2166–2172. [Google Scholar] [CrossRef]
Duarte, M.F.; Baraniuk, R.G. Kronecker Compressive Sensing. IEEE Trans. Image Process. 2012, 21, 494–504. [Google Scholar] [CrossRef]
Buyssens, P.; Daisy, M.; Tschumperle, D.; Lezoray, O. Exemplar-Based Inpainting: Technical Review and New Heuristics for Better Geometric Reconstructions. IEEE Trans. Image Process. 2015, 24, 1809–1824. [Google Scholar] [CrossRef]
Duarte, M.F.; Sarvotham, S.; Baron, D.; Wakin, M.B.; Baraniuk, R.G. Distributed compressed sensing of jointly sparse signals. In Proceedings of the Conference Record of the Thirty-Ninth Asilomar Conference onSignals, Systems and Computers, Pacific Grove, CA, USA, 30 October–2 November 2005; pp. 1537–1541. [Google Scholar]
Candes, E.J.; Wakin, M.B. An introduction to compressive sampling. IEEE Signal Proc. Mag. 2008, 25, 21–30. [Google Scholar] [CrossRef]
Bhattacharjee, T.; Maity, S.P. Progressive and hierarchical share-in-share scheme over cloud. J. Inf. Secur. Appl. 2019, 46, 108–120. [Google Scholar] [CrossRef]
Bian, L.H.; Suo, J.L.; Dai, Q.H.; Chen, F. Experimental comparison of single-pixel imaging algorithms. J. Opt. Soc. Am. A 2018, 35, 78–87. [Google Scholar] [CrossRef] [PubMed]
Guo, W.; Qin, J.; Yin, W. A New Detail-Preserving Regularization Scheme. SIAM J. Imaging Sci. 2014, 7, 1309–1334. [Google Scholar] [CrossRef]
Boyd, S.; Parikh, N.; Chu, E.; Peleato, B.; Eckstein, J. Distributed optimization and statistical learning via the alternating direction method of multipliers. Found. Trends^® Mach. Learn. 2011, 3, 1–122. [Google Scholar] [CrossRef]
Li, C. An Efficient Algorithm for Total Variation Regularization with Applications to the Single Pixel Camera and Compressive Sensing. Master’s Thesis, Department of Computational and Applied Mathematics, Houston, TX, USA, 2009. [Google Scholar]
Li, F.; Chen, H.; Pediredla, A.; Yeh, C.; He, K.; Veeraraghavan, A.; Cossairt, O. CS-ToF: High-resolution compressive time-of-flight imaging. Opt. Express 2017, 25, 31096–31110. [Google Scholar] [CrossRef] [Green Version]
Kirmani, A.; Colaço, A.; Wong, F.N.; Goyal, V.K. Exploiting sparsity in time-of-flight range acquisition using a single time-resolved sensor. Opt. Express 2011, 19, 21485–21507. [Google Scholar] [CrossRef]
Czajkowski, K.M.; Pastuszczak, A.; Kotynski, R. Single-pixel imaging with Morlet wavelet correlated random patterns. Sci. Rep. 2018, 8, 8. [Google Scholar] [CrossRef] [Green Version]
Marques, E.C.; Maciel, N.; Naviner, L.; Cai, H.; Yang, J. A Review of Sparse Recovery Algorithms. IEEE Access 2019, 7, 1300–1322. [Google Scholar] [CrossRef]
Wang, Z.; Bovik, A.C.; Sheikh, H.R.; Simoncelli, E.P. Image quality assessment: From error visibility to structural similarity. IEEE T. Image Process. 2004, 13, 600–612. [Google Scholar] [CrossRef] [Green Version]

Figure 1. Schematic diagram of the proposed compressive sensing RGI system.

Figure 2. The time sequence of gain modulation for the RGI system: (a) time sequence of receiver gate, applied voltage of the EOM, gain function, laser pulse, and received pulse of APD; (b) received pulse and its subpulses with gain function

g_{i n} (t)

; (c) received pulse and its subpulses with gain function

g_{d e} (t)

g_{i n} (t)

; (c) received pulse and its subpulses with gain function

g_{d e} (t)

Figure 3. Schematic diagram of electro-optic modulation.

Figure 4. 3D scene ranged from 20.3 m to 20.9 m: (a) the conceptual graph of the “U & R” target; (b) the original intensity map; (c) the original depth map.

Figure 5. The 3D reconstruction results of the “U & R” target scene by various methods: (a) the reconstructed intensity map by CS-TOF; (b) the reconstructed depth map by CS-TOF; (c) the reconstructed intensity map by M-CS-TOF; (d) the reconstructed depth map by M-CS-TOF; (e) the reconstructed intensity map by the proposed method; (f) the reconstructed depth map by the proposed method.

Figure 6. The objective assessment of reconstruction qualities of the “U & R” target scene using different approaches: (a) the plots of PSNR for reconstructed intensity maps as a function of subrates; (b) the plots of SSIM for reconstructed intensity maps as a function of subrates; (c) the plots of NMSE for reconstructed depth maps as a function of subrates; (d) the plots of SSIM for reconstructed depth maps as a function of subrates.

Figure 7. 3D scene ranged from 100 m to 115 m: (a) a 3D model of a tank T80 in the target scene; (b) the original intensity map; (c) the original depth map.

Figure 8. The 3D reconstruction results of the T80 target scene by various methods: (a) the reconstructed intensity map by CS-TOF; (b) the reconstructed depth map by CS-TOF; (c) the reconstructed intensity map by M-CS-TOF; (d) the reconstructed depth map by M-CS-TOF; (e) the reconstructed intensity map by the proposed method; (f) the reconstructed depth map by the proposed method.

Figure 9. The objective assessment of reconstruction qualities of the T80 target scene using different approaches: (a) the plots of PSNR for reconstructed intensity maps as a function of subrates; (b) the plots of SSIM for reconstructed intensity maps as a function of subrates; (c) the plots of NMSE for reconstructed depth maps as a function of subrates; (d) the plots of SSIM for reconstructed depth maps as a function of subrates.

Figure 10. The 3D reconstruction results of the “U & R” target with varying reflectance: (a) the reconstructed intensity map; (b) the reconstructed depth map.

Table 1. Parameters of the simulation experiments.

Expression	Value
Wavelength	905 nm
Sampling Rate of APD	1 GHz
Peak Power of Transmitter Pulse	70 W
FWHM	10 ns
Efficiency of Optical Transmitting System	0.9
Efficiency of Optical Receiving System	0.9
Single-pass Atmospheric Transmittance	0.98

Table 2. Performance of several CS-based 3D imaging methods.

Imaging Methods	CS-TOF	M-CS-TOF	Proposed
FWHM	resolution related	free	free
Sampling rate	higher	higher	lower
Sets of measurements	some	dozens	2
Number of frames	some	dozens	2
Execution time	medium	long	short

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

An, Y.; Zhang, Y.; Guo, H.; Wang, J. Compressive Sensing Based Three-Dimensional Imaging Method with Electro-Optic Modulation for Nonscanning Laser Radar. Symmetry 2020, 12, 748. https://doi.org/10.3390/sym12050748

AMA Style

An Y, Zhang Y, Guo H, Wang J. Compressive Sensing Based Three-Dimensional Imaging Method with Electro-Optic Modulation for Nonscanning Laser Radar. Symmetry. 2020; 12(5):748. https://doi.org/10.3390/sym12050748

Chicago/Turabian Style

An, Yulong, Yanmei Zhang, Haichao Guo, and Jing Wang. 2020. "Compressive Sensing Based Three-Dimensional Imaging Method with Electro-Optic Modulation for Nonscanning Laser Radar" Symmetry 12, no. 5: 748. https://doi.org/10.3390/sym12050748

APA Style

An, Y., Zhang, Y., Guo, H., & Wang, J. (2020). Compressive Sensing Based Three-Dimensional Imaging Method with Electro-Optic Modulation for Nonscanning Laser Radar. Symmetry, 12(5), 748. https://doi.org/10.3390/sym12050748

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Compressive Sensing Based Three-Dimensional Imaging Method with Electro-Optic Modulation for Nonscanning Laser Radar

Abstract

1. Introduction

2. System Description

2.1. Proposed System Setup

2.2. Proposed System Model

3. CS-Based Electro-Optic Modulation Method for 3D Imaging

3.1. Compressive Sensing Theory

3.2. Peak Value Obtained for 3D Reconstruction

3.3. 3D Image Reconstruction

4. Simulation Results and Discussion

4.1. Reconstruction Performance of the Discrete Target

4.2. Reconstruction Performance of the Continuous Target

5. Conclusions

Author Contributions

Funding

Acknowledgments

Conflicts of Interest

Abbreviations

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI