Three-Dimensional Imaging with Bistatic Vortex Electromagnetic Wave Radar

: The vortex electromagnetic wave radar based on monostatic geometry can realize the imaging of a radar target in both range dimension and azimuth dimension by adjusting vortex electromagnetic wave modes, but it has no ability to distinguish the elevation of the radar target. To achieve three-dimensional imaging of the radar target, a three-dimensional imaging method of vortex electromagnetic wave based on bistatic geometry is proposed. The echo model of the target based on vortex electromagnetic wave in bistatic geometry is constructed ﬁrstly. The elevation resolution of the target is realized, and then the three-dimensional coordinate of the target is reconstructed. Moreover, the elevation resolution and the computational complexity of the proposed method are also analyzed. Simulation results demonstrate the effectiveness of the proposed method. The results also illustrate that the array element requirements and the computational complexity of presented algorithm are less than Multiple Signal Classiﬁcation (MUSIC) algorithm.


Introduction
With the diverse applications for radar target detection, the traditional imaging technology cannot meet the new requirements, the imaging ability of synthetic aperture radar (SAR) or inverse synthetic aperture radar (ISAR), which is limited in the situation of no relative motion between the targets and radar. Driven by these actual requirements, new radar imaging technologies have emerged in past decades. In particular, the radar imaging technology based on vortex electromagnetic (EM) wave has attracted extensive attention. The vortex EM wave has the spiral structure along the propagation direction and the phase distribution is helical [1]. The different modes can be obtained by changing the orbital angular momentum (OAM); thus, the vortex EM wave provides a new freedom for radar imaging.
In the last ten years, a large number of research results have appeared regarding vortex EM wave imaging technology [2][3][4][5][6][7][8][9][10][11][12][13][14][15][16]. Ref. [2] shows the ability of vortex EM wave to distinguish the azimuth of radar targets by using the dual relationship between OAM modes and azimuth. Thereafter, Refs. [3][4][5][6][7][8][9] realized the high-resolution imaging of azimuth and reconstructed the two-dimensional imaging results of the target by using fast Fourier transform (FFT) or sparse Bayesian learning (SBL). Refs. [10,11] realized that radar coincidence imaging (RCI) depends on spatial phase distribution characteristics of the vortex EM wave in the different modes. The vortex EM wave radar with fractional modes has good anti-noise performance [12,13]. Ref. [14] reduces the requirement of the azimuth resolution on the amount of mode data by using compressed sensing (CS) theory. Ref. [15] only obtains elevation of the radar target by CS through the location relationship between the target and elements of uniform circular array (UCA). Ref. [16] uses Multiple Signal Classification (MUSIC) algorithm to further improve the azimuth resolution of radar target and obtain the imaging results in both azimuth dimension and elevation dimension through two-dimensional spatial search. However, it has no ability to achieve three-dimensional imaging of radar target because of the correlation problem of scattering points. In summary, the above research [2][3][4][5][6][7][8][9][10][11][12][13][14][15][16] reconstructs the imaging of radar targets with the monostatic geometry by using UCA. The transmitter is usually set at the center of the UCA, and the receiver uses the UCA in this observation geometry, thus only two-dimensional imaging can be obtained with this monostatic configuration.
In fact, it is valuable to obtain more characteristics of the target from three-dimensional (3-D) imaging results, such as the target size, the relative position of scatters and so on. In order to obtain these characteristics of the target, lots of research apply the vortex EM wave to the SAR or other conventional radars [17][18][19][20][21]. Ref. [22] based on interferometric SAR processing solves the height information with the designed observation geometry.
Ref. [23] utilizes the relative motion between radar and targets and obtains the 3-D imaging of the targets. Refs. [24][25][26] applies the vortex EM wave to SAR, the three-dimensional information is obtained in which the term of fast time, slow time and the OAM mode, respectively. It means that the cross range is solved by virtual antenna aperture and the azimuth is solved by the OAM modes. Refs. [23][24][25][26][27] describe the 3-D imaging of the target using monostatic geometry based on the vortex EM wave, but the 3-D imaging method still needs the aperture in the term of space.
In order to realize the three-dimensional imaging of stationary or non-relative moving targets based on vortex electromagnetic wave radar, a three-dimensional imaging method with bistatic vortex electromagnetic wave radar is proposed in this paper. Compared with monostatic geometry, the bistatic geometry is applied, where two transmitters are set outside of the UCA. The angle between transmitting path and receiving path is relatively large, which means a new observation perspective for the target. Thus, in addition to the range-azimuth imaging, the elevation of radar target is expected to be obtained by utilizing the geometric relationship between the transmitters and the receivers.
In this paper, the echo model of bistatic vortex electromagnetic wave radar is constructed, and the capability of the elevation resolution of radar target is analyzed in Section 2. In Section 3, the three-dimensional imaging method is given. The elevation resolution and the calculation complexity of proposed algorithm are discussed in detail. Compared with [9], the number of array elements of receiving unit and the calculation complexity are significantly reduced. Finally, the effectiveness of the proposed method is verified by simulation experiments in Section 4. Conclusions are made in Section 5.

