Three-Dimensional Imaging of Vortex Electromagnetic Wave Radar with Integer and Fractional Order OAM Modes

: Vortex electromagnetic (EM) waves, with different orbital angular momentum (OAM) modes, have the ability to distinguish the azimuth of radar targets, and then the two-dimensional reconstruction of the targets can be achieved. However, the vortex EM wave imaging methods in published research have no ability to obtain the elevation of the targets, and thus, the three-dimensional spatial structure and richer feature information of the radar target cannot be obtained. Therefore, a three-dimensional imaging method of vortex EM waves with integer- and fractional-order OAM modes is proposed in this paper, which can realize a three-dimensional reconstruction of a radar target based on a uniform circular array (UCA) with two-step imaging. First, the vortex EM wave with integer- and fractional-order OAM modes is generated, and the echo model with different OAM mode types is established. Thereafter, the echo with integer order is processed to obtain the range-azimuth image by fast Fourier transform (FFT). Then, in order to realize the three-dimensional reconstruction, the echo with fractional order is processed by utilizing the butterﬂy operation and analyzing the characteristics of the fractional Bessel function. Moreover, the resolution and reconstruction precision of the azimuth and elevation are analyzed. Finally, the effectiveness of the proposed method is veriﬁed by simulation experiments.


Introduction
Since L. Allen first revealed the optical vortex phenomenon in 1992 [1], this physical characteristic has received wide attention. According to classical electrodynamic theory [2], the electromagnetic (EM) wave also has spin angular momentum (SAM) and orbital angular momentum (OAM). Compared with the optical vortex, the EM wave with these two angular momentums is called the vortex EM wave [3]. The OAM corresponds to the front phase of the EM wave, and the wavefront phase can be modulated by changing the OAM modes; moreover, the integer orders thus satisfy orthogonality. This modulation capability of the wavefront phase provides a new degree of freedom in addition to the time, frequency and polarization domains. It has attracted attention in the fields of communication and radar, and a large number of research results have emerged [4][5][6][7][8][9][10][11][12], especially in radar imaging technology.
The azimuthal reconstruction of the target depends on the length of the synthetic or real aperture with a planar wave. It has limitations in the imaging time and imaging quality. Different from the planar wave imaging technology, the azimuthal reconstruction with a vortex EM wave no longer depends on the aperture's size but on the range of OAM modes. It means that the radar with a vortex EM wave can meet the needs of specific types of scatterers, and the resolution of the azimuth and elevation are also discussed. In Section 4, numerous simulations are carried out. Finally, in Section 5, conclusions are drawn.

Echo Model Based on UCA
The geometry of vortex EM wave detection based on a UCA is illustrated in Figure 1. The vortex EM waves with different orders are generated by a UCA, which comprises N antennas with the radius a, and the center of the UCA is set at the coordinate origin O. The nth element of the UCA is denoted by S n , and the transmitting signal from S n is denoted by s n (t). As known from Refs. [2,3], S n is modulated with the phase φ n and φ n = 2π N n; therefore, the transmitting signal s n (t) can be expressed as s n (t) = rect[t/T] • exp(jπKt 2 + j2π f c t) • exp(jφ n α) = s(t) exp(jφ n α) (1) where rect [•] is the function of the pulse signal, s(t) is a linear frequency modulation (LFM) signal and s(t) = rect[t/T] • exp(jπKt 2 + j2π f c t). K is the LFM rate, T is the transmitted pulse width, f c is the center frequency of the signal, and α is the order of the OAM mode. Thus, the total signal transmitted by the UCA can be expressed below.
s T (t) = N ∑ n=1 s(t) exp(jφ n α) (2) the presented method. The structure of this paper is organized as follows. In Section 2, echo models with integer and fractional order are built. In Section 3, the imaging flow is given for differen types of scatterers, and the resolution of the azimuth and elevation are also discussed. In Section 4, numerous simulations are carried out. Finally, in Section 5, conclusions are drawn.

