## 1.

## Introduction

Laser range-gated imaging can penetrate a strong scattering medium and capture the image of an object at a certain distance by applying active laser illumination, changing the gate delay time, and varying the gate width. It has been widely used in areas such as military reconnaissance, monitoring, underwater detection, and search and rescue. For cases when the line of sight is blocked by corners of walls or other building obstructions within the environment of city streets and buildings (as shown in Fig. 1), traditional imaging techniques and equipment are unable to achieve imaging beyond the direct field of view. Non-line-of-sight imaging can bypass the corners of walls or other obstructions to achieve imaging beyond the direct field of view, which can be used to reduce casualties in modern urban street fighting, and to increase battlefield awareness and operational efficiency.

In 2006 and 2007, researchers at the Research Institute for Optronics and Pattern Recognition of Germany and Swedish Defence Research Agency (FOI)^{1}^{,}^{2} of Sweden carried out joint research on short-wavelength infrared laser range-gated imaging for ground applications in monostatic and bistatic configurations. A number of test targets (test panels, vehicles, and soldiers with weapons) were included in the tests, and windows in buildings and cars, traffic signs, and vehicle surfaces were utilized as the intermediary reflective surfaces. In 2011, the FOI group^{3} proposed the loss factor model through further study of the effect of the characteristics (including the shape and the reflection characteristics) of the intermediary reflective surface on laser illuminance attenuation and measured the loss factors with different intermediary reflective surfaces. Though their work was insightful, their model is not effective for describing the imaging process. Since 2007, the Media Lab at the Massachusetts Institute of Technology (MIT)^{4} has been studying transient imaging based on femtosecond laser range-gated imaging technology and proposed an algorithm for extracting object information from multireflection, realizing a means of seeing around the corner. In 2012, researchers at MIT^{5} realized a three-dimensional reconstruction of hidden targets utilizing diffuse reflections of ultrashort laser pulses using an intermediary reflective surface. Since 2012, researchers at the Beijing Institute of Technology (BIT)^{6}^{,}^{7} have been studying non-line-of-sight imaging. The BIT group uses a nanosecond laser range-gated imaging system to realize non-line-of-sight imaging of an USAF 1951 standard target at different distances with glass and ceramic tile as the respective intermediary reflective surfaces. These existing researches show the potential application prospects of non-line-of-sight imaging in areas such as urban warfare, counter-terrorism, and disaster relief.

Since non-line-of-sight imaging based on laser range-gated imaging is affected not only by the parameters of the laser range-gated imaging system but also by the reflection characteristics of the intermediary reflective surface and the scene, construction of an optical imaging model is important to the development and application of non-line-of-sight imaging. This paper describes the construction of an image contrast model of non-line-of-sight imaging from the viewpoint of a two-dimensional illumination distribution on the photosensitive surface, in order to analyze the performance of a non-line-of-sight imaging system utilizing a laser range-gated imaging system and to provide a theoretical basis for actual non-line-of-sight imaging system design.

## 2.

## Principle of Non-Line-of-Sight Imaging

The key to non-line-of-sight imaging is to make it possible to see a street-corner target or an indoor target by utilizing an intermediary reflective surface, as illustrated in Fig. 2. A is the observer, B is the target under observation, and C is the intermediary reflective surface. B is outside the range of objects that can be directly imaged by A, so that C is utilized to realize non-line-of-sight imaging of target B.

A laser range-gated imaging system, consisting of a nanosecond laser and an intensified charged coupled device (ICCD) imaging system, is located at location A. The ICCD imaging system is focused on location B. The beam from the nanosecond laser from location A experiences first-order reflections (consisting of first-order specular reflection and first-order diffuse reflections, where the first-order diffuse reflection contains first-order backscattering and first-order forward-scattering) at the intermediary reflective surface C. The first-order specular reflection and first-order forward-scattering illuminate target B and then are reflected at target B in a Lambertian fashion (the second-order reflection). First-order backscattering is received by the ICCD imaging system, forming a defocused image and reducing the image contrast of target B. The second-order Lambertian reflection strikes the intermediary reflective surface C and experiences third-order reflections (consisting of third-order specular reflection, third-order backscattering, and third-order forward-scattering) at location C. Third-order specular reflection is received by the ICCD imaging system, forming a focused image of target B. Like the first-order backscattering, the third-order forward-scattering is also received by the ICCD imaging system, reducing the contrast of the image of target B.

