Next Article in Journal
A Multivariate Drought Index for Seasonal Agriculture Drought Classification in Semiarid Regions
Next Article in Special Issue
A Wideband Noise Radar System Using a Phased Array with True Time Delay
Previous Article in Journal
Automatic Detection of Pothole Distress in Asphalt Pavement Using Improved Convolutional Neural Networks
Previous Article in Special Issue
Calibration of MIMO Radar Transmitting and Receiving Array Using Scene Object Measurement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Ship Discrimination Method Based on High-Frequency Electromagnetic Theory

1
Department of Electronic Engineering, Tsinghua University, Beijing 100084, China
2
College of Missile Engineering, Rocket Force Engineering University, Xi’an 710025, China
3
School of Computer and Communication Engineering, University of Science and Technology Beijing, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(16), 3893; https://doi.org/10.3390/rs14163893
Submission received: 2 July 2022 / Revised: 4 August 2022 / Accepted: 9 August 2022 / Published: 11 August 2022
(This article belongs to the Special Issue Recent Advances in Signal Processing and Radar for Remote Sensing)

Abstract

:
Ship target detection using radar has important applications in the military and civilian fields. As a decoy, the corner reflector (CR) can successfully deceive a radar by its strong radar cross-section (RCS) to protect a ship. In order to discriminate between a CR and ship, this paper proposes a discrimination method based on three-dimensional characteristics. First, we deduce the basic scattering of CR by the high-frequency electromagnetic theory, and propose a CR decomposition which can solve the problem that the Krogager decomposition has terrible errors in clutter. Then, we introduce the definition of the main scattering polarization and give the multi-dimensional characteristic of CR. Subsequently, we analyze the spatial-time characteristic of radar based on the three-dimensional proportional guidance. With the CR mean square error (MSE), a CR discrimination method is proposed based on the time-spatial-polarization (TSP) joint domains. Finally, the proposed method is analyzed and compared using the fully polarimetric data of Feko software, which can achieve 95% discrimination probability and 4.1% false alarm probability.

1. Introduction

Ship target detection is essential for many applications in marine target monitoring [1], maritime rescue [2], and port security [3]. Especially in the military field, the ability to effectively detect a ship target whilst being interfered with by false targets is very important. The corner reflector (CR) can be used for polarimetric synthetic aperture radar (SAR) calibration [4], false targets [5], and polarimetric radar interferometer [6] due to many advantages [7]. In the application of a false target, CR can successfully deceive a radar by its strong radar cross-section (RCS) for protecting important targets, i.e., ships. In this paper, our goal was to discriminate the CR and ship for target detection. In the existing literature, there are usually three types of discrimination methods, i.e., high-range resolution profile (HRRP), micro-Doppler, and polarization, as shown in the following.
For the discrimination method based on HRRP, the difference between the scattering points between the ship and single CR was analyzed in [8], and then the support vector machine (SVM) was utilized for discrimination. Subsequently, to discriminate the array CRs, six parameters of HRRP were introduced, which could find that the CR HRRP was easily affected by the arrangement and spacing in [9]. A novel HRRP-based method by changing the modulation waveform of a transmitted wave was proposed in [10]. The discriminations of the dilute jamming and centroid jamming of CR were analyzed in [11], which ignored the relationship between the number of sub-strong scattering points and the detection threshold. After extracting multiple geometric characteristics of HRRP, many references also proposed different discrimination methods based on data mining. For example, the sparse dictionary learning was utilized to reconstruct the target, and then a discrimination method based on the error ratio of ships and CRs was proposed in [12]. Nine parameters of HRRP were introduced for training the deep neural network in [13]. The extreme learning machine algorithm was utilized in [14]. Although the difference between the ship and single CR is obvious, the array CR can realistically simulate a ship’s radial distribution through a reasonable arrangement, so that there is no significant difference between the array CR and ship in HRRP [15]. Additionally, the polarization angle [16], RCS [17], and HRRP [18] all have very similar characteristics to ships, by optimizing the arrangement, spacing, and type of array CRs. The above conclusions can effectively improve the similarity between the CR and ship’s HRRP, thus reducing the effectiveness of the HRRP-based methods.
The second type of discrimination is based on the micro-Doppler. Chen first introduced the concept of the micro-Doppler, which could be used in many different applications in [19,20], e.g., the vehicle discrimination method [21] and the false warhead discrimination method [22]. Subsequently, some references also proposed the micro-Doppler-based methods to discriminate CR and ship. For example, the micro-motion model of CR with the effect of ocean waves was established, and the difference in roll direction was utilized to discriminate in [23,24]. According to that, the mass and volume of CR were much smaller than that of the ship, and it could find that the CR micro-Doppler phenomenon was more obvious in [25]. Zhu et al. expanded the micro-Doppler characteristic with three parameters in [26]. In summary, the aforementioned references all analyzed the difference between a CR and ship in micro-Doppler. However, these assumed that the ocean wave only propagates in one fixed direction, and ignored the influence of wind speed or relative motion in [23]. Furthermore, the micro-motion model of a ship should include 6 degrees of freedom, but there are only two degrees of freedom analyzed to simplify the model in [23,24,25,26].
Last but not least, there are lots of discrimination methods based on the polarization information. The Krogager decomposition was used to decompose the polarization scattering matrix (PSM) of ships and CRs, and then the changes of basic scattering were compared in [27]. The Krogager decomposition was also utilized to obtain the coefficient of each basic scattering, and a discrimination method by SVM was proposed in [15]. Fang et al. optimized the odd scattering and even scattering matrices, and then utilized the polarimetric similarity parameter to discriminate CR in [28,29]. Wang et al. introduced five polarization invariants and one polarization shape factor of CR, and then proposed an SVM-based method to discriminate CR in [30,31]. Considering the measurement accuracy of the polarization phase in a sea clutter environment, two polarization invariants were ignored, and thus were not used as training parameters for SVM in [32]. Similarly, four characteristics of HRRP and two characteristics of polarization were analyzed in [33]. In brief, the aforementioned references based on polarization could be divided into two categories. (1) By extracting the polarization characteristics of the CR and ship, the discrimination methods are proposed by SVM. Such methods are not only susceptible to the data’s authenticity, but also to the arrangement of array CRs. (2) The coefficient of each basic scattering is obtained based on polarization decomposition, which could be utilized to discriminate the ship and CR. The theoretical basis of a secondary category is that CR is mainly odd scattering, whilst that of the ship is mainly even scattering [27,28]. However, the CR might mainly be even scattering and the ship might also have odd scattering, e.g., when the incident wave is far from the axis of symmetry. Therefore, the discrimination method would have the problem of high false alarm probability, which is based solely on the coefficient of each basic scattering.
The main goal of this paper was to discriminate the CR and ship. In contrast to the above three types of methods (i.e., HRRP, micro-Doppler, and polarization), we will analyze the scattering characteristics of the CR and ship based on high-frequency electromagnetic theory, which can obtain the characteristics of the time–spatial–polarization joint domains of the target. Regarding the analysis of the target characteristics using high-frequency electromagnetic theory, scholars have performed much pioneering work. Aiming at the CR formed by many surfaces on ships, Reference [34] developed a criterion that gave the required angle as a function of the desired RCS reduction and the electrical size of CR.
In [35], the PO and physical theory of diffraction (PTD) were used to determine the RCS of dihedral CRs in the azimuthal plane. In contrast to the method of PTD in [35], the uniform theory of diffraction (UTD) plus an imposed edge diffraction extension was used to predict the RCS of dihedral CRs in [36]. In comparison to shooting-bouncing-ray (SBR) with measurements and other rigorous theoretical methods, the authors in [37] verified that SBR yielded good results for calculating RCS from a three-dimensional target.
In order to analyze the PSM of a target, the authors in [38] proposed the analytic physical optics method based on PO-geometrical optics (GO). Since numerical methods cannot be applied to some complex structures, an asymptotic method based on the PO combined with the GO approximations was proposed for octahedral reflectors and icosahedral CRs in [39]. In [40], an analytic scattering model for 3D bistatic scattering was derived from a dihedral CR using GO and PO. Compared with the aforementioned high-frequency electromagnetic methods, we proposed a joint method of PO-GO-SBR in [41], which could derive a closed-form electromagnetic formula for a target with any structure. Based on the high-frequency electromagnetic theory, the present method in this paper independently analyzes the single reflection (SR), double reflection (DR), and triple reflection (TR) of PSM by controlling the number of scattering in the SBR, which is important in the wake of obtaining the odd and even scattering characteristics of target. In order to discriminate between the CR and ship, the polarization characteristic and the spatial–time characteristic of a target were developed. The main innovations of this paper are as follows.
  • The CR decomposition based on the main polarization. For the basic scattering of CR, the Krogager decomposition [27,28,29] has three assumptions, as presented in Section 3, which would increase the coefficient errors of odd scattering and even scattering as SCR decreases. To cure the above problem, this paper deduces the SR, DR, and TR of CR with the electromagnetic theory, and proposes a novel CR decomposition based on the main polarization, which could accurately obtain the coefficients of odd scattering and even scattering in clutter.
  • The multi-dimensional characteristic of CR. Aiming at the one-dimensional characteristic of CR and the ship, the characteristics of HRRP [8,9,10,11,12,13,14], micro-Doppler [23,24,25,26], polarization [27,28,29,30,31,32] are separately analyzed. For the multi-dimensional characteristic, only the authors in [33] comprehensively analyzed HRRP and polarization. However, its multi-dimensional characteristic is only a superposition of characteristic quantities, and does not consider the interaction between characteristics. Therefore, this paper introduces the definition of the main scattering polarization (MSP), and gives the multi-dimensional characteristic of CR based on the polarization and angle, i.e., the amplitude of MSP is a single-peak curve when the azimuth or pitch angle of an incident wave is monotonic within 35° without clutter. Based on the above property, the multi-characteristic between the ship and CR would not be affected by the angle of the incident wave, which is more applicable.
  • The CR discrimination based on the TSP joint domains. In contrast to References [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,42], we focus on analyzing the characteristic changes of a target during the radar movement, rather than being limited to a fixed pulse. That is, this paper utilizes multiple pulses for discrimination. Based on the three-dimensional proportional guidance method with the angle constraint, this paper gives the spatial-time characteristic of the radar, i.e., the azimuth and pitch angles of a radar’s line of sight (LOS) change monotonically as the radar approaches the ship. Subsequently, we introduce a definition of the CR-MSE parameter, and propose a novel CR discrimination method based on the TSP joint domains, which can achieve 95% discrimination probability and a 4.1% false alarm probability.