Echo Model Based on UCA
The geometry of vortex EM wave detection based on a UCA is illustrated in Figure  1. The vortex EM waves with different orders are generated by a UCA, which comprise N antennas with the radius a, and the center of the UCA is set at the coordinate origin O . The nth element of the UCA is denoted by n S , and the transmitting signal from n S is denoted by ( ) n s t . As known from Refs. [2,3], n S is modulated with the phase n φ and . K is the LFM rate, T is the transmitted pulse width, c f is the center frequency of the signal, and α is the order o the OAM mode. Thus, the total signal transmitted by the UCA can be expressed below. Assume ( , , ) p p p P r θ ϕ is an arbitrary scatterer in space, where p r is the range be tween point P and the coordinate origin O , and p θ and p ϕ represent the elevation Assume P(r p , θ p , ϕ p ) is an arbitrary scatterer in space, where r p is the range between point P and the coordinate origin O, and θ p and ϕ p represent the elevation and azimuth of P, respectively. Based on the geometry, as shown in Figure 1, the echo of the scatterer P received by the UCA can be expressed as [9]: where τ p is the time delay and τ p ≈ 2r p c ; σ is the normalized scattering coefficient; c is the speed of light; k is the wavenumber; J α (•) is the Bessel function of the first kind of order α; and α is an integer number. According to Refs. [19,20], the signal transmitted by elements of the UCA can be phase-modulated and also amplitude-modulated; then, the vortex EM wave with fractional order is obtained. Similarly, the echo of scatterer P for fractional order can be given by where the value of α is the fraction in Equation (4). Assuming the target comprises M independent scatterers, the echo of the target can be uniformly expressed as below from the description in Equations (3) and (4) Assuming the range of the integer-order OAM mode is [−l, l], and the fractional order OAM modes separately are α 1/2 = 1 2 and α −1/2 = − 1 2 . Then, the echo signals with different kinds of orders in matrix form can be described as where [α 1 , α 1 , · · · α 2l ] corresponds to the range of the OAM-mode order [−l, l].

Three-Dimensional Imaging Algorithm
According to Equation (6), the echo signal in S can be expressed as a matrix, and the imaging of the range and azimuth can be obtained from S by the two-dimensional FFT method. However, it has no ability to obtain the elevation of the target. In order to obtain the elevation resolution and realize the three-dimensional reconstruction of the target, the signals with integer order in S are operated by two-dimensional FFT first. The information on the range and azimuth of the target can be obtained, which is called coarse imaging. Secondly, the scatterers are classified into two types according to the range distribution from the results of coarse imaging. Finally, the elevation of two kinds of scatterers can be solved by utilizing the echo with ± 1 2 OAM orders. As a result, the three-dimensional reconstruction is realized.

Coarse Imaging
The echo signals with integer order are selected to form a new matrix S in , and S in is defined as S in = [s M (t, α 1 ), . . . . ., s M (t, α 2l )] T . The center point P c of the target is initially set as the reference point, and then the echo of P c with integer order can be expressed as where s P c (t, α n ) is the echo of P c with nth order. The range of the target can be obtained through range imaging with matched filtering for S in , and the processed signal can be expressed as where s Match (t, α n ) can be expressed as where ps f t (•) is the point-spread function in the time domain. Then, S Match can be processed in terms of the OAM mode by FFT to obtain the azimuth image. With the above processing, the results of coarse imaging can be obtained. The scatterers are classified into two types by the distribution of scatterers in the range bin; that is, if there is only one scatterer in the range bin, the scatterer is defined as the first kind scatterer and denoted by P ; if more than one scatterer exists in the same range bin, these scatterers are defined as the second kind scatterers denoted by P . Assume the number of P is L 1 and the number of P is L 2 , then, L 1 and L 2 must satisfy L 1 + L 2 = M.