For certain intermediary reflective surfaces (like ceramic tile), the effective target signal illuminance (third-order specular reflection) may be overwhelmed by the first-order backscattering and third-order forward-scattering, resulting in the target being undetectable. The first-order backscattering is apparently larger than the third-order forward-scattering, so the first-order backscattering is the dominant factor reducing the contrast of the image of target B. Thus, a better non-line-of-sight image of target B can be obtained by filtering out the first-order backscattering, which can be realized by laser range-gated imaging. Meanwhile, an image intensifying technique is also needed because the effective target signal illuminance (third-order specular reflection) is typically weak. The intermediary reflective surface is also a very important factor in non-line-of-sight imaging, which can be generalized to include a full specular reflector like glass (1), a partial specular reflector like ceramic tile (2), and an ideal diffuse surface (3). The study in this paper is mainly carried out on intermediary reflective surface types 1 and 2.

Construction of an optical imaging model of the non-line-of-sight imaging process, which is based on the parameters of the intermediary reflective surface, the parameters of the laser range-gated imaging system, and the properties of the target, is vital for the theoretical analysis and system design of non-line-of-sight imaging based on laser range-gated imaging. This model assumes that the imaging will occur at night or in low-light conditions to reduce the impact of environmental illuminance.

## 3.

## Image Contrast Model of Non-Line-of-Sight Imaging

As shown in Fig. 3, the laser beam is projected onto an intermediary reflective surface forming the light spot ${S}_{1}$ and experiences first-order reflections; then the first-order reflections are projected onto the target forming light spot ${S}_{2}$. The distance from the laser to the geometric center of ${S}_{1}$ is ${l}_{1}$. The distance from the geometric center of ${S}_{1}$ to the geometric center of ${S}_{2}$ is ${l}_{2}$. The vector tracing the direction of the incident light from ${S}_{1}$ to ${S}_{2}$ makes an angle $\alpha $ with ${\mathbf{n}}_{1}$ (the vector normal to ${S}_{1}$) and an angle $\beta $ with ${\mathbf{n}}_{2}$ (the vector normal to ${S}_{2}$). The optical axis of the laser and the optical axis of the ICCD are assumed to be coincident, and the divergence angle of the laser and the field view angle of the ICCD are also assumed to be coincident. In this section, the laser transmission characteristics (attenuation and distribution) in the non-line-of-sight imaging process will be analyzed.

## 3.1.

### Pulsed Laser and Range-Gated Imaging System

## (1)

$${I}_{0}(t)=\frac{P(t)}{{\mathrm{\Omega}}_{\text{laser}}}=\frac{\mathrm{\Phi}}{{t}_{p}\xb7{\mathrm{\Omega}}_{\text{laser}}}\xb7\mathrm{rect}\left[\frac{1}{{t}_{p}}\right(t+\frac{{t}_{p}}{2}\left)\right]\phantom{\rule{0ex}{0ex}}={I}_{0}\xb7\mathrm{rect}\left[\frac{1}{{t}_{p}}\right(t+\frac{{t}_{p}}{2}\left)\right],$$## (2)

$${\mathrm{\Omega}}_{\text{laser}}=\pi \text{\hspace{0.17em}}{\mathrm{sin}}^{2}\text{\hspace{0.17em}}\theta ,$$## 3.2.

### Reflection Characteristics of the Intermediary Reflective Surface

The reflection characteristics (consisting of specular reflection characteristics and diffuse reflection characteristics) of the intermediary reflective surface are among the important factors that can influence the non-line-of-sight imaging results. Without loss of generality, the intermediary reflective surface is assumed to be free of absorption, so that its reflection characteristics can be described by the specular reflectivity ${\rho}_{s}$ and the diffuse reflectivity ${\rho}_{d}$. ${\rho}_{s}$ is a constant, ${\rho}_{d}$ obeys Lambertian reflection, and ${\rho}_{s}=1-{\rho}_{d}$. The reflectance distribution function is described by

with## 3.3.

### Target Signal Illuminance and Disturbance Caused by the Intermediary Reflective Surface under CW-Laser Illumination

## 3.3.1.

#### Illuminance on the target plane ${S}_{2}$

The illuminance on the target plane ${S}_{2}$ is denoted as ${E}_{2}$, which consists of two parts: ${E}_{2s}$ is caused by the first-order specular reflection of the laser and ${E}_{2d}$ is caused by the first-order forward-scattering. The illuminance on ${S}_{1}$ should be obtained first.

### First-order specular reflection

The area illuminated on ${S}_{2}$ by the laser first-order specularly reflected by ${S}_{1}$ can be written as

## (7)

$${s}_{2}=\frac{{({l}_{1}+{l}_{2})}^{2}\xb7{\mathrm{\Omega}}_{\text{laser}}}{\mathrm{cos}\text{\hspace{0.17em}}\beta}.$$Therefore, ${E}_{2s}$ can be calculated as

## (8)