This paper is organized as follows. Section 2 analyzes the scattering model of the CR: Section 3 proposes a CR decomposition; Section 4 proposes a CR discrimination method based on TSP joint domains; Section 5 performs a mathematical simulation; and Section 6 concludes this paper.

2. The Scattering Model of CR

For the scattering model, the numerical algorithms usually have high accuracy, i.e., the method of moment (MoM) [43], and finite difference time domain (FDTD) [44], which request huge computational complexity. Herein, the radar frequency belongs to a higher frequency (i.e., X, Ku, K, Ka), and the size of CR and resolution are greater than 10 times the wavelength. Therefore, the high-frequency approximate algorithms, i.e., PO, GO, and SBR, have the advantages of high accuracy and computational complexity at the same time, which are suitable for analyzing the scattering model of CR [40]. In this paper, we will use the high-frequency approximate algorithms to obtain the scattering model of CR, and then analyze the polarization characteristic of CR, as shown below.
There is a target V which is a distance away from the target observation point O. The target center is O , and the O O vector is r . The incident wave source is Q, and the Q O vector is s . We randomly chose a triangular facet A B C on the target surface, where its normal vector is n ^ . The vector from O to the center of A B C is r . The electric field and magnetic field of the incident wave are E i , H i , respectively, and the electric field and magnetic field of the reflected wave are E r , H r , respectively. The scattering relationship is displayed in Figure 1.
The relationship between the incident wave and reflected wave can be obtained through the Stratton–Chu equation, i.e.,
E r r = [ j ω μ n ^ × H G r , r + n ^ × E × G r , r + n ^ · E G r , r ] d s H r r = [ j ω ε n ^ × E G r , r + n ^ × H × G r , r + n ^ · H G r , r ] d s ,
where ω is the angular frequency, ε is the dielectric constant of the propagation medium, μ is the magnetic permeability, and s is the area of A B C . E is the total electric field, i.e., E = E i + E r . H is the total magnetic field, i.e., H = H i + H r . G r , r is the Green’s function as follows.
G r , r = e j k r r 4 π r r ,
where k = 2 π / λ .
Subsequently, we can simplify the reflected wave’s electric field E r by the Huygens principle, viz.,
E r r = × n ^ × E G r , r d s + 1 j ω ε × × G r , r n ^ × H d s .
In the same way, the reflected wave’s magnetic field H r is
H r r = × n ^ × H G r , r d s 1 j ω μ × × G r , r n ^ × E d s .
Since the integral terms in Equations (3) and (4) are difficult to perform, we utilize high-frequency approximation methods for simplification.

2.1. Single Reflection

For the SR of CR, its electromagnetic scattering is analyzed using the Physical Optics (PO) [45]. First, the distance D R between O and A B C in Figure 1 is
D R = | r r | = | r | 1 r · r | r | 2 + r 2 2 | r | 2 ( r · r ) 2 2 | r | 4 + | r | r ^ · r ,
where O is far away from O , i.e., | r | > > | r | , and r ^ = r / | r | .
In the same way, the distance D S between Q and A B C is
D S = | s + r | | s | + s ^ · r ,
where s ^ = s / | s | .
Then, Green’s function is simplified by Equation (5), i.e.,
G r , r = exp j k r r 4 π r r = exp j k | r | exp j k r ^ · r 4 π | r | r ^ · r = exp j k | r | exp j k r ^ · r 4 π 1 | r | 1 + r ^ · r | r | + exp j k | r | exp j k r ^ · r 4 π | r | .
Therefore, E r can be further simplified [46], i.e.,
E r r = j k e j k | r | 4 π | r | { μ ε r ^ × r ^ × n ^ × H r ^ × n ^ × E } e j k r ^ · r d s .
On a perfect electric conductor, the area can be approximated as a plane. Then, the property is exhibited in Equation (9) and therefore the radius of curvature on the target surface and the nearby surface is much larger than the wavelength.
n ^ × E = 0 n ^ × H = 2 n ^ × H i .
Subsequently, E r can be obtained using Equations (8) and (9), viz.,
E r r = j k e j k | r | 2 π | r | μ ε r ^ × r ^ × n ^ × H i e j k r ^ · r d s ,
where H i is the incident wave’s magnetic field, i.e., H i = H 0 e j k · D S + 2 π f 0 t / D S . Furthermore, H i can be simplified by the far-field approximation in Equation (11).
H i = H 0 | s | e j k | s | + 2 π f 0 t e j k s ^ · r .
Then, let the substitute Equation (11) into (10), where E r is
E r r = j k e j k ( | r | + | s | ) + j 2 π f 0 t 2 π | r | | s | μ ε r ^ × r ^ × n ^ × H 0 e j k r ^ s ^ · r d s .
The properties of the plane electromagnetic waves are
H 0 = ε μ s ^ × E 0 E 0 = μ ε s ^ × H 0 .
Therefore, E r can be further simplified in Equation (14).
E r r = j k e j k ( | r | + | s | ) + j 2 π f 0 t 2 π | r | | s | I n ^ · E 0 r ^ · s ^ r ^ s ^ n ^ · s ^ r ^ · E 0 r ^ E 0 ,
where the integral part I [47] is
I = e j k r ^ s ^ · r d s = 1 j k r ^ s ^ · n ^ 2 | r ^ s ^ | 2 · m = 1 M I { r ^ s ^ × n ^ · a m a m 1 · exp j k 2 r ^ s ^ · a m + a m 1 · sin k 2 r ^ s ^ · a m a m 1 k 2 r ^ s ^ · a m a m 1 } ,
where M I is determined according to the number of sides of A B C , i.e., M I = 3 . a m is the vertice of A B C , m = 1 , 2 , , M I .
In brief, Equation (14) is the electric field of a bistatic radar. In the case of a monostatic radar, r = s , and the property holds, i.e., r ^ · E 0 = s ^ · E 0 = 0 . Therefore, the SR of CR can be obtained by simplifying Equation (14), as displayed in Equation (16).
E single r = j k e j 2 π f 0 ( t 2 | r | / c ) 2 π | r | 2 I single n ^ · r ^ E 0 ,
where c is the light velocity and the integral part I single is
I single = j r ^ × n ^ 2 k 1 r ^ · n ^ 2 · m = 1 M { a m a m 1 exp j k r ^ · a m + a m 1 · sin k r ^ · a m a m 1 k r ^ · a m a m 1 } .

2.2. Double Reflection