Three-Dimensional Imaging
Through coarse imaging, only two-dimensional information can be obtained. In order to achieve the elevation reconstruction, the signal matrix S f r with ± 1 2 orders is established, which is defined as S f r = [s M (t, α 1/2 ), s M (t, α −1/2 )] T . According to Refs. [40,41], when α equals ± 1 2 , the first type of Bessel function can be expressed as After the operation of matched filtering, the fractional order echo S f r can be further expressed as The transmitting signal is set as the reference signal. De-chirping and FFT are used to process S f r ; therefore, the processed S f r can be expressed as where R ∆ = r p − r re f and r re f is the range between the reference point and transmitters. Due to the transmitting signal being set as the reference signal, therefore, r re f = 0 and R ∆ = r p . From the properties of the sinc function, the information of the target range can be obtained by the operation with de-chirping and FFT successively. If the scatterers are in different range bins, each peak can represent only one scatterer. According to Equations (14) and (15), the amplitude of each scatterer on a one-dimensional range profile with ± 1 2 order can be expressed as Divide the amplitude A P,α −1/2 by A P,α 1/2 to yield the following: Equation (17) can be further expressed as The target elevation information corresponding to different range bins can be obtained through Equation (19). Although the two-dimensional imaging result of the target is obtained by coarse imaging, the target azimuth still needs to be solved from fine imaging by utilizing the echo signal with ± 1 2 OAM modes. The signals expressed in Equations (14) and (15) are processed with the butterfly operation, which is shown in Figure 2.
Divide the amplitude The target elevation information corresponding to d tained through Equation (19). Although the two-dimension is obtained by coarse imaging, the target azimuth still need ing by utilizing the echo signal with 1 2 ± OAM modes. T tions (14) and (15) The butterfly operation is based on fractional-order OAM modes.
In Figure 2, the sum of The signal S add can be further expressed as Assume s p = exp(jϕ p ) sin 2 (ka sin θ p ) + exp(−jϕ p ) cos 2 (ka sin θ p ), and exponential term exp(jϕ p ) can be expanded by utilizing the Euler formula. Then, s p can be further expressed as According to Equation (20), Equation (19) can be further represented as From Equation (21), if only the amplitude of the range image is considered, then the amplitude in different range bins can be expressed as Similarly, S sub can be expressed as The amplitude of signal S sub in different distance units can be expressed as Then, we have The elevation has been solved by Equation (18); thus, the azimuth of the target can be derived from Equation (25) as Therefore, according to Equations (18) and (26), the azimuth and elevation of the target can be obtained by processing the echo signal with fractional OAM modes. However, the above signal processing flow is based on the premise that the scatterers are distributed in different range bins. It means that the above three-dimensional imaging processing is applicable for the first kind of scatterers. If the scatterers belong to the second kind, there will be no resolution of the azimuth and elevation with the above processing. In order to obtain the elevation of the second kind of scatterers, the azimuth ϕ p from the result of the coarse imaging needs to be used. Equation (27) can be further described as Equation (28) shows that the elevation of different scatterers in the same range bin can be solved by using the information from coarse imaging, and then the second kind of scatterers can be three-dimensional reconstructed. Therefore, the three-dimensional reconstruction of the target can be realized from Equations (19), (27) and (28). It must be mentioned that the amplitude of the range imaging may be easily disturbed and unstable by many factors, such as noise. In fact, the range information can be obtained from just one pulse, and the radar transmits many pulses per second; thus, many groups of range images can be obtained per second. The average amplitude of range images can be used to solve the elevation and azimuth in order to reduce the effect of noise.
To sum up, the flow of the three-dimensional imaging method based on the vortex EM wave with integer-and fractional-order OAM modes are shown in Figure 3, and the steps can be described as follows: Step (6): by using the azimuth information of coarse imaging, the elevation of the second kind of scatterer is solved through Equation (28), and the three-dimensional reconstruction of this kind of scatterer is realized.
Step (7): the results in Steps 5 and 6 are fused to realize the three-dimensional reconstruction of the radar target.  According to the above steps, the three-dimensional reconstruction of the radar target can be realized. Next, let us discuss the resolution of the target range, azimuth and elevation of the proposed method. The proposed method realizes range imaging by using FFT. Therefore, according to Ref. [26], the range resolution depends on the point-spread function (PSF) and is inversely proportional to the signal bandwidth B , which can be expressed as Due to the different kinds of scatterers, there are also differences in solving the azimuth and elevation of the scatterers. For the second kind of scatterers, the azimuth imaging is achieved by FFT in terms of the OAM mode. Then, the resolution can be expressed as Step (1): the matrix of the echo signal S is established based on a UCA with integer and fractional orders.
Step (2): the echo signal S in with integer orders is processed by FFT, and coarse imaging, in terms of the range and azimuth, can be obtained.
Step (3): according to the result of coarse imaging, the scatterers are classified as the first or second kind of scatterers.
Step (4): the echo signal S f r with fractional orders is processed with the butterfly operation, then S add and S sub can be obtained.
Step (5): the elevation and azimuth of the first kind of scatterers are solved through Equations (18) and (26), and the three-dimensional reconstruction of this kind of scatterer is realized.
Step (6): by using the azimuth information of coarse imaging, the elevation of the second kind of scatterer is solved through Equation (28), and the three-dimensional reconstruction of this kind of scatterer is realized.
Step (7): the results in Steps 5 and 6 are fused to realize the three-dimensional reconstruction of the radar target.
According to the above steps, the three-dimensional reconstruction of the radar target can be realized. Next, let us discuss the resolution of the target range, azimuth and elevation of the proposed method. The proposed method realizes range imaging by using FFT. Therefore, according to Ref. [26], the range resolution depends on the point-spread function (PSF) and is inversely proportional to the signal bandwidth B, which can be expressed as Due to the different kinds of scatterers, there are also differences in solving the azimuth and elevation of the scatterers. For the second kind of scatterers, the azimuth imaging is achieved by FFT in terms of the OAM mode. Then, the resolution can be expressed as where ∆α is the range of the OAM mode order. The elevations of the second kind of scatterers are obtained by Equation (27), which is solved through the arcsine function.
Assuming the resolution of elevation is ρ θ , the calculation errors of elevation ∆x satisfy: According to the properties of trigonometric functions, the value of ∆x is different in the areas of [0, π 2 ], which is shown in Figure 4. Figure 4 shows the relationship of the calculation errors and ∆θ, and ∆θ is the elevation resolution by the reconstruction algorithm in this paper. It is obvious that ∆x decreases with increasing elevation θ. π , which is shown in Figure 4. Figure 4 shows the relationshi calculation errors and θ Δ , and θ Δ is the elevation resolution by the reconstructio rithm in this paper. It is obvious that x Δ decreases with increasing elevation θ . For the first kind of scatterers, the elevation is solved from Equation (18) by u the properties of the Bessel function with fractional order. By comparing Equat with Equation (28), the elevation of the first kind of scatterer is also obtained by an function, and thus it must satisfy Equation (31). However, the azimuth still need solved by Equation (27), which is obtained by an arctan function. Assuming the a resolution is ϕ ρ , the calculation error of elevation y Δ will satisfy π . Figure 5 shows the relationship between the calcula rors and ϕ Δ , and ϕ Δ is the azimuth resolution by the reconstruction algorithm paper. It is obvious that y Δ increases with the increasing azimuth ϕ .
x For the first kind of scatterers, the elevation is solved from Equation (18) by utilizing the properties of the Bessel function with fractional order. By comparing Equation (18) with Equation (28), the elevation of the first kind of scatterer is also obtained by an arcsine function, and thus it must satisfy Equation (31). However, the azimuth still needs to be solved by Equation (27), which is obtained by an arctan function. Assuming the azimuth resolution is ρ ϕ , the calculation error of elevation ∆y will satisfy According to the properties of trigonometric functions, the value of ∆y is also different in the areas of [0, π 2 ]. Figure 5 shows the relationship between the calculation errors and ∆ϕ, and ∆ϕ is the azimuth resolution by the reconstruction algorithm in this paper. It is obvious that ∆y increases with the increasing azimuth ϕ.
The requirement of ∆x and ∆y are compared with the same resolution, which is shown in Figure 6, and ∆ is defined as a calculation error of the azimuth or elevation. It can be found that the requirement of ∆y is significantly higher than the requirement of ∆x.
According to the imaging flow for the first kind of scatterers, the azimuth can be obtained not only from the coarse imaging but also from Equation (27). The resolution of azimuth depends on different factors in these two methods. In the coarse imaging method, the resolution satisfies Equation (30) and depends on the range of OAM modes; while in the second method, it satisfies Equation (32) and depends on the calculation errors. In certain conditions, especially when the azimuth angle is relatively large, the azimuth resolution in the coarse imaging method will be lower than the second method. The requirement of x Δ and y Δ are compared with the same resolution, w shown in Figure 6, and Δ is defined as a calculation error of the azimuth or eleva can be found that the requirement of y Δ is significantly higher than the requirem x Δ . According to the imaging flow for the first kind of scatterers, the azimuth can tained not only from the coarse imaging but also from Equation (27). The resolu azimuth depends on different factors in these two methods. In the coarse imaging m the resolution satisfies Equation (30) and depends on the range of OAM modes; w the second method, it satisfies Equation (32) and depends on the calculation er certain conditions, especially when the azimuth angle is relatively large, the azimu olution in the coarse imaging method will be lower than the second method.
The accuracy of target reconstruction is not only related to the resolution of el  The requirement of x Δ and y Δ are compared with the same resolution, w shown in Figure 6, and Δ is defined as a calculation error of the azimuth or elev can be found that the requirement of y Δ is significantly higher than the require x Δ . According to the imaging flow for the first kind of scatterers, the azimuth can tained not only from the coarse imaging but also from Equation (27). The resolu azimuth depends on different factors in these two methods. In the coarse imaging m the resolution satisfies Equation (30) and depends on the range of OAM modes; w the second method, it satisfies Equation (32) and depends on the calculation er certain conditions, especially when the azimuth angle is relatively large, the azim olution in the coarse imaging method will be lower than the second method.
The accuracy of target reconstruction is not only related to the resolution of el and azimuth but also related to the signal-to-noise ratio (SNR). The elevation and a of the target are obtained by calculating the amplitude term of the one-dimensiona profile of the echo signal; therefore, the echo signal with fractional order can be ex as: The accuracy of target reconstruction is not only related to the resolution of elevation and azimuth but also related to the signal-to-noise ratio (SNR). The elevation and azimuth of the target are obtained by calculating the amplitude term of the one-dimensional range profile of the echo signal; therefore, the echo signal with fractional order can be expressed as: where n(t) is Gaussian white noise. According to the properties of Gaussian white noise, Gaussian white noise has values in all frequencies, and thus, the amplitude term can be expressed in Equations (34) and (35) after de-chirping and the FFT operation.
A P,α −1/2 = A P,α −1/2 + N P,α −1/2 ( f ) (36) where N P,α 1/2 ( f ) and N P,α −1/2 ( f ) are the amplitudes of the Gaussian white noise of each scatterer with ± 1 2 order; thus, the v p is expressed as: where v randn is a random value that satisfies the Gaussian distribution, and it is known that the SNR can be approximated as: It can be found that the value of v randn will increase as the SNR decreases from Equation (38). It also leads to a decrease in the algorithm's reconstruction accuracy as the SNR decreases.