$${E}_{2s}=\frac{P\xb7{\rho}_{s}}{{s}_{2}}=\frac{P\xb7{\rho}_{s}}{{({l}_{1}+{l}_{2})}^{2}\xb7{\mathrm{\Omega}}_{\text{laser}}}\xb7\mathrm{cos}\text{\hspace{0.17em}}\beta .$$The target distance $l$ is equal to the sum of ${l}_{1}$ and ${l}_{2}$. If $l$ is much larger than the square root of ${s}_{2}$, and ${S}_{2}$ is illuminated by the vertically incident laser beam ($\mathrm{cos}\text{\hspace{0.17em}}\beta =1$), Eq. (8) can be further simplified.

### First-order forward-scattering

To facilitate the analysis and calculations, Fig. 3 is expanded along the optical path. Assuming that the optical axis of the laser coincides with the optical axis of the ICCD imaging system, a Cartesian coordinate system is established with the geometric center of the lens as the Cartesian coordinate origin. The optical axis is along the $z$ axis, and the $y$ axis and the $x$ axis are, respectively, parallel to the detector array’s length and width directions, as is shown in Fig. 4. $\phi $ is the angle between the $x$ axis and the projection of ${\mathbf{n}}_{1}$ on the $xoy$ plane, and ${\mathbf{n}}_{2}$ points in the direction opposite to the direction of the $z$ axis. ${S}_{1}$ and ${S}_{2}$ are, respectively, projections of the laser beam on the intermediary reflective surface and on the target.

According to the geometric relationship, the projection equations for ${S}_{1}$ and ${S}_{2}$ on the $xoy$ plane are, respectively, given by

## (9)

$${x}^{2}+{y}^{2}-{\mathrm{tg}}^{2}\theta \xb7{[\mathrm{tg}\alpha \xb7(x\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\phi -y\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\phi )+{l}_{1}]}^{2}\le 0,$$${P}_{1}({x}_{1},{y}_{1},{z}_{1})$ and ${P}_{2}({x}_{2},{y}_{2},{z}_{2})$ are random points, respectively, on ${S}_{1}$ and ${S}_{2}$, and ${d}_{S1}$ is a microsurface element around ${P}_{1}({x}_{1},{y}_{1},{z}_{1})$ (as shown in Fig. 5).

## (11)

$$d{E}_{2d}=\frac{{L}_{1}\xb7d{s}_{1}}{{l}_{p2}^{2}}\xb7\mathrm{cos}\text{\hspace{0.17em}}{\gamma}_{2}\xb7\mathrm{cos}\text{\hspace{0.17em}}{\eta}_{2}\phantom{\rule{0ex}{0ex}}=\frac{{E}_{1}\xb7{\rho}_{d}\xb7\mathrm{cos}\text{\hspace{0.17em}}{\eta}_{2}\xb7\mathrm{cos}\text{\hspace{0.17em}}{\gamma}_{2}}{\pi \xb7{l}_{p2}^{2}\xb7\mathrm{cos}\text{\hspace{0.17em}}\alpha}\mathrm{d}{x}_{1}\mathrm{d}{y}_{1}.$$## (12)

$$\mathrm{cos}\text{\hspace{0.17em}}{\gamma}_{2}=\frac{\mathrm{tg}\alpha \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\phi ({x}_{2}-{x}_{1})+\mathrm{tg}\text{\hspace{0.17em}}\alpha \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\phi ({y}_{2}-{y}_{1})-({z}_{2}-{z}_{1})}{\sqrt{{({x}_{2}-{x}_{1})}^{2}+{({y}_{2}-{y}_{1})}^{2}+{({z}_{2}-{z}_{1})}^{2}}\sqrt{1+{\mathrm{tg}}^{2}\alpha}}(0\le {\gamma}_{2}\le \frac{\pi}{2}),$$## (13)

$$\mathrm{cos}\text{\hspace{0.17em}}{\eta}_{2}=\frac{-({z}_{2}-{z}_{1})}{\sqrt{{({x}_{2}-{x}_{1})}^{2}+{({y}_{2}-{y}_{1})}^{2}+{({z}_{2}-{z}_{1})}^{2}}}(0\le {\eta}_{2}\le \frac{\pi}{2}),$$Therefore,

## (15)

$${E}_{2d}={\iint}_{{S}_{1}}d{E}_{2d}\phantom{\rule[-0.0ex]{2em}{0.0ex}}(\alpha +\theta <\frac{\pi}{2}).$$In summary,

## 3.3.2.

#### Target signal illuminance on the photosensitive surface

## (17)