In this section, we combine PO and SBR [48] together, which can separately analyze DR and TR for the CR of any shape. This method based on PO-SBR contains two steps, i.e., ray tracing and electromagnetic calculation. In the process of ray tracing, the CR is decomposed into a series of triangular facets, where ray tracing is performed on the triangular facets. Let A 2 be the triangular facet of DR, and A 1 be the triangular facet of SR. When the incident wave with the unit vector s 0 ^ illuminates the triangular facet A 1 with the normal vector n 1 ^ , the unit vector of the reflected wave is r 0 ^ , as presented in Figure 2. In addition, a s , a r are the incident and reflected angles of A 1 , respectively.
First, Equation (18) can be given according to specular reflection.
r 0 ^ × n 1 ^ = s 0 ^ × n 1 ^ .
By multiplying Equation (18) by ( r 0 ^ + s 0 ^ ) and n 1 ^ , respectively, the properties of r 0 ^ and s 0 ^ are
r 0 ^ = s 0 ^ 2 s 0 ^ · n 1 ^ n 1 ^ s 0 ^ = r 0 ^ 2 r 0 ^ · n 1 ^ n 1 ^ .
Then, by substituting Equations (19) into (14), E 1 ( r 0 ) is
E 1 r 0 = j k e j k ( | r 0 | + | s 0 | ) + j 2 π f 0 t 2 π | r 0 | | s 0 | I · [ 2 r 0 ^ · n 1 ^ 2 E 0 · n 1 ^ r 0 ^ + 2 r 0 ^ · n 1 ^ E 0 · n 1 ^ n 1 ^ + r 0 ^ · n 1 ^ r 0 ^ · E 0 r 0 ^ r 0 ^ · n 1 ^ E 0 ] .
Since the incident wave is a plane wave, i.e., s ^ · E 0 = 0 , E 1 ( r 0 ) can be further simplified in Equation (21).
E 1 r 0 = j k e j k ( | r 0 | + | s 0 | ) + j 2 π f 0 t 2 π | r 0 | | s 0 | I r 0 ^ · n 1 ^ 2 E 0 · n 1 ^ n 1 ^ E 0 .
If the phase part j e j k ( | r 0 | + | s 0 | ) + j 2 π f 0 t and the amplitude part k I r 0 ^ · n 1 ^ / 2 π | r 0 | | s 0 | are ignored in Equation (21), the polarization information after SR is
E 1 = 2 E 0 · n 1 ^ n 1 ^ E 0 .
In Figure 2, A 1 is the triangular facet of SR, and Q 1 is the incident wave source. Based on the equivalent image source method, we let Q 2 be the mirror image source. Then, Q 2 can be regarded as the incident wave source of DR, and s 1 can also be regarded as the incident source to the target center. Using the geometric relationship, we find that the vector from Q 2 to the center of A 1 and the vector from the center of A 1 to the center of A 2 are the same. Therefore, s + r A 1 = s 0 , and s 1 + r A 1 = r 0 . According to Equation (19), s 1 after the SR is
s 1 = s 2 ( s + r A 1 ) · n 1 ^ n 1 ^ .
Subsequently, multiply the left and right sides of Equation (23) by ( I 2 n 1 ^ n 1 ^ T ) , and | s 1 | can be obtained by utilizing | n 1 ^ | = 1 , i.e.,
| s 1 | = | s 2 ( E I 2 n 1 ^ n 1 ^ T ) n 1 ^ n 1 ^ T r A 1 | | s | 2 s ^ · ( E I 2 n 1 ^ n 1 ^ T ) n 1 ^ n 1 ^ T r A 1 ,
where E I is unit matrix, and the unit vector of s 1 is
s 1 ^ = s 1 | s 1 | ( E I 2 n 1 ^ n 1 ^ T ) 1 s | s | = ( E I 2 n 1 ^ n 1 ^ T ) s ^ .
Finally, the DR of CR can be obtained by substituting Equations (22), (24) and (25) into (14), i.e.,
E double r = j k e j k ( | r | + | s | ) + j 2 π f 0 t 2 π | r | | s | · I double e 2 j k s ^ · ( E I 2 n 1 ^ n 1 ^ T ) n 1 ^ n 1 ^ T r A 1 · { n 2 ^ · E 1 r ^ · s ^ 1 r ^ s ^ 1 n 2 ^ · s ^ 1 r ^ · E 1 r ^ E 1 } ,
where n 2 ^ is the unit normal vector of the DR’s triangle facet, and I double = e j k r ^ s ^ 1 · r d s .

2.3. Triple Reflection

Similarly to DR, let A 1 be the triangular facet of SR with the normal vector n 1 , where the incident wave and reflected wave are s and s 1 , respectively. The triangular facet of DR is A 2 with the normal vector n 2 , where the incident wave and reflected wave are s 1 and s 2 , respectively, as exhibited in Figure 3.
First, the polarization information of DR is displayed in Equation (27) by referring to Equation (22).
E 2 = 2 E 1 · n 2 ^ n 2 ^ E 1 .
Similarly referring to Equation (24), s 2 after the DR is
| s 2 | | s | 2 s ^ · m = 1 2 i = 1 m ( I 2 n i ^ n i ^ T ) n m ^ n m ^ T r A m .
Subsequently, the unit vector of s 2 is
s 2 ^ = s 2 | s 2 | i = 1 2 ( I 2 n i ^ n i ^ T ) 1 s | s | = ( I 2 n 2 ^ n 2 ^ T ) ( I 2 n 1 ^ n 1 ^ T ) s ^ .
Finally, the TR of CR can be obtained by substituting Equations (27)–(29) into (14), i.e.,
E triple r = j k e j k ( | r | + | s | ) + j 2 π f 0 t 2 π | r | | s | I triple × e 2 j k s ^ · m = 1 2 i = 1 m ( I 2 n i ^ n i ^ T ) n m ^ n m ^ T r A m × n 3 ^ · E 2 r ^ · s ^ 2 r ^ s ^ 2 n 3 · s ^ 2 r ^ · E 2 r ^ E 2 ,
where n 3 ^ is the unit normal vector of TR’s triangular facet, and I triple = s e j k r ^ s ^ 2 · r d s .

3. The CR Decomposition Based on the Main-Polarization

The SR, DR, and TR model of CR are shown in Equations (16), (26) and (30), respectively. When the incident wave’s polarization changes, both the phase and amplitude do not change. Thus, by ignoring the phase and amplitude, the polarizations of SR, DR, and TR are only related to
E single r E 0 E double r n 2 ^ · E 1 r ^ · s ^ 1 r ^ s ^ 1 n 2 ^ · s ^ 1 r ^ · E 1 r ^ E 1 E triple r n 3 ^ · E 2 r ^ · s ^ 2 r ^ s ^ 2 n 3 ^ · s ^ 2 r ^ · E 2 r ^ E 2 .
First, let the horizontal polarization be E ϕ and the vertical polarization be E θ . Taking SR in Equation (31) as an example, the reflected wave is E ϕ when the incident wave is E ϕ , and the reflected wave is E θ when the incident wave is E θ . Therefore, the PSM of SR is [ 1 0 ; 0 1 ] . Similarly, the PSMs of SR, DR, and TR are
E single = 1 0 0 1 E double = U p i p j 1 0 χ p i p j 1 U p i p j 1 E triple = 1 0 0 1 ,
where χ p i p i is determined by the ratio of DR from plane p i to plane p j and from plane p j to plane p i , U p i p j is
U p i p j = cos ψ p i p j sin ψ p i p j sin ψ p i p j cos ψ p i p j ,
where ψ p i p j is the rotation angle for DR from plane p i to plane p j .
Then, we denote the odd scattering as the sum of SR and TR, and symbolize the even scattering as the DR. References [27,28,29] obtained the odd scattering and even scattering by the Krogager decomposition [49], which decomposes the scattering matrix into spherical scattering (odd scattering), dihedral scattering (even scattering), and helix scattering, i.e.,
S = e j ϕ 0 { e j ϕ S k S 1 0 0 1 + k D cos 2 θ D sin 2 θ D sin 2 θ D cos 2 θ D + k H e j 2 θ D 1 j j 1 } ,
where k S , k D and k H are the coefficients of spherical scattering, dihedral scattering, and helix scattering, respectively. ϕ 0 is the absolute phase of scattering matrix, ϕ S is the relative phase of spherical scattering to dihedral scattering, and θ D is the rotation of the target relative to the horizontal polarization.
Comparing Equations (32) and (34), we find that the Krogager decomposition has the following assumptions when performing the odd scattering and even scattering of CR.
  • It assumes that the helix scattering is a noise after removing the odd scattering and even scattering, whereas the helix scattering is poorly representative of clutter;
  • It assumes that the relative phase of dihedral scattering and helix scattering is 0;
  • It assumes that the dihedral scattering only has a single dihedral scattering. Since there are three kinds of dihedral scatterings in trihedral CR, the dihedral scattering in Equation (34) cannot characterize the even scattering of CR.