Echo of Bistatic Vortex EM Wave Radar
As shown in Figure 1, the UCA is used as the receiving unit and the center of UCA is located at the coordinate origin O. The UCA is consisting of N omnidirectional elements placed equidistantly on the circle with radius b. The location of the nth elements of the UCA is denoted as S n . The transmitting unit is an independent omnidirectional antenna located at O, and the transmitted linear frequency modulation (LFM) signal is denoted as s(t). Assume P(r, θ, φ) is an arbitrary independent scatterer in the space. According to [3], the echo of the scatterer P received by the UCA can be expressed as: where σ is the scattering coefficient of P, τ is the delay of the signal and τ = 2R n c , c is the speed of the light, R n is the propagation distance in the space, and l is the mode number of the vortex EM wave. When the angle between the transmission path OP and the receiving path S n P is small, the radar geometry shown in Figure 1 can be approximated to a monostatic radar. Therefore, R n can be approximately calculated as 2r thus the time delay τ = 2r c Although the information of the range r and azimuth φ can be obtained from (1) by FFT method in terms of fast time and OAM mode, respectively, the elevation θ still Remote Sens. 2022, 14, 2972 3 of 13 cannot be calculated from the amplitude term of the echo. Therefore, in the monostatic geometry, only a two-dimensional image can be obtained. mode number of the vortex EM wave. When the angle between the transmission path OP and the receiving path n S P is small, the radar geometry shown in Figure 1 can be approximated to a monostatic radar. Therefore, n R can be approximately calculated as 2r , thus the time delay 2 = r c τ . Although the information of the range r and azimuth φ can be obtained from (1) by FFT method in terms of fast time and OAM mode, respectively, the elevation θ still cannot be calculated from the amplitude term of the echo. Therefore, in the monostatic geometry, only a two-dimensional image can be obtained. Let the transmitter move from point O to point O′ that located in the positive direction of X-axis, and the distance between O and O′ is d . The echo of scatterer P received by element n S can be expressed as: where n t′ is the fast time and = Compared with the monostatic geometry, the angle between the transmission path O P ′ and the receiving path n S P is relatively large, so this kind of radar geometry can be considered as bistatic. The propagation distance n R need be expressed as  Let the transmitter move from point O to point O that located in the positive direction of X-axis, and the distance between O and O is d. The echo of scatterer P received by element S n can be expressed as: where t n is the fast time and t n = t − τ n , τ n is the propagation delay and τ n = R n /c, T p is the signal pulse width, K is the rate of frequency modulation, ϕ n is the angle between the element S n and positive direction of X-axis in the XOY plane. The total echo received by the N elements of UCA can be expressed as: Compared with the monostatic geometry, the angle between the transmission path O P and the receiving path S n P is relatively large, so this kind of radar geometry can be considered as bistatic. The propagation distance R n need be expressed as R n = R 1 + R 2 , where R 1 is the distance from element O to point P and R 2 is the distance from point P to element S n Assume the unit vectorr is along the direction of the line OP, it can be written as:r = sin θ cos φ·X + sin θ sin φ·Ŷ + cos θ·Ẑ Then, the vector expression of line OS n can be given as r n = b cos ϕ n ·X + sin ϕ n ·Ŷ , and the vector expression of line OO can be given as r O = d·X. Thus, the distance R 1 and R 2 can be separately expressed as: According to (5) and (6), the R n with bistatic is different from monostatic by adding the term of d sin θ cos φ, thus the total echo can be further expressed as: , and according to the Janoci-Anger formula, s(n) can be expressed as: where J m (·) is the first kind of mth order Bessel function. It can be easily calculated that the expression of N ∑ n=1 e j(l−m)ϕ n equals N. When the value of N is large, the amplitude of lth order Bessel function is much larger than other orders. Thus, (8) can be approximately expressed as: Finally, the echo received by UCA of P can be expressed as: where σ = σNj l . Then, the 2-D Fourier transform can be used in processing echo data to obtain the range and azimuth of the target. Compare the expression of echo in (1) and (10), the information of azimuth is same in two different geometries, but the information of the range in bistatic geometry is expressed as r + d sin θ cos φ 2 . It means the elevation θ of target can be obtained from the range information.