Simulation
In this section, the performance of the proposed method is validated by simulations. The UCA comprises 50 elements with a radius of 0.1 m, and the center of the UCA is set at the origin. The integer OAM mode during this is [−20, 20], and the fractional OAM modes are 1 2 and − 1 2 , respectively. The LFM signal is set as a transmitting signal; the center of frequency is 10 GHz; and the bandwidth of the signal is 100 MHz. There are seven independent scatterers on the target, and the coordinates of these scatterers are listed in Table 1 and shown in Figure 7.

Processing of Coarse Imaging
The echo signal in S is processed by FFT in terms of the frequency and OAM mode, respectively, and the results of the range and azimuth are shown in Figure 8. According to the conversion formula of the Cartesian coordinate system to a Spherical coordinate system, the scatterers in Table 1 can be written in the Spherical coordinate system in Table  2.

Processing of Coarse Imaging
The echo signal S in is processed by FFT in terms of the frequency and OAM mode, respectively, and the results of the range and azimuth are shown in Figure 8. According to the conversion formula of the Cartesian coordinate system to a Spherical coordinate system, the scatterers in Table 1 can be written in the Spherical coordinate system in Table 2.

Processing of Signal with Fractional OAM Mode
The scatterers are partitioned into two kinds after coarse imaging. In order to obtain the elevation information, the echo signal with 1 2 ± OAM modes is processed. The dechirping and the butterfly operation are used successively, and then the results of add S and sub S are shown in Figure 9, where Figure 9a,b show the range imaging of add S and sub S , respectively. It can be found that the distribution of the scatterers in Figure 9a is the same as that in Figure 8a; however, the amplitudes of the range image peaks are different. By utilizing this difference and the characteristic of the Bessel function, the elevation of the scatterers can be obtained.   Figure 8 shows the range and azimuth imaging of scatterers. It means the distribution of the scatterers can be obtained by coarse imaging. According to the results of coarse imaging, the scatterers P 6 and P 7 are located in the same range bin, and the scatterers P 1~P5 distribute in different range bins, respectively. According to the classification criteria for the scatterers defined in this article, the scatterers P 1~P5 are classified as the first kind of scatterers and P 6 and P 7 are classified as the second kind of scatterers.