In brief, the coefficient errors of odd scattering and even scattering after Krogager decomposition would increase as SCR decreases. To resolve the aforementioned problem, we propose a novel CR decomposition by the main-polarization as follows.
First, we denote the CR and clutter scattering matrix as S C and S N , respectively. The scattering matrix of echo is
S = S C + S N = S H H S H V S V H S V V ,
where S H V means the horizontal transmitting and vertical receiving polarization, and the others are defined similarly.
For the scattering matrix S of i-th pulse, we transform S to Kennaugh matrix [50] as shown in Equation (36).
K i = A 0 + B 0 C k H k F k C k A 0 + B k E k G k H k E k A 0 + B k D k F k G k D k A 0 + B 0 ,
where A 0 , B 0 , B k , C k , D k , E k , F k , G k are the Huynen parameters [51].
Then, the average Kennaugh matrix during N pulses is
K = 1 N i = 1 N K i .
Subsequently, the average Kennaugh matrix K can be decomposed by Yang decomposition [52], i.e.,
K = K 0 + K n ,
where K 0 is the Kennaugh matrix of CR, and K n is the Kennaugh matrix of clutter after pulse accumulation.
Then, K 0 is converted into a coherent matrix T 0 , viz.,
T 0 = 2 A 0 + B 0 C k j D k H k + j G k C k + j D k B 0 + B k E k + j F k H k j G k E k j F k B 0 B k .
Then, the Pauli vector k p = [ k 1 k 2 k 3 ] T is given using the coherence matrix T 0 , as shown in Equation (40).
T 0 = | k 1 | 2 k 1 k 2 * k 1 k 3 * k 2 k 1 * | k 2 | 2 k 2 k 3 * k 3 k 1 * k 3 k 2 * | k 3 | 2 ,
where k i * is the conjugate value of k i .
Therefore, the scattering matrix S k of CR after denoising can be obtained, i.e.,
S k = k 1 S p 1 + k 2 S p 2 + k 3 S p 3 ,
where S p 1 , S p 2 , and S p 3 are the fundamental scatterings of Pauli decomposition [50].
In the above process, S k is obtained using the Pauli vector k p = [ k 1 k 2 k 3 ] T , in which the coefficients are from the coherence matrix T 0 . By comparing the original scattering matrix S C , S k only has an absolute phase difference, which does not affect the coefficients of the fundamental scatterings.
The trihedral CR mainly includes four basic scatterings, i.e., SR, DR, TR, and edge diffraction. When the radar incident wave irradiates the central axis of CR, TR is dominant. However, SR, DR, and edge diffraction are dominant when the radar incident wave irradiates the edge or vertical plane of CR. Since the edge diffraction is still 20 dB lower than SR and DR [53,54], the decomposition of CR is exhibited in Equation (42) after ignoring the influence of edge diffraction.
S k = S S R + S T R + S D R = S odd + S even ,
where S S R , S D R , S T R are the SR, DR, and TR, respectively. S odd , S even are odd scattering and even scattering, respectively.
For the CR odd scattering, it consists of SR and TR. Based on Equation (32), the odd scattering is
S odd = S S R + S T R = k S R e j θ S R 1 0 0 1 + k T R e j θ T R 1 0 0 1 = k odd e j θ odd 1 0 0 1 ,
where k S R , k T R are the coefficients of SR and TR, respectively. θ S R , θ T R are the phases of SR and TR, respectively. k odd and θ odd are the coefficient and phase of odd scattering, respectively.
For the CR even scattering, it consists of three kinds of dihedral scatterings. Based on Equation (32), the even scattering is
S even = S D R 12 + S D R 23 + S D R 31 = k D R 12 e j θ D R 12 U 12 1 0 χ 12 1 U 12 1 + k D R 23 e j θ D R 23 U 23 1 0 χ 23 1 U 23 1 + k D R 31 e j θ D R 31 U 31 1 0 χ 31 1 U 31 1 = k even e j θ even 1 k cross 1 e j Δ θ 1 k cross 2 e j Δ θ 2 1 ,
where S D R p i p j is the DR of planes p i and p j , { p i , p j = 1 , 2 , 3 } . k D R p i p j , θ D R p i p j are the coefficient and phase of S D R p i p j , respectively. k even , θ even are the even scattering’s coefficient and phase, respectively. k cross i , Δ θ i are the cross-polarization’s coefficient and phase, respectively.
After analyzing the basic scatterings of CR, the CR DR is formed by three kinds of dihedral scattering as shown in Equation (44). By comparing Equations (34) and (44), we find that the dihedral scattering in Krogager decomposition cannot accurately represent the CR DR. To resolve the above problem, we optimize the decomposition of CR by utilizing Equations (43) and (44), and then propose the CR decomposition method based on main polarization as follows.
By combining Equations (43) and (44), the scattering matrix S k of CR after denoising is
S k = S odd + S even = k odd e j θ odd 1 0 0 1 + k even e j θ even 1 k cross 1 e j Δ θ 1 k cross 2 e j Δ θ 2 1 .
Subsequently, we utilize the main polarization of S k to analyze the odd scattering and even scattering, i.e.,
k odd e j θ odd + k even e j θ even = S k 1 , 1 k odd e j θ odd k even e j θ even = S k 2 , 2 ,
where S k x , y is the element of the x-th row and the y-th column of the matrix S k .
Based on Equations (36)–(41), (45), and (46), the coefficients k odd and k even of CR odd scattering and even scattering can be quickly obtained, thus completing the CR decomposition. The algorithm flow chart of CR decomposition is shown in Figure 4. First, we transform the Sinclair matrix into a Kennaugh matrix, and perform denoising by utilizing pulse accumulation and Yang decomposition. Then, we transform K 0 to the Sinclair matrix by the Pauli vector. Finally, the coefficients of odd scattering and even scattering can be obtained by the main-polarization decomposition. With the CR decomposition, we cannot only cure the shortcomings of Krogager decomposition, but also accurately obtain the odd scattering and even scattering of CR in clutter.

4. The CR Discrimination Method Based on TSP Joint Domains

A new CR decomposition method is proposed in Section 3, and we can obtain the polarization characteristic of CR according to the above new CR decomposition method. On this basis, we integrate the spatial–time characteristic of the radar and target into the polarization characteristic, thereby proposing a discrimination method of CR, the details of which are as follows.

4.1. The Polarization Characteristic

With the coefficients of the CR odd scattering and even scattering, this section will analyze the polarization characteristic of CR, thus facilitating to the distinction between CR and ship. First, we introduce the definition of MSP as follows.
Definition 1.
The CR mainly includes SR, DR, and TR. The SR and TR belong to the odd scattering, and DR belongs to the even scattering. At a certain azimuth and pitch angle, the even scattering is the main scattering polarization (MSP) if the even scattering is stronger than odd scattering. On the contrary, the odd scattering is MSP.
By analyzing the basic scatterings of odd scattering and even scattering in trihedral CR, we find that the SR, DR, and TR are all single-peak functions, the analysis and simulation of which are shown in Appendix A. Combining the changing rule of odd scattering and even scattering, the polarization characteristic is given in Appendix A, i.e., the MSP coefficient of CR is a single-peak curve, when the azimuth or pitch angle of radar incident wave changes monotonically within 35°.

4.2. The Spatial–Time Characteristic

Section 4.1 gives the polarization characteristic of CR, which is a characteristic which has come from the mutual motion of the radar and target. For a moment of stillness, the aforementioned polarization characteristic does not exist. Therefore, we need to consider the mutual movement of the radar and target, and then ensure the existence condition of the polarization characteristic in Section 4.1. In the spatial–time characteristic of the radar and target, the radar usually detects a ship on the sea surface from the air, and finally attacks the ship in the form of diving to increase the effect and speed. Actually, a large pitch angle attack is required when hitting the target, i.e., striking a ship from a nearly vertical direction, thus improving the strike capability. This section will take the large pitch angle attack as an example to analyze the spatial-time characteristic of the radar and target, which is based on the three-dimensional proportional guidance method.
The spatial relationship of the radar and target is shown in Figure 5. The radar’s position is O, the target is T 0 , and the LOS between the radar and target is O T 0 . O X g Y g Z g , O X t Y t Z t is the radar and target’s inertial coordinate system, respectively. φ L , θ L are the LOS azimuth and pitch angle, respectively. O X l Y l Z l is the LOS coordinate system, where O X l is consistent with the LOS direction, O Z l is downward in the vertical plane and also orthogonal to O X l , O Y l is determined based on the right-hand rule. The velocity of radar is V O , and φ V , θ V are the velocity azimuth and pitch angle, respectively. O X v Y v Z v is the radar velocity’s coordinate system, where O X v is consistent with the velocity direction, O Z v is downward in the vertical plane as well as orthogonal to O X v , and O Y v is determined by the right-hand rule. The velocity of the target is V t , and its azimuth and pitch angle are φ T and θ T , respectively.
There are three coordinate systems in Figure 5, i.e., the inertial coordinate system O X g Y g Z g , the LOS coordinate system O X l Y l Z l , and the velocity coordinate system O X v Y v Z v . The rotation matrices between the three coordinate systems are exhibited in Equation (47), which can be used to transform O X g Y g Z g into O X l Y l Z l or O X v Y v Z v .
C g l = cos θ L cos φ L cos θ L sin φ L sin θ L sin φ L cos φ L 0 sin θ L cos φ L sin θ L sin φ L cos θ L C g v = cos θ V cos φ V cos θ V sin φ V sin θ V sin φ V cos φ V 0 sin θ V cos φ V sin θ V sin φ V cos θ V .
For a brief description, we take the pitch plane as an example in the following. First, an offset term of angular rate is
w θ = η θ θ d θ L N θ t = η θ v x L O S θ d θ L N θ r ,
where η θ and N θ both are guidance constants, θ d is the constrained pitch angle when striking target, and t is the remaining time.
Subsequently, the LOS velocity v x L O S can be obtained as shown in Equation (49), i.e., by transferring the radar velocity from the velocity coordinate system to an LOS coordinate system.
v x L O S v y L O S v z L O S = C g l C g v T v x o v y o v z o ,
where v x o , v y o , v z o are the components of radar velocity in a velocity coordinate system, respectively.
Therefore, the injunction of a pitch rate in the proportional guidance is
Ω θ = N θ w θ L O S w θ ,
where w θ L O S is the actual change rate of the LOS pitch angle.
In the same way, the injunction of an azimuth rate in the proportional guidance is
Ω φ = N φ w φ L O S η φ v x L O S φ d φ L N φ r ,
where η φ and N φ both are navigation constants, φ d is the constrained azimuth angle when striking the target, and w φ L O S is the actual change rate of the LOS azimuth angle.
In brief, the accelerated velocity of a radar can be controlled by iteratively updating Ω θ and Ω φ with Equations (50) and (51). Therefore, we can successfully perform the three-dimensional proportional guidance with angle constraints, which can control the radar-to-strike target at a constrained angle. Now, we take an example to briefly explain the spatial–time characteristic in a terminal guidance process.
Example 1.
The LOS distance between radar and ship is 20 km, the initial LOS azimuth angle φ L is 30 ° , and the initial LOS pitch angle θ L is 45 ° . The radar velocity V O is 1000 m / s , the radar velocity’s azimuth angle φ V is 10 ° , and the radar velocity’s pitch angle θ V is 40 ° . The ship velocity V t is 15 m / s , the ship velocity’s azimuth angle φ T is 20 ° , and the ship velocity’s pitch angle θ T is 0 ° .
In the three-dimensional proportional guidance, the pitch constraint angle θ d is 90 ° , and the azimuth constraint angle φ d is not limited. The trajectory of the radar is displayed in Figure 6, and the changes in the LOS azimuth angle and pitch angle are shown in Figure 7. We find that both the azimuth angle and pitch angle are monotonic changing functions. To sufficiently demonstrate the above conclusion, we randomly adjusted the initial LOS azimuth angle φ L , the initial LOS pitch angle θ L , the radar velocity’s azimuth angle φ V , the radar velocity’s pitch angle θ V , and the ship velocity’s azimuth angle φ T . After 10,000 Monte Carlo simulations, we also found that the LOS azimuth angle and pitch angle both change monotonically when the radar approaches the target, which is the spatial–time characteristic of the radar in the terminal guidance process.