Three-Dimensional Imaging Method in Bistatic Geometry
Suppose the radar target consists of M scatterers and the location of the mth scatterer is denoted by P m (r m , θ m , φ m ). In the bistatic geometry, the echo of radar target can be written as: where τ m is the time delay of the mth scatter and τ m = R m c , σ m is the scattering coefficient of the mth scatterer. The target can be imaged in range dimension by using matched filtering. The center point of the target is set as P 0 (r 0 , θ 0 , φ 0 ), then the echo of the point P 0 is set as the reference signal and denoted by s 0 (t, l). The signal output by the matched filtering can be expressed as: where σ = σ 2 m N 2 J l (kb sin θ m )J l (kb sin θ 0 ). The information of the range R m can be obtained from the peak locations of (12).
The FFT in terms of l is subsequently used to get the azimuth φ m . In fact, the Bessel function is a part of the amplitude of the echo, and its effects on imaging result of azimuth must be considered. As it is known that when kb sin θ 1, the Bessel function can be approximately expressed as: From (13), it can be found that, in the FFT result of the Bessel function, there will be a null at zero frequency and a pair of symmetrically distributed narrowband spectrum shift components at ± π 4 , which has no relation with the elevation θ. Therefore, the spectrum shift by J l (kb sin θ m )J l (kb sin θ 0 ) equals the effect by the square of Bessel function. The square of Bessel function is shown in Figure 2. It will exhibit symmetric spectral components at ± π 2 after FFT. According to [3], the frequency shift can be eliminated through Hilbert transformation, so the echo signal is firstly processed through Hilbert transformation before the FFT, then the azimuth information of the scatterer is obtained. π π π π × + + + From (13), it can be found that, in the FFT result of the Bessel function, there will be a null at zero frequency and a pair of symmetrically distributed narrowband spectrum shift components at 4 π ± , which has no relation with the elevation θ . Therefore, the spectrum shift by  Figure 2. It will exhibit symmetric spectral components at 2 π ± after FFT. According to [3], the frequency shift can be eliminated through Hilbert transformation, so the echo signal is firstly processed through Hilbert transformation before the FFT, then the azimuth information of the scatterer is obtained. can be given by: where r is the range vector and After the above processing, the azimuth φ m and the range R m have been obtained. However, R m is not the real range information of the target but includes the three-dimensional information of the target. From the expression R m = r m + d sin θ m cos φ m 2 , where d is known, the azimuth φ m can be solved. For M scatterers, the expression can be given by: where r is the range vector and r = [r 1 , r 2 , · · · r m ] T , θ is the elevation vector and θ = [sin θ 1 , sin θ 2 , · · · sin θ m ] T , A is the equation coefficient and A = d 2 [cos φ 1 , cos φ 2 · · · , cos φ m ] T , • is the Hadamard product. The elevation and the range of M scatterers need to be solved from (14) to achieve the target three-dimensional reconstruction. However, in fact, M groups of (r m , φ m ) cannot be solved just from (14).
Let an additional transmitting element be added at point O along the positive direction of X-axis, which is shown in Figure 3, and the distance between O and O is d , the signal transmitted by element O is s (t), which is orthogonal with s(t). It means the transmitters and receivers form two mutually independent channels, which is recorded as TR1 and TR2, respectively. Similarly, the range information R m can be obtained from the echo data of TR2, and the expression R m = r m + d sin θ m cos ϕ m 2 can also be written as: where A is the equation coefficient vector, A = d 2 [cos φ 1 , cos φ 2 , · · · , cos φ m ] T , R is the range vector of TR2, R = [R 1 , R 2 , · · · R m ] T . Equations (14) and (15) can be united to solve the elevation of the target, and the expression of elevation of the mth scatterer can be given by:  (14) and (15) can be united to solve the elevation of the target, and the expression of elevation of the m th scatterer can be given by: Figure 3. The radar geometry with two transmitters.
According to (14)- (16), the information of elevation is obtained. Thus, the elevation resolution is determined by the range resolution when the azimuth resolution is enough. Suppose the elevation resolution is no less than θ Δ , following must be satisfied: However, the value of x Δ is different in the various extent from 0 to / 2 π as shown in Figure 4. Figure 4 shows the value of x Δ is small when the elevation is relatively large, which means the bandwidth of signal should be set according to the maximum elevation of the target. According to (14)- (16), the information of elevation is obtained. Thus, the elevation resolution is determined by the range resolution when the azimuth resolution is enough. Suppose the elevation resolution is no less than ∆θ, following must be satisfied: where ρ R = c 4B and B is the bandwidth of the transmitting signal. The elevation resolution is in inverse proportion to the bandwidth of signal. If the elevation resolution is set as ∆θ, the bandwidth needs to be B ≥ c 4∆x , where ∆x = max[sin(θ m + ∆θ) − sin θ m ]. However, the value of ∆x is different in the various extent from 0 to π/2 as shown in Figure 4. Figure 4 shows the value of ∆x is small when the elevation is relatively large, which means the bandwidth of signal should be set according to the maximum elevation of the target.  is solved from (16), r can be obtained from (14) or (15). In summary, the 3-D imaging method proposed in this paper is described as follows: Step (1): a Cartesian coordinate system is constructed with the center of the UCA array as the origin, then two transmitters are set up at different positions along the X-axis, which composes TR1 and TR2 with the UCA, respectively; Step (2): the echo data of TR1 and TR2 is processed by matched filtering in terms of the fast time; Step (3): Hilbert transform and FFT are successively applied in terms of the OAM mode; Step (4): the solution of R and φ are united to solve the real range and elevation of each scatterer, then the 3-D imaging is obtained.
The flow of the three-dimensional imaging method is summarized in Figure 5. After θ= [sin θ 1 , sinθ 2 , · · · , sin θ m ] T is solved from (16), r can be obtained from (14) or (15). In summary, the 3-D imaging method proposed in this paper is described as follows: Step (1): a Cartesian coordinate system is constructed with the center of the UCA array as the origin, then two transmitters are set up at different positions along the X-axis, which composes TR1 and TR2 with the UCA, respectively; Step (2): the echo data of TR1 and TR2 is processed by matched filtering in terms of the fast time; Step (3): Hilbert transform and FFT are successively applied in terms of the OAM mode; Step (4): the solution of R and φ are united to solve the real range and elevation of each scatterer, then the 3-D imaging is obtained.
The flow of the three-dimensional imaging method is summarized in Figure 5.
Step (2): the echo data of TR1 and TR2 is processed by matched filtering in terms of the fast time; Step (3): Hilbert transform and FFT are successively applied in terms of the OAM mode; Step (4): the solution of R and φ are united to solve the real range and elevation of each scatterer, then the 3-D imaging is obtained.
The flow of the three-dimensional imaging method is summarized in Figure 5. The computational complexity of the proposed method is obviously less than the MUSIC method. The computational complexity of the standard MUSIC is given by [28]: The computational complexity of the proposed method is obviously less than the MUSIC method. The computational complexity of the standard MUSIC is given by [28]: where J is the number of spatial search points and J = π 2∆θ , ∆θ is the elevation resolution. The number of the target must be known in advance for MUSIC algorithm, and the number of the OAM mode needs to be more than the target scatterer number [16]. Thus, the number of the array elements N is assumed one more than the number of target M, then it yields: It is obvious that the computational complexity of MUSIC is mainly decided by the number of array elements and the elevation resolution. The computational complexity of the proposed method depends on 2-D FFT in terms of the fast time and the OAM mode, and it can be given by: (20) where N range is the number of sampling points in the fast time domain, 6N range log 2 (N range ) is the amount of calculation in range imaging and 2N log 2 (N) is the amount of calculation in azimuth imaging. When N and M are set as the same as the MUSIC method, (20) can be further expressed as: It can be found that the computational complexity of the proposed method is mainly decided by N and N range . N range is in direct proportion with the bandwidth of transmitting 8 of 13 signal such as N range = BX/c, where X is the range of detection area. The range resolution has the relation with the bandwidth of transmitting signal, so N range can be given by: The elevation resolution of the proposed method is decided by the range resolution, so N range will be increased with the rise of the elevation resolution and lead to higher computational complexity.
By comparing (19) and (21), the term related to N is the cubic term and the residual term is the square of N and J for the MUSIC method, but the term in (21) related to N and N range is all approximate linear relation. With the increased resolution in elevation dimension, the amount of calculation of MUSIC will increase exponentially, but will increase only linearly for the proposed method.