$${E}_{su}=\frac{1}{4}\xb7\pi \xb7{L}_{u}\xb7{\tau}_{o}\xb7{\rho}_{s}\xb7{\left(\frac{{D}_{0}}{{f}^{\prime}}\right)}^{2}/[1+\frac{1}{4}\xb7{\left(\frac{{D}_{0}}{{f}^{\prime}}\right)}^{2}]\phantom{\rule{0ex}{0ex}}=\frac{\pi \xb7{L}_{u}\xb7{\tau}_{o}\xb7{\rho}_{s}}{1+4{F}^{2}}.$$## 3.3.3.

#### Disturbance illuminance on the photosensitive surface caused by scattering of the intermediary reflective surface

The disturbance illuminance on the photosensitive surface caused by scattering of the intermediary reflective surface consists of two parts: ${E}_{d}$ (third-order forward-scattering) and ${E}_{\text{back}}$ (first-order backscattering). Since the intermediary reflective surface is not on the target plane, each point on ${S}_{1}$ corresponds to a diffuse plaque on the photosensitive surface, and the diffuse plaque can be expressed by a point spread function. A Cartesian coordinate system is established with the geometric center of the lens as a Cartesian coordinate origin. $P(0,0,{l}_{1})$ is the intersection of the intermediary reflective surface and the $z$ axis, and $P$’ is the image point of $P$, as shown in Fig. 6.

According to geometrical optics, the center of the diffuse plaque for point $({x}_{1},{y}_{1})$ on ${S}_{1}$ is given by

## (18)

$$({x}_{0}^{\prime},{y}_{0}^{\prime})=(-n{l}^{\prime}{x}_{1}/{l}_{1},-n{l}^{\prime}{y}_{1}/{l}_{1}),$$The radius of the diffuse plaque can be written as

## (20)

$${r}_{0}(z)=\frac{1}{2}\left\{\right[1-\frac{(z-n\xb7{f}^{\prime})\xb7{l}^{\prime}}{z{f}^{\prime}}]\xb7{D}_{0}\}=\frac{1}{2}\left\{\right[1-\frac{({l}_{1}+\mathrm{tg}\alpha \xb7{x}_{1}-n\xb7{f}^{\prime})\xb7{l}^{\prime}}{({l}_{1}+\mathrm{tg}\alpha \xb7{x}_{1}){f}^{\prime}}]\xb7{D}_{0}\},$$The point spread function of the diffuse plaque for point $({x}_{1},{y}_{1})$ can be expressed as

## (21)

$$h({x}^{\prime},{y}^{\prime},{x}_{1},{y}_{1},{l}_{1})=\mathrm{circ}\left[\frac{\sqrt{{({x}^{\prime}-{x}_{0}^{\prime})}^{2}+{({y}^{\prime}-{y}_{0}^{\prime})}^{2}}}{{r}_{0}(z)}\right],$$## (22)