4.3. The Discrimination Method

With the analyses of the polarization characteristic and the spatial–time characteristic, we find that the MSP of CR is a single-peak curve when radar approaches the ship. In contrast, the ship’s odd scattering and even scattering will change irregularly as the LOS changes, and because of this, the ship’s strong scattering points are unstable. Accordingly, the MSP of a ship is not a single-peak curve, which is different to that of CR. In order to quantitively describe the MSP amplitude curve, this paper introduces a CR-MSE parameter as follows.
First, we carry out the MSP amplitude within a period of time, and then utilize the quadratic function for curve fitting. The MSE between the MSP and quadratic function is
S MSE = 1 N i = 1 N P i f i 2 ,
where N is the number of target echoes after pulse accumulation, P i is the MSP amplitude, and f i is the quadratic function of MSP.
Due to fact that the MSP of CR is a single-peak curve and that the MSP of a ship changes irregularly, the S MSE of CR is much smaller than that of a ship. Therefore, S MSE can be utilized as a parameter to discriminate CRs and ships, i.e.,
Target = Corner Reflector , S MSE < T value Ship , S MSE T value ,
where T value is the judgment threshold. Since the judgment threshold is related to the polarization characteristic of CR, it is difficult to give a numerical expression from the perspective of probability distribution. Therefore, this paper analyzes the judgment thresholds under different discrimination probabilities and different false alarm probabilities by means of a large number of simulations. The details are given in Section 5.3.
Taking the influence of clutter into account, we propose a CR discrimination method, which is composed of the CR decomposition and Equations (52) and (53). The algorithm flow chart is shown in Figure 8, and the detailed steps are as follows.
Step 1: For the scattering matrices of M pulses in a period of time, we perform pulse accumulation every N pulses;
Step 2: The M / N MSPs are given by the CR decomposition using Equations (36)–(41), (45), and (46);
Step 3: The MSP is fitted by the quadratic function, so as to calculate S MSE using Equation (52);
Step 4: The ship and CR can be discriminated by their S MSE using Equation (53), where it is judged as a CR if S MSE < T value , otherwise it is a ship.

5. Simulation

5.1. The CR Polarization Characteristic

In Section 4, we analyzed that the CR MSP is a single-peak curve when the azimuth or pitch angle changes monotonically within 35 ° without clutter. Therefore, this section will demonstrate the above conclusion by analyzing the MSP amplitude under any angle change.
Example 2.
The length, width, and height of CR are all 0.25 m. The frequency of incident wave is 12 GHz. Therefore, the wavelength λ is 0.025 m, and the dimension of CR satisfies 10 λ . Subsequently, the surface of CR is divided by a triangular surface, and the length of triangular surface is less than λ / 10 . The range of azimuth and pitch angle are both 0 ° , 90 ° , and the angular interval is 1 ° . With the change in azimuth and pitch angle, the coefficients of odd scattering and even scattering are calculated.
First, we randomly choose the initial pitch angle θ and azimuth angle φ , i.e., θ = 50 ° and φ = 20 ° . The monotonic change curves of the pitch angle and azimuth angle are randomly given and are satisfied within 35 ° , as exhibited in Figure 9. Subsequently, we calculate the CR amplitude changes of different scatterings, i.e., the odd scattering, even scattering, and MSP, as shown in Figure 10. However, the MSP is obtained by Equation (54).
P MSP ( i ) = odd i , i f odd i > even i P MSP ( i ) = even i , i f odd i even i ,
where P MSP ( i ) is MSP. odd i , even i are the coefficients of odd scattering and even scattering, respectively.
Then, we analyze the number of MSP changes by odd scattering or even scattering, viz.,
P change = i = 2 M η ( i ) η ( i ) = 1 , i f k i 1 k i < 0 η ( i ) = 0 , e l s e ,
where M is the number of pulses required for discrimination and k i represents whether MSP is composed of odd scattering or even scattering, as shown in Equation (56).
k i = 1 , i f odd i > even i k i = 1 , i f odd i even i .
Finally, the monotonic changing curves of the pitch angle and azimuth angle are randomly generated by utilizing 10,000 Monte Carlo simulations, which can satisfy the changing range within 35 ° . By Equations (54)–(56), the numbers of MSP changing P change can be obtained as displayed in Table 1. The results show that the changes of MSP are mainly 0 and 1. If P change = 0 , MSP is only odd scattering or even scattering. If P change = 1 , MSP is composed of “odd scattering + even scattering” or “even scattering + odd scattering”. Because of this, the odd scattering and even scattering are both monotonic functions, and the changing trend of MSP is mainly a monotonic function (i.e., a special single-peak function) or single-peak function, thus demonstrating the polarization characteristic of CR.

5.2. The CR Discrimination Method

Focusing on the CR discrimination method based on TPS joint domains, this section will analyze the discrimination probability under different observing time and different SCRs.
Example 3.
The simulation model of CR is the same as that of Example 2. Now, the same triangular facet division and incident wave frequency are used to simulate the ship model. The ship’s length is 116 m, width is 28 m, and height is 39 m, as shown in Figure 11. Subsequently, the trajectories of radar and target are simulated. The LOS distance between the radar and ship is 20 k m , the radar velocity | V O | is 1000 m / s , the ship velocity | V t | is 15 m / s , and the direction parameters are presented in Table 2. Among them, the radar velocity’s azimuth angle φ V is randomly selected within 30 ° more or less than the initial LOS azimuth angle φ L . The radar strikes the ship with a three-dimensional proportional guidance method based on angle constraints, where the pitch constraint angle θ d is 90 ° , and the azimuth constraint angle φ d is not limited.
First, let the observing time be 10 s, and randomly select the direction parameters. Then, we perform 10,000 Monte Carlo simulations to calculate the CR-MSEs of CR and ship according to Equation (52), as displayed in Figure 12. The result shows that the average CR CR-MSE is 0.0003, whereas that of the ship is 0.0041. Therefore, the CR CR-MSE is much smaller than that of ship, which demonstrates that the CR and ship can be discriminated based on CR-MSE. By changing the judgment threshold with a large number of simulations, the discrimination probability and false alarm probability are exhibited in Figure 13. As the discrimination probability increases, the false alarm probability increases. Therefore, we can determine the judgment threshold based on Figure 13, when the discrimination probability and false alarm probability both are fixed. Regarding the discrimination of corner reflectors, all methods (i.e., HRRP, micro-Doppler, and polarization methods in the Introduction) need to obtain a threshold level using simulation or actual data. The method in this paper is the same as the above methods, which fails to give a closed-form formula for the threshold T v a l u e . However, compared with other methods, the threshold in this paper is obtained by combining the spatial–time characteristic and polarization characteristic. Therefore, the adaptability of the threshold T v a l u e is more extensive. For example, for different scenarios with different parameters, i.e., the initial LOS, the velocity of the radar and the velocity of the ship, the discrimination probability, and false alarm probability of the present method are better than those of other methods, as shown in Example 5.
Subsequently, the observing time is chosen to be between 1–15 s, respectively. Based on the 10,000 Monte Carlo simulations, the judgment threshold is counted under the premise that the discrimination probability is 95 % . Therefore, the false alarm probability can be calculated, as shown in Figure 14. As the observing time increases, the false alarm probability first decreases and then increases. The result can indicate that the fluctuation of the ship’s MSP is not obvious when the observing time is short. At this time, the judgment threshold is low, which will increase the false alarm probability of CR. Nevertheless, the changing trend of the CR MSP will fluctuate when the observing time is long, which also increases the false alarm probability of the CR. Therefore, the difference between the CR and ship is the most obvious when the observing time is 7 s under the premise of ensuring the discrimination probability of 95 % . At this time, the false alarm probability is 4.1 % .
After determining the observation time and threshold, we will analyze the discrimination probability and false alarm probability under different SCRs. Let the SCR change in the scope of 0 dB , 20 dB , whilst the discrimination probability and false alarm probability are shown in Figure 15 and Figure 16. The discrimination probability and false alarm probability under different SCRs are mainly closely related to the denoising effect of CR decomposition. In Figure 17, we will find that the CR decomposition has a better denoising effect when the SCR is higher than 7 dB. Therefore, the discrimination probability is higher than 94 % and the false alarm probability is lower than 7 % when the SCR is higher than 7 dB, as exhibited in Figure 15 and Figure 16. As the SCR increases, the discrimination probability will tend towards 95 % and the false alarm probability will tend to 4.1 % , which is consistent with the theoretical analysis of the CR decomposition, thus demonstrating that the present method is effective in clutter.