Simulation
Assume the UCA consists of 50 elements, the center of array is located at the origin and the radius of UCA is 0.3 m. The two transmitters are set at the direction of X-axis in different position and the distance between them is 2 m. The LFM signal is set as transmitting signal, the center carrier frequency is 10 GHz and the bandwidth is 2 GHz. The range of the OAM mode is [−20, 20]. Four scatterers are set in the space, which are P 1 (298, 0.1π, 0.1π), P 2 (300, 0.15π, 0.3π), P 3 (303, 0.3π, 0.15π) and P 4 (294, 0.35π, 0.3π) respectively.

Range Imaging in Bistatic Geometry
The result of the scatterers in the range dimension can be obtained through the matched filtering as shown in Figure 6. Figure 6a is the range imaging in the monostatic geometry based on [4][5][6], and the result is r= [298 300 303 294] T and it is same as the real range of target. Figure 6b,c show the range imaging in the monostatic geometry of TR1 and TR2, respectively. According to (10), the information of R m contains the 3-D coordinate information of target and R m = r m + d sin θ m cos ϕ m 2 ; thus, from the expression of the R m , the range result of TR1 and TR2 is listed in Table 1. It is verified by simulation results in Figure 6b,c, respectively. range of target. Figure 6b,c show the range imaging in the monostatic geometry of TR1 and TR2, respectively. According to (10), the information of m R contains the 3-D  Figure 6b,c, respectively.