Processing of Signal with Fractional OAM Mode
The scatterers are partitioned into two kinds after coarse imaging. In order to obtain the elevation information, the echo signal with ± 1 2 OAM modes is processed. The dechirping and the butterfly operation are used successively, and then the results of S add and S sub are shown in Figure 9, where Figure 9a,b show the range imaging of S add and S sub , respectively. It can be found that the distribution of the scatterers in Figure 9a is the same as that in Figure 8a; however, the amplitudes of the range image peaks are different. By utilizing this difference and the characteristic of the Bessel function, the elevation of the scatterers can be obtained.

Reconstruction of the First Kind of Scatterers
The elevation and azimuth of the first kind of scatterers can be obtained from Equations (18) and (26). In order to improve the reconstruction accuracy and reduce the effect of noise, multi-pulses are used to obtain the average of the amplitude of range imaging. Assume that the signal-to-noise ratio (SNR) is 5 dB, the pulse repetition frequency (PRF) is 1000, and the average amplitude of 1000-times-range imaging is used. The results of two-dimensional imaging, in terms of the elevation and azimuth, are shown in Figure 10, where the red color represents the ideal model of scatterers, and the blue color represents the rebuilt scatterers. Figure 10a shows the imaging of all the scatterers, and Figure 10b shows the imaging result of the first kind of scatterers. From the results shown in Figure  10, the imaging method for the first kind of scatterers can realize the accurate reconstruction of 1 P~5 P . However, for the second kind of scatterers 6 P and 7 P , only one scatterer can be rebuilt, and the solved information is also not the real coordinate of the scatterers.