5.3. The Comparisons of the Proposed Work

In this paper, we proposed a new polarization decomposition and a CR discrimination method. Therefore, we will give the comparisons of the proposed work in this section.
First, as for the new polarization decomposition, we will compare the coefficient errors of odd scattering and even scattering, which are based on the Krogager decomposition, main polarization decomposition, and the present method in this paper under different SCRs.
Example 4.
The simulation parameters are the same as that of Example 2. The range of the azimuth and pitch angle are both 0 ° , 90 ° , and the angular interval is 1 ° . The scattering matrices of SR, DR, and TR were obtained by the MoM in Feko software.
After calculating the true values of odd scattering and even scattering, the clutter is added into four polarization channels, respectively. Subsequently, the coefficients of odd scattering and even scattering are performed by Krogager decomposition, main-polarization decomposition, and the present method, respectively. The coefficient errors are compared by Equation (57) as follows.
error = i = 1 90 j = 1 90 C odd i , j D odd i , j 2 + C even i , j D even i , j 2 3600 ,
where C odd i , j , C even i , j are the true values of odd scattering and even scattering under different azimuth and pitch angles, respectively. D odd i , j , D even i , j are the odd scattering and even scattering coefficients by the decomposition method, respectively.
In order to analyze the influence of clutter, SCR will change in the scope of 0 dB , 20 dB . Then, the error coefficients of Krogager decomposition, main-polarization decomposition, and the present method are displayed in Figure 17. By comparing the Krogager decomposition and main-polarization decomposition, the errors of two decomposition methods are not much different when SCR is small. However, the main-polarization decomposition has a smaller error than Krogager decomposition when the SCR is higher than 5 dB. This result shows that the dihedral scattering of Krogager decomposition cannot accurately represent the CR even scattering due it being constituted by three different DRs as shown in Equation (44). By observing the present method in Figure 17, its error is obviously better than that of Krogager decomposition and main-polarization decomposition when the SCR is higher than 7 dB, due to the CR decomposition utilizing the denoise function of the Yang decomposition. When SCR > 7 dB, the Yang decomposition can better remove noise, but its denoising effect is not obvious when the SCR is too low. Therefore, the compared results can demonstrate that the present method has a better denoising effect and is more adaptive to decompose the CR odd scattering and even scattering in clutter.
Secondly, we will compare the CR discrimination method in this paper with the HRRP method [14] and the polarization method [31].
Example 5.
The ship model, ship velocity, radar velocity, and relative range are the same as that of Example 3. The array CRs are composed by four CRs, for which the relative locations are [ 0 m , 29 m , 63 m , 117 m ] and the relative RCS is [ 1 , 0.8 , 2 , 1 ] . For the parameters of the radar, its frequency is 12 GHz, its bandwidth is 30 MHz, its pulse width is 2 μs, sampling frequency is 100 MHz, and SCR is 18 dB.
In [14], the extreme learning machine algorithm is proposed to discriminate between ship and CR based on the HRRP differences. For example, the HRRPs of the ship and array CRs are displayed in Figure 18 and Figure 19. In [31], six polarization parameters are utilized to discriminate between the ship and CR, i.e., the determinant value of PSM, trace of PSM, trace of power matrix, determinant value of power matrix, polarization shape factor, and combined invariant coefficient. In 10,000 Monte Carlo trials, we randomly selected the space and velocity parameters of the radar and target, i.e., the initial LOS, the velocity of the radar and the velocity of the ship. By performing 10,000 Monte Carlo simulations, the discrimination probability and false alarm probability of the three methods are exhibited in Table 3. For the HRRP method, the results are very susceptible to the arrangement between the CR and type of ship, thus having the highest false alarm probability. For the polarization method, it only utilizes the polarization information at an independent moment, which has an error when the incident wave is far from the LOS symmetry axis. In contrast, the method presented in this paper has comprehensively utilized the time, spatial, and polarization information. With the increases in the use of discriminative information and the expansion of the incident angle, the present method has higher discrimination probability and lower false alarm probability.

6. Conclusions

In this paper, we proposed a novel discrimination method based on high-frequency electromagnetic theory to discriminate the CR and ship. First, we proposed CR decomposition based on the main polarization, which could accurately obtain the coefficients of CR odd scattering and even scattering in clutter. Then, we gave the polarization characteristic of CR, i.e., the MSP amplitude of CR was a single-peak curve under ideal situation without clutter, when the azimuth or pitch angle of an incident wave changed monotonically within 35 ° . Subsequently, we analyzed the spatial–time characteristic based on the three-dimensional proportional guidance. Finally, we proposed a novel CR discrimination method based on the TSP joint domains by constructing the CR-MSE parameter, which could achieve 95 % discrimination probability and 4.1 % false alarm probability.

Author Contributions

Conceptualization, H.Y.; Investigation, H.H.; Methodology, Y.H. and J.Y. (Jian Yang); Writing—review & editing, J.Y. (Junjun Yin). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The Polarization Characteristic

For the monostatic radar detecting CR, the incident wave is r, the azimuth angle is φ , and the pitch angle is θ , as shown in Figure A1. Let the planes O A B , O B C , and O C A be planes 1, 2, and 3, respectively. Then, S S R 1 , S S R 2 , and S S R 3 are the SR of planes 1, 2, and 3, respectively. Among planes 1 and 2, planes 2 and 3, and planes 3 and 1, the DR in between are taken as S D R 12 , S D R 23 , and S D R 31 , respectively. The TR is S T R .
Figure A1. Schematic diagram of the incident wave.
Figure A1. Schematic diagram of the incident wave.
Remotesensing 14 03893 g0a1
First, the odd scattering matrix of CR is
S odd = S S R 1 + S S R 2 + S S R 3 + S T R .
As for the S S R 1 of plane 1, the coefficient of S S R 1 is the largest when the incident wave is perpendicular to plane 1. Moreover, S S R 1 is a monotonically decreasing function as the incident wave moves away from the normal line of the plane 1. In the same way, we find that both S S R 2 and S S R 3 are monotonically decreasing functions. Additionally, the maximum value of S S R 2 is located on the normal line of plane 2, and the maximum value of S S R 3 is located on the normal line of plane 3. As for S T R , the coefficient of S T R is the largest when the incident wave is located on the LOS symmetry axis of CR (i.e., φ = 45 ° and θ = 54 ° ). When the incident wave is far from the LOS symmetry axis, S T R is also a monotonically decreasing function.
Subsequently, the even scattering matrix of CR is
S even = S D R 12 + S D R 23 + S D R 31 .
For S D R 12 between planes 1 and 2, the coefficient of S D R 12 is the largest when the incident wave is on the symmetry plane of dihedral CR (i.e., A B C O A ) and is also parallel to the bottom surface O C A . For the sake of brevity, we define the symmetrical parallel axis of dihedral CR, i.e., it is located on the symmetry plane of dihedral CR and is also parallel to the bottom surface at the same time. As the incident wave moves away from the symmetrical parallel axis of planes 1 and 2, S D R 12 is a monotonically decreasing function. In the same way, S D R 23 and S D R 31 are both monotonically decreasing functions. Meanwhile, the maximum value of S D R 23 is located on the symmetrical parallel axis of planes 2 and 3, and the maximum value of S D R 31 is located on the symmetrical parallel axis of planes 3 and 1.
In summary, by analyzing the basic scatterings of odd scattering and even scattering in trihedral CR, we find that the SR (i.e., S S R 1 , S S R 2 , and S S R 3 ), DR (i.e., S D R 12 , S D R 23 , and S D R 31 ), and TR (i.e., S T R ) are all single-peak functions. To further analyze the changing rule of odd scattering and even scattering, we give a simulation in the following.
Example A1.
The length, width, and height of CR are all 0.25 m. The frequency of the incident wave is 12 G H z . Therefore, the wavelength λ is 0.025 m, and the dimension of CR satisfies 10 λ . Subsequently, the surface of CR is divided by a triangular surface, and the length of a triangular surface is less than λ / 10 . We utilize the method of moments (MoM) in Feko software to simulate CR, and the range of azimuth and pitch angle both are 0 ° , 90 ° where the angular interval is 1 ° .
Aiming at the fully polarimetric data, the coefficients of CR odd scattering and even scattering are obtained by CR decomposition, as shown in Figure A2 and Figure A3.
Figure A2 exhibits that how the coefficient of CR odd scattering changes in function of pitch and azimuth angles. The SR S S R 1 is the largest in the normal line of plane 1, corresponding to the areas φ = 0 90 ° and θ = 0 ° in Figure A2. The SR S S R 2 is the largest in the normal line of plane 2, corresponding to the area φ = 0 ° and θ = 9 0 ° . The SR S S R 3 is the largest in the normal line of plane 3, corresponding to the area φ = 9 0 ° and θ = 9 0 ° . The TR S T R is the largest in the LOS symmetry axis of CR, corresponding to the area φ = 4 5 ° and θ = 54 ° . For the odd scattering containing S S R 1 , S S R 2 , S S R 3 , and S T R at the same time, the odd scattering is a single-peak curve when the azimuth or pitch angle changes within 35 ° , which is the minimum of [ 35 ° , 45 ° , 90 ° ] , as displayed in Figure A2.
Figure A2. The odd scattering of CR.
Figure A2. The odd scattering of CR.
Remotesensing 14 03893 g0a2
Figure A3 presents that how the coefficient of CR even scattering changes with pitch and azimuth angles. The DR S D R 12 is the largest among the symmetrical parallel axes of planes 1 and 2, corresponding to the area φ = 0 ° and θ = 45 ° in Figure A3. The DR S D R 23 is the largest among the symmetrical parallel axes of planes 2 and 3, corresponding to the area φ = 4 5 ° and θ = 9 0 ° . The DR S D R 31 is the largest among the symmetrical parallel axis of planes 3 and 1, corresponding to the area φ = 9 0 ° and θ = 45 ° . For the even scattering containing S D R 12 , S D R 23 , and S D R 31 at the same time, the even scattering is a single-peak curve when the azimuth or pitch angle changes within 45 ° , which is a minimum of [ 45 ° , 90 ° ] , as shown in Figure A3.
Figure A3. The even scattering of CR.
Figure A3. The even scattering of CR.
Remotesensing 14 03893 g0a3
Combining the changing rule of odd scattering and even scattering, we find that the MSP coefficient of CR is a single-peak curve, when the azimuth or pitch angle of radar incident wave changes monotonically within 35 ° .