$$\mathrm{circ}(r)=\mathrm{circ}\left(\sqrt{{x}^{2}+{y}^{2}}\right)=\{\begin{array}{cc}1& \sqrt{{x}^{2}+{y}^{2}}\le 1\\ 0& \text{otherwise}\end{array}.$$### First-order backscattering illuminance on the photosensitive surface

According to Eq. (21), the first-order backscattering ${E}_{\text{back}}$ can be expressed as

## (23)

$${E}_{\text{back}}=\frac{{E}_{1}\xb7{\rho}_{d}\xb7{\tau}_{o}\xb7{\mathrm{\Omega}}_{p}\text{\hspace{0.17em}}\mathrm{cos}\alpha \text{\hspace{0.17em}}}{{\pi}^{2}}\phantom{\rule{0ex}{0ex}}\xb7{\iint}_{{S}_{1}}\frac{1}{{r}_{0}^{2}(z)}\mathrm{circ}\left[\frac{\sqrt{{({x}^{\prime}-{x}_{0})}^{2}+{({y}^{\prime}-{y}_{0})}^{2}}}{{r}_{0}(z)}\right]\mathrm{d}{s}_{1},$$## (24)

$${\mathrm{\Omega}}_{p}=\pi {D}_{0}^{2}/4{({l}_{1}+{x}_{1}\xb7\mathrm{tan}\text{\hspace{0.17em}}\alpha )}^{2}.$$${\mathrm{\Omega}}_{p}$ is the solid angle defined by the rays emanating from point ${P}_{1}({x}_{1},{y}_{1},{z}_{1})$ to the aperture defined by the lens. In fact, since the distances between the camera and every point on ${S}_{1}$ are not the same, the value of ${\mathrm{\Omega}}_{p}$ for each of the points is not the same. As the laser divergence angle $\theta $ is very small, an approximate treatment for ${\mathrm{\Omega}}_{p}$ is made here.

### Third-order forward-scattering illuminance on the photosensitive surface

## (25)

$$d{E}_{1}^{\prime}=\frac{{L}_{2}\xb7d{s}_{2}}{{l}_{{p}_{1}}^{2}}\xb7\mathrm{cos}\text{\hspace{0.17em}}{\eta}_{1}\xb7\mathrm{cos}\text{\hspace{0.17em}}{\gamma}_{1}=\frac{{E}_{2}\xb7{\rho}_{\mathrm{max}}\xb7\mathrm{cos}\text{\hspace{0.17em}}{\eta}_{1}\xb7\mathrm{cos}\text{\hspace{0.17em}}{\gamma}_{1}}{\pi \xb7{l}_{{p}_{1}}^{2}}\mathrm{d}{x}_{2}\mathrm{d}{y}_{2}.$$In Fig. 5, ${d}_{S2}$ is a microsurface element around ${P}_{2}({x}_{2},{y}_{2},{z}_{2})$. $d{E}_{1}^{\prime}$ is the illuminance on ${P}_{2}({x}_{2},{y}_{2},{z}_{2})$ caused by ${d}_{S2}$. The illuminance is inversely proportional to the square of the distance as expressed in Eq. (25), where ${\gamma}_{1}$ is the angle between vectors $\overrightarrow{{P}_{1}{P}_{2}}$ and ${\mathbf{n}}_{1}$ and ${\eta}_{1}$ is the angle between vectors $\overrightarrow{{P}_{1}{P}_{2}}$ and ${\mathbf{n}}_{2}$. Then we obtain

## (26)

$${E}_{1}^{\prime}={\iint}_{{S}_{2}}d{E}_{1}^{\prime}\phantom{\rule[-0.0ex]{2em}{0.0ex}}(\alpha +\theta <\pi /2).$$Since ${I}_{p}$ can be calculated as

## (27)

$${I}_{p}(x,y,z)={L}_{p}\xb7\mathrm{cos}\text{\hspace{0.17em}}\alpha \xb7d{s}_{1}^{\prime}={E}_{1}^{\prime}\xb7{\rho}_{d}\xb7\mathrm{cos}\text{\hspace{0.17em}}\alpha \xb7d{s}_{1}^{\prime}/\pi ,$$## (28)

$$d{E}_{d}=\frac{{E}_{1}^{\prime}\xb7{\rho}_{d}\xb7{\tau}_{o}\xb7{\mathrm{\Omega}}_{p}\xb7\mathrm{cos}\text{\hspace{0.17em}}\alpha}{{\pi}^{2}\xb7{r}_{0}^{2}(z)}d{s}_{1}^{\prime}.$$In summary,

## 3.4.

### Image Contrast Model of Non-Line-of-Sight Imaging

Non-line-of-sight imaging can be achieved with laser range-gated imaging with a pulsed laser; therefore it is necessary to make a distinction between third-order reflection and first-order backscattering.

## 3.4.1.

#### Illuminance on the photosensitive surface under pulsed laser illumination

Under pulsed laser illumination, the illuminance on the photosensitive surface caused by a single pulse is a function of time $t$.

## (30)

$${E}_{su}(t)={E}_{su}\xb7\mathrm{rect}\left\{\frac{1}{{t}_{p}}\right[t+\frac{{t}_{p}}{2}-\frac{2({l}_{1}+{l}_{2})}{c}\left]\right\}\phantom{\rule[-0.0ex]{2em}{0.0ex}}(u=t,b),$$## 3.4.2.

#### Signal and noise under pulsed laser illumination

Considering the time sequence of the laser range-gated imaging system and the discrete nature of the detector unit, the illuminance on the detector unit $(i,j)$ caused by the target and background within a pulse period can be calculated as

## (33)

$${Q}_{u}(i,j)=ab\xb7{E}_{su}\xb7{\int}_{{t}_{d}}^{{t}_{d}+{t}_{g}}\mathrm{rect}\left\{\frac{1}{{t}_{p}}\right[t+\frac{{t}_{p}}{2}-\frac{2({l}_{1}+{l}_{2})}{c}\left]\right\}\mathrm{d}t,$$The illuminance on the detector unit $(i,j)$ caused by scattering of the intermediary reflective surface can be expressed, respectively, as

## (34)

$${Q}_{\text{back}}(i,j)={\iint}_{ds}{E}_{\text{back}}[(i+1/2)a,(j+1/2)b]\mathrm{d}{x}^{\prime}\mathrm{d}{y}^{\prime}\phantom{\rule{0ex}{0ex}}{\int}_{{t}_{d}}^{{t}_{d}+{t}_{g}}\mathrm{rect}\left[\frac{1}{{t}_{p}}\right(t+\frac{{t}_{p}}{2}-\frac{2{l}_{1}}{c}\left)\right]\mathrm{d}t,$$## (35)

$${Q}_{d}(i,j)={\iint}_{ds}{E}_{d}[(i+1/2)a,(j+1/2)b]\mathrm{d}{x}^{\prime}\mathrm{d}{y}^{\prime}\phantom{\rule{0ex}{0ex}}\xb7{\int}_{{t}_{d}}^{{t}_{d}+{t}_{g}}\mathrm{rect}\left\{\frac{1}{{t}_{p}}\right[t+\frac{{t}_{p}}{2}-\frac{2({l}_{1}+{l}_{2})}{c}\left]\right\}\mathrm{d}t.$$We simplify the imaging system noise and environmental noise by combining them into an equivalent incident illuminance ${E}_{\mathrm{eq}}$, which is used as shown in the following equation:

The total disturbance of the non-line-of-sight imaging system can be expressed as

## 3.4.3.

#### Apparent contrast of the image on a photosensitive surface

The apparent contrast of the image on a photosensitive surface can be described by

## (38)

$$C(i,j,{l}_{1},{l}_{2})=\left|\frac{({Q}_{t}+{Q}_{n})-({Q}_{b}+{Q}_{n})}{{Q}_{t}+{Q}_{b}+2{Q}_{n}}\right|=\left|\frac{{\rho}_{t}-{\rho}_{b}}{{\rho}_{t}+{\rho}_{b}}\right|\phantom{\rule{0ex}{0ex}}\xb7\left|\frac{1}{1+{Q}_{n}/{Q}_{tb}}\right|={C}_{0}(i,j){T}_{sh}(i,j,{l}_{1},{l}_{2}),$$## 3.4.4.

#### Analysis and discussion

In addition to systematic noise, the scattered disturbance illuminance from the intermediary reflective surface is the main factor in image quality degradation.

1. For photoelectric imaging systems based on human eyes, detection and identification of the extended source target depend on the target contour brought out by the apparent contrast of the image. Although people simply believe that 3 to 5% is the contrast threshold for target detection with human eyes, the more effective approaches are based on minimum resolvable contrast (MRC)

^{8}^{,}^{9}or target task performance (TTP).^{10}Enhancing the apparent contrast of the image is conducive to improving the performance of non-line-of-sight imaging system.2. It is impossible to separate ${E}_{d}(t)$ from ${E}_{st}(t)$ because they are synchronized with each other. However, ${E}_{\text{back}}$ may be reduced or eliminated by time gating, as shown in Fig. 7.

• ${l}_{2}\ge c\xb7{t}_{p}/2$: ${E}_{\text{back}}$ is separable from ${E}_{st}(t)$ [as in Fig. 7(a)]. ${E}_{\text{back}}$ may be filtered out completely without reducing the laser illuminance efficiency, and thus, the highest contrast image is obtained.

• ${l}_{2}<c\xb7{t}_{p}/2$: When the gate opens upon the return of the target signal, part of the first-order backscattering is accepted by the detector, and the image contrast is decreased to some extent. To this end, the image contrast can be improved by adjusting ${t}_{d}$ [as in Fig. 7(b)].

${E}_{\text{back}}$ can be totally filtered out by setting the strobe gating delay time to ${t}_{d}=2{l}_{2}/c+{t}_{p}$. An optimum contrast image can be obtained if the gate width is simultaneously reduced to ${t}_{g}=2{l}_{2}/c$. However, utilization rate of the laser illuminance decreases with increase of ${t}_{d}$. The laser range-gated imaging system cannot receive the target signal illuminance with ${t}_{d}=2({l}_{1}+{l}_{2})/c+{t}_{p}$, and the efficiency decreases to 0. Since the laser range-gated imaging itself is low-signal-to-noise ratio (SNR) imaging, sacrificing the laser illuminance efficiency significantly reduces the SNR of the image. Therefore, practical application of non-line-of-sight imaging with laser range-gated imaging requires careful balance of the overall system.

## 4.

## Imaging Characteristic Analysis of Non-Line-of-Sight Imaging

Based on the image contrast model of non-line-of-sight imaging constructed above, this section analyzes the non-line-of-sight imaging results with different intermediary reflective surfaces. Typical parameters of non-line-of-sight imaging with a laser range-gated imaging system are shown in Table 1. The imaging depth of field is 36 m with a laser pulse width of 120 ns and a gate width of 120 ns; 40 laser pulses are received by the photodetector within a frame when the imaging frame rate is 25 fps and the laser pulse frequency is 1000 Hz, so that a laser illuminance equivalent to 4.8 *μ*J is received by the photodetector per frame.

## Table 1

Parameters of non-line-of-sight imaging with the laser range-gated imaging system.

Parameter | Value |
---|---|

ρt (target reflectivity) | 1 |

ρb (background reflectivity) | 0 |

Φ (energy of a single pulse in μJ) | 0.12 |

tp (laser pulse width in ns) | 120 |

Laser pulse frequency in Hz | 1,000 |

θ (divergence angle of the laser beam) | 2 deg |

D0 (diameter of objective lens in mm) | 20 |

f’ (focal length of objective lens in mm) | 80 |

τ0 (lens transmittance) | 0.9 |

Size of photodetector (mm2) | 14×10.5 |

tg (gate width in ns) | 120 |

Image frame rate (Hz) | 25 |

## 4.1.

### Impacts of Reflection Characteristics of the Intermediary Reflective Surface on Non-Line-of-Sight Imaging

The impacts of the reflection characteristics of the intermediary reflective surface on non-line-of-sight imaging are considered with ${l}_{1}=5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}$, ${l}_{2}=20\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}$, and $\alpha =45\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$. Since ${l}_{2}\ge c$·${t}_{p}/2=18\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}$, the target signal illuminance and first-order backscattering are separable.