Reconstruction of the First Kind of Scatterers
The elevation and azimuth of the first kind of scatterers can be obtained from Equations (18) and (26). In order to improve the reconstruction accuracy and reduce the effect of noise, multi-pulses are used to obtain the average of the amplitude of range imaging. Assume that the signal-to-noise ratio (SNR) is 5 dB, the pulse repetition frequency (PRF) is 1000, and the average amplitude of 1000-times-range imaging is used. The results of two-dimensional imaging, in terms of the elevation and azimuth, are shown in Figure 10, where the red color represents the ideal model of scatterers, and the blue color represents the rebuilt scatterers. Figure 10a shows the imaging of all the scatterers, and Figure 10b shows the imaging result of the first kind of scatterers. From the results shown in Figure 10, the imaging method for the first kind of scatterers can realize the accurate reconstruction of P 1~P5 . However, for the second kind of scatterers P 6 and P 7 , only one scatterer can be rebuilt, and the solved information is also not the real coordinate of the scatterers.
The reconstruction method for the first kind of scatterers proposed in this paper utilizes the amplitude of the range imaging; therefore, it is necessary to consider the performance of the proposed method under the different SNRs. Assuming the range of SNR is from −10 dB to 15 dB, the average deviation and average variance of elevation and azimuth are calculated and recorded as shown in Figure 11, where the red color represents the error of elevation and the blue color represents the error of azimuth. When the SNR is higher than 5 dB, the reconstruction method has good performance, and the errors are stable. The average deviation and average variance of elevation are no more than 0.06π and 0.005π, respectively. The average deviation and average variance of azimuth are no more than 0.03π and 0.001π, respectively. Compared to the reconstruction error of elevation and azimuth under the same SNR, it can be found that the estimation accuracy of the azimuth is obviously better than the elevation.
where the red color represents the ideal model of scatterers, and the blue color represents the rebuilt scatterers. Figure 10a shows the imaging of all the scatterers, and Figure 10b shows the imaging result of the first kind of scatterers. From the results shown in Figure  10, the imaging method for the first kind of scatterers can realize the accurate reconstruction of 1 P~5 P . However, for the second kind of scatterers 6 P and 7 P , only one scatterer can be rebuilt, and the solved information is also not the real coordinate of the scatterers. The reconstruction method for the first kind of scatterers proposed in this paper utilizes the amplitude of the range imaging; therefore, it is necessary to consider the performance of the proposed method under the different SNRs. Assuming the range of SNR is from −10 dB to 15 dB, the average deviation and average variance of elevation and azimuth are calculated and recorded as shown in Figure 11, where the red color represents the error of elevation and the blue color represents the error of azimuth. When the SNR is higher than 5 dB, the reconstruction method has good performance, and the errors are stable. The average deviation and average variance of elevation are no more than 0.06π and 0.005π , respectively. The average deviation and average variance of azimuth are no more than 0.03π and 0.001π , respectively. Compared to the reconstruction error of elevation and azimuth under the same SNR, it can be found that the estimation accuracy of the azimuth is obviously better than the elevation. When the elevation and azimuth have been solved, the three-dimensional reconstruction can be realized. Figure 12 shows the three-dimensional imaging of 1 P~5 P when the SNR is 5 dB. The effectiveness of the proposed method for the first kind of scatterers is verified.
angle/ When the elevation and azimuth have been solved, the three-dimensional reconstruction can be realized. Figure 12 shows the three-dimensional imaging of P 1~P5 when the SNR is 5 dB. The effectiveness of the proposed method for the first kind of scatterers is verified.
tion and azimuth under different SNRs.
When the elevation and azimuth have been solved, the three-dimensional reconstruction can be realized. Figure 12 shows the three-dimensional imaging of 1 P~5 P when the SNR is 5 dB. The effectiveness of the proposed method for the first kind of scatterers is verified.