References

  1. Cerutti-Maori, D.; Klare, J.; Brenner, A.R.; Ende, J.H.G. Wide-area traffic monitoring with the SAR/GMTI system PAMIR. IEEE Trans. Geosci. Remote Sens. 2008, 46, 3019–3030. [Google Scholar] [CrossRef]
  2. Renga, A.; Graziano, M.D.; Moccia, A. Segmentation of marine SAR images by sublook analysis and application to sea traffic monitoring. IEEE Trans. Geosci. Remote Sens. 2019, 57, 1463–1477. [Google Scholar] [CrossRef]
  3. Gao, G.; Wang, X.; Lai, T. Detection of moving ships based on a combination of magnitude and phase in along-track interferometric SAR—Part II: Statistical modeling and CFAR detection. IEEE Trans. Geosci. Remote Sens. 2015, 53, 3582–3599. [Google Scholar] [CrossRef]
  4. Shi, L.; Li, P.X.; Yang, J.; Zhang, L.P.; Ding, X.L.; Zhao, L.L. Polarimetric SAR calibration and residual error estimation when corner reflectors are unavailable. IEEE Trans. Geosci. Remote Sens. 2020, 58, 4454–4471. [Google Scholar] [CrossRef]
  5. Luo, Y.; Guo, L.X.; Zuo, Y.C.; Liu, W. Time-domain scattering characteristics and jamming effectiveness in corner reflectors. IEEE Access 2020, 9, 15696–15707. [Google Scholar] [CrossRef]
  6. Baffelli, S.; Frey, O.; Werner, C.; Hajnsek, I. Polarimetric calibration of the Ku-band advanced polarimetric radar interferometer. IEEE Trans. Geosci. Remote Sens. 2018, 56, 2295–2311. [Google Scholar] [CrossRef]
  7. Sarabandi, K.; Chiu, T.C. Optimum corner reflectors for calibration of imaging radars. IEEE Trans. Antennas Propag. 1996, 44, 1348–1361. [Google Scholar] [CrossRef]
  8. Tao, Z.Y.; Yu, G.; Fu, Q. Discrimination method of ship and corner reflector based on SVM. In Proceedings of the IEEE Sixth International Conference on Instrumentation and Measurement, Harbin, China, 21–23 July 2016. [Google Scholar]
  9. Shi, Z.G.; Fu, Q. Study on the electromagnetic signature analysis and discrimination method of ship targets and outboard jamming. J. Hunan Univ. Sci. Technol. 2004, 19, 70–73. [Google Scholar]
  10. Zhang, R.W.; Li, Y.P.; Jiao, Y.F. Waveform selection of active radar seeker for anti-corner-reflector-interference. Electron. Opt. Control 2015, 22, 20–24. [Google Scholar] [CrossRef]
  11. Cui, K.B.; Wei, W.; Xi, C.; Yuan, N.C. A kind of method of anti-corner reflector interference for millimeter wave high resolution radar system. In Proceedings of the 2016 Progress in Electromagnetic Research Symposium, Shanghai, China, 8–11 August 2016. [Google Scholar]
  12. Yuan, W. Passive fake target recognition based on one-dimensional distance image sparse representation. Telecommun. Eng. 2018, 58, 798–804. [Google Scholar]
  13. Zhang, J.; Hu, S.L.; Yang, Q.; Fan, X.M. RCS statistical features and recognition model of air-floating corner reflector. Syst. Eng. Electron. 2019, 41, 780–786. [Google Scholar]
  14. Yuan, H.D.; Fu, X.J.; Zhao, C.X.; Xie, M.; Gao, X.Y. Ship and corner reflector identification based on extreme learning machine. In Proceedings of the IEEE International Conference on Signal, Information and Data Processing, Chongqing, China, 11–13 December 2019. [Google Scholar]
  15. Zhu, Z.Z.; Tang, G.F.; Cheng, Z.; Huang, P.K. Discrimination method of ship and corner reflector based on polarization decomposition. Shipboard Electron. Countermeas. 2010, 33, 15–21. [Google Scholar]
  16. Zhang, J.; Hu, S.L.; Fan, X.M.; He, J.B. Optimization of the position situation of the air-floating corner reflector based on HRRP and PA. Tactical Missile Technol. 2018, 3, 105–109. [Google Scholar]
  17. Zhang, J.; Hu, S.L.; Yang, Q.; Fan, X.M.; He, J.B. Optimization of position situation of air-floating corner reflectors based on similarity of RCS amplitude characteristic. J. Nav. Univ. Eng. 2019, 31, 32–36. [Google Scholar]
  18. Zhang, J.; Hu, S.L.; Wu, L.G.; Fan, X.M.; Yang, Q. Air-floating corner reflectors dilution jamming placement position. In Proceedings of the IEEE 8th Data Driven Control and Learning Systems Conference, Dali, China, 24–27 May 2019. [Google Scholar]
  19. Chen, V.C. Micro-doppler effect in radar: Phenomenon, model, and simulation. IEEE Trans. Aerosp. Electron. Syst. 2006, 42, 2–21. [Google Scholar] [CrossRef]
  20. Chen, V.C. High frequency scattering from trihedral corner reflectors and other benchmark targets SBR versus experiment. In The Micro-Doppler Effect in Radar; Artech House: Boston, MA, USA; London, UK, 2011; pp. 35–36. [Google Scholar]
  21. Graeme, E.S.; Karl, W. Radar micro-doppler signature classification using dynamic time warping. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 1078–1096. [Google Scholar]
  22. Wang, L.; Liu, H.W. Method for micro-motion target recognition and motion parameter extraction based on time-frequency analysis. J. Electron. Inf. Technol. 2010, 32, 1812–1817. [Google Scholar] [CrossRef]
  23. Huang, M.J.; Chen, J.J.; Zhao, H.Z.; Fu, Q. Micro-doppler modeling and simulation of corner reflector in sea surface. Syst. Eng. Electron. 2012, 34, 1781–1787. [Google Scholar]
  24. Huang, M.J.; Zhao, H.Z.; Fu, Q.; Feng, G.Y. A sea corner-reflector jamming identification method based on micro-doppler feature. J. Astronaut. 2012, 33, 1486–1491. [Google Scholar]
  25. Tao, Z.Y.; Fu, Q. Discrimination method of ship and corner reflector based on micro-doppler feature. In Proceedings of the 2016 International Conference on Computer Engineering, Information Science and Application Technology, Guilin, China, 24–25 September 2016. [Google Scholar]
  26. Zhu, H.; Wang, Q.P.; Pan, Y.J.; Tai, N.; Yuan, N.C. A sea corner-reflector jamming identification method based on time-frequency feature. In Proceedings of the IEEE International Conference on Signal Processing, Ningbo, China, 19–22 September 2015. [Google Scholar]
  27. Tu, J.H.; Tang, G.F.; Xiao, H.T.; Xu, C. A study of anti-corner reflector based on polarization decomposition. Radar Sci. Technol. 2009, 7, 85–90. [Google Scholar]
  28. Fang, M.J. Study on the Suppression of Compound Jamming Based on Polarization Coherent Radar; National University of Defense Technology: Changsha, China, 2013. [Google Scholar]
  29. Fang, M.J.; Zhu, Y.F.; Huang, M.J.; Fu, Q. Sea surface target polarization feature extraction based on modified odd-time and even-time scattering models. In Proceedings of the IEEE 2nd International Conference on Measurement, Information and Control, Harbin, China, 16–18 August 2013. [Google Scholar]
  30. Wang, W.Q. Radar ECM and ECCM Simulations; Beijing Institute of Technology: Beijing, China, 2016. [Google Scholar]
  31. Wang, W.Q.; Liu, K.; Zou, B.; Huang, Y.; Jiang, L.Z.; Hu, X. An Anti-Interference Method Based on Polarized Radar Corner Reflector. China Patent No. CN201710712673.3, 18 August 2017. [Google Scholar]
  32. Liang, Z.R.; Wang, Y.; Zhang, X.F.; Xie, M.; Fu, X.J. Identification of ship and corner reflector in sea clutter environment. In Proceedings of the IEEE 15th International Conference on Signal Processing, Beijing, China, 6–9 December 2020. [Google Scholar]
  33. Li, B.D.; Zhang, Z.J.; Jiang, J.; Xu, Y.Z.; Zhang, T.J.; Shi, J.N.; Chen, Y. An Interference Countermeasure Method for Sea Surface Corner Reflector Assisted by Polarization Characteristics. China Patent No. 201911227720.0, 4 December 2019. [Google Scholar]
  34. Knoit, E. RCS reduction of dihedral corners. IEEE Trans. Antennas Propag. 1977, 25, 406–409. [Google Scholar] [CrossRef]
  35. Griesser, T.; Balanis, C. Backscatter analysis of dihedral corner reflectors using physical optics and the physical theory of diffraction. IEEE Trans. Antennas Propag. 1987, 35, 1137–1147. [Google Scholar] [CrossRef]
  36. Griesser, T.; Balanis, C. Dihedral corner reflectors backscatter using higher order reflections and differactions. IEEE Trans. Antennas Propag. 1987, 35, 1235–1247. [Google Scholar] [CrossRef]
  37. Baldauf, J.; Lee, S.W.; Lin, L.; Jeng, S.K.; Scarborough, S.M.; Yu, C.L. High frequency scattering from trihedral corner reflectors and other benchmark targets SBR versus experiment. IEEE Trans. Antennas Propag. 1991, 39, 1345–1351. [Google Scholar] [CrossRef]
  38. Jackson, J.A.; Rigling, B.D.; Moses, R.L. Canonical scattering feature models for 3D and bistatic SAR. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 525–541. [Google Scholar] [CrossRef]
  39. Kubicke, G.; Bourlier, C.; Saillard, J. High-Frequency bistatic scattering by depolarizing, nearly omnidirectional reflectors higher order polyhedral reflectors. IEEE Trans. Antennas Propag. 2008, 56, 3029–3035. [Google Scholar] [CrossRef]
  40. Jackson, J.A. Analytic physical optics solution for bistatic, 3D scattering from a dihedral corner reflector. IEEE Trans. Antennas Propag. Lett. 2012, 60, 1486–1495. [Google Scholar] [CrossRef]
  41. He, Y.; He, H.; Hu, C.; Yin, J.; Yang, J. Polarization analysis of trihedral corner reflector with high-frequency approximation. IEEE Trans. Antennas Propag. 2022. ahead of print. [Google Scholar] [CrossRef]
  42. Zhang, J.; Hu, S.L.; Yang, Q.; Wu, L.G. Study on composite electromagnetic scattering characteristics of air-floating corner reflector and chaff. In Proceedings of the 12th International Conference on Intelligent Computation Technology and Automation, Xiangtan, China, 26–27 October 2019. [Google Scholar]
  43. Wu, J.; Khamas, S.K.; Cook, G.G. Analysis of a conformal archimedean spiral antenna printed within layered dielectric cylindrical media using the method of moments. IEEE Trans. Antennas Propag. 2012, 60, 3967–3971. [Google Scholar] [CrossRef]
  44. Sun, Q.; Ren, Q.; Zhan, Q.; Liu, Q.H. 3-D domain decomposition based hybrid finite-difference time-domain/finite-element time-domain method with nonconformal meshes. IEEE Trans. Microw. Theory Tech. 2017, 65, 3682–3688. [Google Scholar] [CrossRef]
  45. Feray, H.; Levent, S.; Pyotr, Y.U. On the modified theory of physical optics. IEEE Trans. Antennas Propag. 2013, 61, 6115–6119. [Google Scholar]
  46. Silver, S. Microwave Antenna Theory and Design; McGraw-Hill Book Company: New York, NY, USA, 1949. [Google Scholar]
  47. Gordo, W.B. Far-field approximations to the kirchoff-helmholtz resentations of scattered fields. IEEE Trans. Antennas Propag. 1975, 23, 590–592. [Google Scholar] [CrossRef]
  48. Ling, H.; Chou, R.C.; Lee, S.W. Shooting and bouncing rays: Calculating the RCS of an arbitrarily shaped cavity. IEEE Trans. Antennas Propag. 1989, 37, 194–205. [Google Scholar] [CrossRef]
  49. Krogager, E.; Czyz, Z.H. Properties of the sphere, diplance, helix decomposition. In Proceedings of the 3rd International Workshop Radar Polarization, Nantes, France, 21–23 March 1995. [Google Scholar]
  50. Lee, J.S.; Pottier, E. Polarimetric Radar Imaging: From Basic to Applications; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  51. Huynen, J.R. Phenomenological Theory of Radar Targets. Ph.D. Thesis, Technical University, Delft, The Netherlands, 1970. [Google Scholar]
  52. Yang, J.; Yamaguchi, Y.; Yamada, H. On modified Huynen’s decomposition of a Kennaugh matrix. IEEE Geosci. Remote Sens. Lett. 2006, 3, 369–372. [Google Scholar] [CrossRef]
  53. Balanis, C.; Polka, L.A. High-Frequency Techniques for RCS Prediction of Plate Geometries; Arizona State University Technical Report; NASA: Washington, DC, USA, 1993. [Google Scholar]
  54. Balanis, C.; Polycarpou, A.C. A physical optics/equivalent currents model for the RCS of trihedral corner reflectors. In Its RCS Analysis of Plate Geometries; NASA: Washington, DC, USA, 1993; Parts 1 and 2; pp. 22–41. [Google Scholar]