Figure 8 shows the curves of the target signal illuminance ${E}_{st}$, the first-order backscattering ${E}_{\text{back}}$, and the third-order forward-scattering ${E}_{d}$ for different values of the diffuse reflectivity ${\rho}_{d}$ of the intermediary reflective surface with a CW laser. Since ${E}_{\text{back}}$ and ${E}_{d}$ are unevenly distributed over the imaging field, the ${E}_{\text{back}}$ and ${E}_{d}$ in Fig. 8 are the respective means of ${E}_{\text{back}}$ and ${E}_{d}$ in the center area of the imaging field. The curves in Fig. 8 indicate the following conclusions:

1. Under a CW laser, the target signal illuminance ${E}_{st}$ monotonically decreases with ${\rho}_{d}$, and first-order backscattering ${E}_{\text{back}}$ monotonically increases with ${\rho}_{d}$. ${E}_{st}$ is larger than ${E}_{\text{back}}$ only when ${\rho}_{d}$ is $<0.05$ (i.e., diffuse reflection takes up only 5% of the total reflection illuminance, and the specular reflection illuminance dominates), that is to say that non-line-of-sight imaging can be realized with only a CW laser when the intermediary reflective surface has very strong specular reflection characteristics. This conclusion indicates that practical realization of non-line-of-sight imaging must use a pulsed laser range-gated imaging mode in order to filter out ${E}_{\text{back}}$. In the instance indicated in Fig. 8, ${E}_{\text{back}}$ can be totally filtered out without reducing the laser illuminance efficiency.