Reconstruction of the Second Kind of Scatterers
From the simulation in Section 4.3, it can be found that the method for the first kind of scatterers is not suitable for the second kind of scatterers. The azimuth of P 6 and P 7 , obtained from the coarse imaging result, should be fully utilized for the three-dimensional reconstruction; thus, the azimuth results are used to solve the elevation of P 6 and P 7 through Equation (27). Assuming the SNR is 5 dB, the two-dimensional imaging result in the elevation and azimuth domains are shown in Figure 13. P 6 and P 7 can be distinguished in the elevation domain in Figure 13.

Reconstruction of the Second Kind of Scatterers
From the simulation in Section 4.3, it can be found that the method for the fir of scatterers is not suitable for the second kind of scatterers. The azimuth of 6 P a obtained from the coarse imaging result, should be fully utilized for the three-dime reconstruction; thus, the azimuth results are used to solve the elevation of 6 P through Equation (27). Assuming the SNR is 5 dB, the two-dimensional imaging r the elevation and azimuth domains are shown in Figure 13. 6 P and 7 P can be guished in the elevation domain in Figure 13. The reconstruction method for the second kind of scatterers also utilizes the tude of range imaging result; thus, the effect of the SNR on the performance of the m must be considered. Assuming the range of the SNR is from 0 dB to 15 dB, the a deviation, average variance and average mean square error (MSE) of elevation ar lated and recorded, as shown in Figure 14. When the SNR is higher than 5 dB, the struction method has good performance, and the errors are stable. The average de and average variance of elevation are no more than 0.06π and 0.005π , respective average MSE is no more than 0.002 π .  The reconstruction method for the second kind of scatterers also utilizes the amplitude of range imaging result; thus, the effect of the SNR on the performance of the method must be considered. Assuming the range of the SNR is from 0 dB to 15 dB, the average deviation, average variance and average mean square error (MSE) of elevation are calculated and recorded, as shown in Figure 14. When the SNR is higher than 5 dB, the reconstruction method has good performance, and the errors are stable. The average deviation and average variance of elevation are no more than 0.06π and 0.005π, respectively. The average MSE is no more than 0.002π.