Figure 1. The relationship between the incident wave, reflected wave, and target.
Figure 1. The relationship between the incident wave, reflected wave, and target.
Remotesensing 14 03893 g001
Figure 2. Double scattering of triangular facets.
Figure 2. Double scattering of triangular facets.
Remotesensing 14 03893 g002
Figure 3. Triple scattering of triangular facets.
Figure 3. Triple scattering of triangular facets.
Remotesensing 14 03893 g003
Figure 4. The CR decomposition.
Figure 4. The CR decomposition.
Remotesensing 14 03893 g004
Figure 5. The spatial relationship between radar and target.
Figure 5. The spatial relationship between radar and target.
Remotesensing 14 03893 g005
Figure 6. The trajectory of radar terminal guidance.
Figure 6. The trajectory of radar terminal guidance.
Remotesensing 14 03893 g006
Figure 7. The angular variation of radar terminal guidance.
Figure 7. The angular variation of radar terminal guidance.
Remotesensing 14 03893 g007
Figure 8. The CR discrimination method based on TSP joint domains.
Figure 8. The CR discrimination method based on TSP joint domains.
Remotesensing 14 03893 g008
Figure 9. The changes of pitch and azimuth angle.
Figure 9. The changes of pitch and azimuth angle.
Remotesensing 14 03893 g009
Figure 10. The amplitude changes of odd scattering, even scattering, and MSP.
Figure 10. The amplitude changes of odd scattering, even scattering, and MSP.
Remotesensing 14 03893 g010
Figure 11. The model of a ship.
Figure 11. The model of a ship.
Remotesensing 14 03893 g011
Figure 12. The CR-MSE of the CR and ship.
Figure 12. The CR-MSE of the CR and ship.
Remotesensing 14 03893 g012
Figure 13. The discrimination probability and false alarm probability.
Figure 13. The discrimination probability and false alarm probability.
Remotesensing 14 03893 g013
Figure 14. Probability of false alarm under different discrimination times.
Figure 14. Probability of false alarm under different discrimination times.
Remotesensing 14 03893 g014
Figure 15. The discrimination probability under different SCRs.
Figure 15. The discrimination probability under different SCRs.
Remotesensing 14 03893 g015
Figure 16. The false alarm probability under different SCRs.
Figure 16. The false alarm probability under different SCRs.
Remotesensing 14 03893 g016
Figure 17. The errors of different decompositions.
Figure 17. The errors of different decompositions.
Remotesensing 14 03893 g017
Figure 18. The HRRP of a ship.
Figure 18. The HRRP of a ship.
Remotesensing 14 03893 g018
Figure 19. The HRRP of array CRs.
Figure 19. The HRRP of array CRs.
Remotesensing 14 03893 g019
Table 1. The number of MSP changes.
Table 1. The number of MSP changes.
P change 012
Number960737914
Table 2. The direction parameters of the radar and ship.
Table 2. The direction parameters of the radar and ship.
TypeThe Initial LOSThe Velocity of RadarThe Velocity of Ship
Azimuth φ L = random φ V = random φ T = random
Pitch θ L = 45 ° θ V = random θ T = 0 °
Table 3. The comparison of three methods.
Table 3. The comparison of three methods.
ProbabilityHRRP MethodPolarization MethodTSP Method
Discrimination 91.5 % 92.3 % 94.7 %
False alarm 14.4 % 8.8 % 4.1 %
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

He, Y.; Yang, H.; He, H.; Yin, J.; Yang, J. A Ship Discrimination Method Based on High-Frequency Electromagnetic Theory. Remote Sens. 2022, 14, 3893. https://doi.org/10.3390/rs14163893

AMA Style

He Y, Yang H, He H, Yin J, Yang J. A Ship Discrimination Method Based on High-Frequency Electromagnetic Theory. Remote Sensing. 2022; 14(16):3893. https://doi.org/10.3390/rs14163893

Chicago/Turabian Style

He, Yaomin, Huizhang Yang, Huafeng He, Junjun Yin, and Jian Yang. 2022. "A Ship Discrimination Method Based on High-Frequency Electromagnetic Theory" Remote Sensing 14, no. 16: 3893. https://doi.org/10.3390/rs14163893

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

Article Metrics

Back to TopTop