2. The target signal illuminance ${E}_{st}$ spans more than three orders of magnitude from $4.6\times {10}^{-3}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{W}/{\mathrm{m}}^{2}$ for a completely specular reflection (${\rho}_{d}=0$) to $5.1\times {10}^{-7}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{W}/{\mathrm{m}}^{2}$ for a nearly completely diffuse reflection (${\rho}_{d}=0.99$). ${E}_{d}$ is several orders of magnitude smaller than ${E}_{st}$, and ${E}_{d}$ is two orders of magnitude smaller than ${E}_{st}$ even when ${\rho}_{d}$ is equal to 0.99. Figure 9 shows that the distribution of ${E}_{d}$ on the photodetector for ${l}_{1}=5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}$, ${l}_{2}=20\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}$, $\phi =0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$, $\theta =2\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$, and ${\rho}_{d}=0.5$ is within the same order of magnitude for $\alpha $ (the tilt angle of the intermediary reflective surface) equals 20 or 60 deg. The image quality is severely disturbed by quantum noise when the target signal illuminance is too low, so ${\rho}_{d}\le 0.7$ is selected as an effective non-line-of-sight imaging condition. ${E}_{st}$ with ${\rho}_{d}=0.7$ is one order of magnitude smaller than ${E}_{st}$ with ${\rho}_{d}=0$. ${E}_{d}$ is always two orders of magnitude smaller than ${E}_{st}$ when ${\rho}_{d}$ changes from 0 to 0.7. In view of the above discussion, ${E}_{d}$ is often negligible.

## 4.2.

### Impacts of the Location of the Intermediary Reflective Surface and the Gate Width on Non-Line-of-Sight Imaging