2-D Imaging in Bistatic Geometry
After the Hilbert transform and 2-D FFT, the imaging result of the range-azimuth can be obtained in the bistatic geometry as shown in Figure 7. Figure 7a,b shows the result of TR1 and TR2, respectively.

2-D Imaging in Bistatic Geometry
After the Hilbert transform and 2-D FFT, the imaging result of the range-azimuth can be obtained in the bistatic geometry as shown in Figure 7. Figure 7a,b shows the result of TR1 and TR2, respectively. According to the 2-D imaging result of TR1 and TR2, the real range r and the elevation θ of the target can be solved by utilizing (14)- (17). The ideal imaging and the simulation imaging by the proposed method are shown in Figure 8. The red points represent the ideal imaging, and the blue points represent the simulation results, respectively. According to the 2-D imaging result of TR1 and TR2, the real range r and the elevation θ of the target can be solved by utilizing (14)- (17). The ideal imaging and the simulation imaging by the proposed method are shown in Figure 8. The red points represent the ideal imaging, and the blue points represent the simulation results, respectively. According to the 2-D imaging result of TR1 and TR2, the real range r and the elevation θ of the target can be solved by utilizing (14)- (17). The ideal imaging and the simulation imaging by the proposed method are shown in Figure 8. The red points represent the ideal imaging, and the blue points represent the simulation results, respectively. It can be seen that the scatterers 1 P , 2 P and 3 P can be reconstructed accurately, but there is some error between the simulation result and actual position of the scatterer 4 P .
The bandwidth of signal decides the range resolution. The range result of some scatterers may shift one range unit, thus the maximum error in range dimension by proposed method is 2-unit. Furthermore, taking scatter 3 P in Figure 8 for analysis, the range difference R Δ between two channels is 0.2668 m in theory, and it corresponds to 8-unit in range. However, the range difference in simulation is 0.2625 m and only corresponds to 7-unit, which causes there to be about 0.007π elevation difference between the simulation and the true value. Figure 9 shows the average calculation error is no more than 0.011π when the elevation resolution is high enough. It can be seen that the scatterers P 1 , P 2 and P 3 can be reconstructed accurately, but there is some error between the simulation result and actual position of the scatterer P 4 . The bandwidth of signal decides the range resolution. The range result of some scatterers may shift one range unit, thus the maximum error in range dimension by proposed method is 2-unit. Furthermore, taking scatter P 3 in Figure 8 for analysis, the range difference ∆R between two channels is 0.2668 m in theory, and it corresponds to 8-unit in range. However, the range difference in simulation is 0.2625 m and only corresponds to 7-unit, which causes there to be about 0.007π elevation difference between the simulation and the true value. Figure 9 shows the average calculation error is no more than 0.011π when the elevation resolution is high enough.  It is known that the elevation of P 4 , the bandwidth of signal and the range resolution is 0.35π, 2 GHz and 0.0375 m, respectively. According to (17), the elevation resolution of P 4 and range resolution ρ r have the relation shown in Figure 10. When ρ r is 0.0375 m and elevation is 0.35π, the elevation error is between [4 • 5 • ]. The elevation of P 4 is 59 • in simulation and the real elevation of P 4 is 63 • , the error is about 4 • , which is consistent with theoretical analysis. is 0.35π , 2 GHz and 0.0375 m, respectively. According to (17), the elevation resolution of 4 P and range resolution r ρ have the relation shown in Figure 10. When r ρ is 0.0375 m and elevation is 0.35π , the elevation error is between [4 5 ]   . The elevation of 4 P is 59  in simulation and the real elevation of 4 P is 63  , the error is about 4  , which is consistent with theoretical analysis.