Three-Dimensional Imaging of the Target
According to the imaging results in Sections 4.3 and 4.4, the two kinds of scatterers can be combined to recover the three-dimensional imaging result of the target, which is shown in Figure 15. The method proposed in this paper based on integer and fractional orders has good performance in the reconstruction of the target, and the three-dimensional coordinates of scatterers 1 P~7 P are accurately reconstructed.

Three-Dimensional Imaging of the Target
According to the imaging results in Sections 4.3 and 4.4, the two kinds of scatterers can be combined to recover the three-dimensional imaging result of the target, which is shown in Figure 15. The method proposed in this paper based on integer and fractional orders has good performance in the reconstruction of the target, and the three-dimensional coordinates of scatterers P 1~P7 are accurately reconstructed. In order to further validate the three-dimensional imaging algorithm proposed in this paper, the three-dimensional imaging algorithm proposed is used to image the simplified scattering model of the aircraft, and the three-dimensional imaging result is shown in Fig  ure 16. It can be seen that the aircraft scattering point model is well reconstructed from Figure 16.

Conclusions
This paper studied the three-dimensional imaging method based on the vortex EM wave with integer-and fractional-order OAM modes. Three-dimensional imaging is real ized through a two-step imaging process. Moreover, the resolution of the elevation and azimuth and the effects of the SNR are analyzed, respectively. The simulation results show that the proposed method can effectively reconstruct the three-dimensional image of the target. The biggest advantage of this method is that the three-dimensional imaging no longer depends on the aperture by relative motion between the radar and target or the antenna array itself. This means that the imaging method proposed in this article can mee the real-time imaging requirements for moving targets, such as security, low-speed aero space target detection and so on. This method provides a new technical approach for three-dimensional imaging based on vortex EM waves. In order to further validate the three-dimensional imaging algorithm proposed in this paper, the three-dimensional imaging algorithm proposed is used to image the simplified scattering model of the aircraft, and the three-dimensional imaging result is shown in Figure 16. It can be seen that the aircraft scattering point model is well reconstructed from Figure 16. In order to further validate the three-dimensional imaging algorithm proposed in this paper, the three-dimensional imaging algorithm proposed is used to image the simplified scattering model of the aircraft, and the three-dimensional imaging result is shown in Figure 16. It can be seen that the aircraft scattering point model is well reconstructed from Figure 16.

Conclusions
This paper studied the three-dimensional imaging method based on the vortex EM wave with integer-and fractional-order OAM modes. Three-dimensional imaging is realized through a two-step imaging process. Moreover, the resolution of the elevation and azimuth and the effects of the SNR are analyzed, respectively. The simulation results show that the proposed method can effectively reconstruct the three-dimensional image of the target. The biggest advantage of this method is that the three-dimensional imaging no longer depends on the aperture by relative motion between the radar and target or the antenna array itself. This means that the imaging method proposed in this article can meet the real-time imaging requirements for moving targets, such as security, low-speed aerospace target detection and so on. This method provides a new technical approach for three-dimensional imaging based on vortex EM waves.

Conclusions
This paper studied the three-dimensional imaging method based on the vortex EM wave with integer-and fractional-order OAM modes. Three-dimensional imaging is realized through a two-step imaging process. Moreover, the resolution of the elevation and azimuth and the effects of the SNR are analyzed, respectively. The simulation results show that the proposed method can effectively reconstruct the three-dimensional image of the target. The biggest advantage of this method is that the three-dimensional imaging no longer depends on the aperture by relative motion between the radar and target or the antenna array itself. This means that the imaging method proposed in this article can meet the real-time imaging requirements for moving targets, such as security, low-speed aerospace target detection and so on. This method provides a new technical approach for three-dimensional imaging based on vortex EM waves.
Author Contributions: J.L. and Y.C. proposed the method and designed the experiments; Y.L. and Q.Z. coordinated the reading, analyzing and categorizing of articles reviewed in this study; J.L. and Y.C. performed the experiments and wrote the paper; J.L. and Y.C. analyzed the data; X.L. and Y.L. revised the paper. All authors have read and agreed to the published version of the manuscript.