The parameter ${t}_{g}$ is chosen to be equal to ${t}_{p}$ to fully use the laser illuminance and to reduce the noise for a typical laser range-gated imaging system, in which case the imaging depth of field is $c$·${t}_{p}$. The gate opens when the echo signal reaches the detector, so that the gate delay time ${t}_{d}$ equals $2({l}_{1}+{l}_{2})/c$. The following discussion is based on the assumption that ${Q}_{\mathrm{eq}}$ is equal to 0.

Figure 10 shows a plot of the image contrast $C$ versus distance ${l}_{2}$ with ${l}_{1}=5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}$ and ${t}_{d}=2({l}_{1}+{l}_{2})/c$. The plot shows that $C$ first decreases with ${l}_{2}$ and then increases with ${l}_{2}$ to a certain constant value at ${l}_{2}=18\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}$ when ${E}_{\text{back}}$ is totally separated from ${E}_{st}$ and that finally $C$ is nearly fixed with increasing ${l}_{2}$. Figure 11 shows a plot of the image contrast $C$ versus distance ${t}_{d}$ with ${l}_{1}=5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}$ and ${l}_{2}=10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}$. The plot shows that $C$ increases with ${t}_{d}$ and that $C$ obtains its maximum value when ${E}_{\text{back}}$ is totally filtered out at ${t}_{d}=140\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{ns}$.

## 5.

## Conclusions

This paper describes the construction of an image contrast model of non-line-of-sight imaging with a laser range-gated imaging system on the basis of the modeling of an intermediary reflective surface. The image contrast model analyzes three factors affecting the image contrast of non-line-of-sight imaging with different intermediary reflective surfaces, where the three factors are target signal illuminance, first-order backscattering, and third-order forward-scattering. A non-line-of-sight imaging experiment simulation was carried out with the parameters of an actual laser range-gated imaging system. The simulation results show that (1) the first-order backscattering is the dominant disturbance to non-line-of-sight images and can be effectively filtered out by laser range-gated imaging and (2) the effect of the third-order forward-scattering on non-line-of-sight images is often negligible since the third-order forward-scattering is insignificant compared to the target signal illuminance in most cases. The image contrast model of non-line-of-sight imaging with a laser range-gated imaging system gives insights into the theoretical analysis and system design of non-line-of-sight imaging.

The simulation results also show that the reflection characteristics of the intermediary reflective surface have a significant effect on non-line-of-sight imaging. The bidirectional reflectance distribution function (BRDF) is a more comprehensive description of surface reflection characteristics; thus, BRDF modeling of typical construction material surfaces is a direction for further research. Meanwhile, since non-line-of-sight imaging with laser range-gated imaging system is usually low-SNR imaging, estimating the range of non-line-of-sight imaging using MRC or TTP is another direction for further research.

## Acknowledgments

This work was supported by the Specialized Research Fund for the Doctoral Program of Higher Education (No. 20111101110013) and the Foundation of the Science and Technology on Low-Light-Level Night Vision Laboratory (No. J20110504). The authors would like to thank all the reviewers for their helpful comments and constructive suggestions. They would also like to thank Mr. Ye Peng and Dr. Xuan Hongwen for their assistance with the programming.

## References

## Biography

**Kaida Xu** received his BS from Mechanical Engineering College, China, in 2001. He is currently pursuing his PhD in optical engineering at Beijing Institute of Technology, China. His technical interests include low-level-light and infrared imaging technology.

**Weiqi Jin** received his PhD in optical engineering from Beijing Institute of Technology of China in 1990. He is currently a professor with the Department of Optical Engineering in the School of Optoelectronics at the Beijing Institute of Technology. His major research interests include low-level-light imaging, thermal imaging, and image processing. He is a fellow of SPIE.

**Shenyou Zhao** received his BS in control technology and instrument from Changchun University of Science and Technology in 2011. He is currently pursuing his master’s degree in optical engineering at Beijing Institute of Technology, China. His main technical interest is optical imaging technique.

**Jing Liu** received her BS in optical engineering from Beijing Institute of Technology of China in 2009. She studied augmented reality in Universidad Politecnica de Madrid of Spain from March to September 2009. She is currently pursuing her PhD in optical engineering at Beijing Institute of Technology of China. Her technical interests include underwater imaging, laser range-gated imaging, polarization imaging, and imaging process.

**Hui Guo** received her PhD in optical engineering from Beijing Institute of Technology in 2005. She is currently a doctoral supervisor in the National Key Laboratory of Science and Technology on Low Light Level Night Vision, Xi’an, China. Her major research direction is low-light-level night vision technology.

**Su Qiu** received his PhD from School of Optoelectronics at Beijing Institute of Technology in 2013. He is currently an assistant professor in School of Optoelectronics at Beijing Institute of Technology. His major research interests include digital image processing, low-level-light and infrared imaging technology.