3-D Imaging in Bistatic Geometry
The azimuth ϕ can be obtained through FFT in terms of the OAM mode, and the real range r and elevation θ can be solved by utilizing (14)- (17). The 3-D imaging result is shown in Figure 11, where the red points are the actual positions, and the blue points are the reconstructed positions. It can be seen that the proposed algorithm can accurately reconstruct the 3-D imaging results of targets.

3-D Imaging in Bistatic Geometry
The azimuth ϕ can be obtained through FFT in terms of the OAM mode, and the real range r and elevation θ can be solved by utilizing (14)- (17). The 3-D imaging result is shown in Figure 11, where the red points are the actual positions, and the blue points are the reconstructed positions. It can be seen that the proposed algorithm can accurately reconstruct the 3-D imaging results of targets. In the next section, the computational complexity of the proposed algorithm is considered. Ref. [16] utilize the MUSIC algorithm through the spatial two-dimensional spectral estimation to distinguish azimuth and elevation of the target. The twodimensional angular resolution of the algorithm depends on the range of the OAM mode when the signal-noise-ratio (SNR) and the number of snapshots are fixed. Thus, the calculation amount of MUSIC is also decided by the range of the OAM mode. In contrast, the computational complexity of the proposed method is not decided by the number of the mode but by the bandwidth of the transmitted signal. The computational complexity of the proposed method is obviously less than [16] with the same azimuth resolution. When the range and azimuth resolution are, respectively, set as 0.0375 m and 0.2π , the calculation amount of proposed method and MUSIC method with different elevation resolution is shown in Table 2.
From Table 2, it can be seen that the proposed method has superiority in the element number and calculation amount especially when the elevation resolution is relatively high. When 3 θ Δ =  , the requirement of element of UCA decline to a quarter and the calculation amount of proposed method is just about 1/8 of the MUSIC method.  In the next section, the computational complexity of the proposed algorithm is considered. Ref. [16] utilize the MUSIC algorithm through the spatial two-dimensional spectral estimation to distinguish azimuth and elevation of the target. The two-dimensional angular resolution of the algorithm depends on the range of the OAM mode when the signal-noiseratio (SNR) and the number of snapshots are fixed. Thus, the calculation amount of MUSIC is also decided by the range of the OAM mode. In contrast, the computational complexity of the proposed method is not decided by the number of the mode but by the bandwidth of the transmitted signal. The computational complexity of the proposed method is obviously less than [16] with the same azimuth resolution. When the range and azimuth resolution are, respectively, set as 0.0375 m and 0.2π, the calculation amount of proposed method and MUSIC method with different elevation resolution is shown in Table 2. From Table 2, it can be seen that the proposed method has superiority in the element number and calculation amount especially when the elevation resolution is relatively high. When ∆θ = 3 • , the requirement of element of UCA decline to a quarter and the calculation amount of proposed method is just about 1/8 of the MUSIC method.

Conclusions
In this paper, a novel method is presented for obtaining the 3-D imaging results with bistatic vortex electromagnetic wave radar. The echo in bistatic vortex electromagnetic wave radar is analyzed and the possibility of elevation resolution is shown. The 3-D imaging method is proposed in order to obtain the spatial information of the target in a spherical coordinate. The simulation results verify the effectiveness and the superiority of the proposed method in array number and computational complexity.