Next Article in Journal
Characteristic Scales of Tropical Convection Based on the Japanese Advanced Himawari-8 Imager Observations
Previous Article in Journal
Robust Object Categorization and Scene Classification over Remote Sensing Images via Features Fusion and Fully Convolutional Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Polarimetric Calibration Technique for a Fully Polarimetric Entomological Radar Based on Antenna Rotation

1
School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China
2
Advanced Technology Research Institute, Beijing Institute of Technology, Jinan 250300, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(7), 1551; https://doi.org/10.3390/rs14071551
Submission received: 4 February 2022 / Revised: 13 March 2022 / Accepted: 22 March 2022 / Published: 23 March 2022

Abstract

:
For entomological radar, the polarization information of the target is usually used to estimate the biological parameters, such as orientation, body length, and mass, of the insect. Thus, the accuracy of polarization measurement directly affects the performance of the biological parameters’ estimation. The polarization measurement error is mainly caused by the imbalance of amplitude and phase between two polarization channels and the crosstalk of the dual-polarization antenna. In order to obtain the correct polarization information of the target, the polarimetric calibration of the entomological radar is required. This paper proposes a new polarimetric calibration technique based on antenna rotation, which does not require the calibrator to have a specific polarization scattering matrix (PSM). Compared with the currently existing calibration techniques, no prior knowledge of the calibrator PSM is required (in fact, any fixed-point target can be used as a calibrator); thus, the errors introduced by the mechanical process can be avoided. Simulations and data measured by radar verify the effectiveness of the method. This method has the potential to be extended to other fully polarimetric radar systems in the future, such as fully polarimetric weather radar, fully polarimetric synthetic aperture radar (SAR), and so on.

1. Introduction

Entomological radar is of great importance in the application field of insect migration monitoring [1,2]. The traditional scanning entomological radar appeared in 1968. The research team led by Schaefer (1969) at the Center for Anti-Locust Research used a specially modified entomological radar for locust migration research in the south of the Saharato. This radar can obtain the number of migrating insects, and only the insect echo amplitude information is used for detection. In 1991, the Natural Resources Institute (NRI) and Commonwealth Scientific and Industrial Research Organisation (CSIRO) jointly developed a zenith-pointing linear-polarised conical-scan (ZLC) configuration for vertical observation radar (VLR), and successfully measured the body length, orientation, displacement direction, and flight speed of flying individuals in Australia [3]. The rotating linear polarization antenna was used to obtain the polarization information of insects, which helps the radar to measure insect body mass and orientation [3,4,5]. In 2016, the fully polarimetric entomological radar developed by the Beijing Institute of Technology can obtain insect targets PSM instantaneously [6], further improving the measurement capability of entomological radar [7,8,9].
The measurement of target polarization information is particularly important for an entomological radar system. Polarization information is used to estimate orientation, body length, and mass of the insect in the entomological radar; thus, the polarization measurement accuracy directly affects the estimation accuracy of insect parameters [10]. Studies have shown that when the polarization direction is parallel to the insect body axis, the echo power reaches the largest value. This conclusion is used as the basic assumption to calculate the orientation of insects [10]. Insect body axis direction can be obtained from the main eigenvector of the insect PSM [11,12]. The insect mass estimation is based on the monotonous mapping relationship between the insect mass and radar cross section (RCS) [13,14]. The estimations of insect body length and body width are based on two invariant target parameters calculated from the insect targets PSM [13]. One represents the target RCS when the polarization direction is perpendicular to the insect body axis, and the other represents the RCS product when the polarization direction is perpendicular to the insect body axis. Therefore, the measurement accuracy of the target PSM affects the estimation accuracy of body length, body mass, and orientation. PSM measurement accuracy is mainly affected by the polarization errors of radar system, including the stability and consistency of amplitude and phase between the horizontal polarization and vertical polarization channels, polarization isolation of dual-polarization antenna, etc. In order to ensure the measurement accuracy of insect PSM, the polarization calibration of an entomological radar is required [15].
The traditional polarization calibration is divided into active and passive techniques. Active polarization calibration refers to the use of a polarimetric active radar calibrator (PARC) to receive, amplify, delay, and transmit the radar radiation signal [16]. The polarization of the receiving antenna and transmitting antenna of PARC can be selected independently, and any PSM can be obtained by combing the polarization directions of the two antennas [17]. Passive polarization calibration refers to the radar calibration that measures the passive reflectors with a specific PSM, including metal sphere, corner reflector, metal wire, etc. [18,19]. Since the passive technique has the advantages of lower cost and simple maintenance, it is suitable for entomological radars in the field. The traditional passive calibration techniques include the isolated antenna calibration technique (IACT) [18], three targets calibration technique (TTCT) [18], single target calibration technique (STCT) [20], generalized calibration technique (GCT) [18,21], and so on. These techniques realize the polarimetric calibration by measuring one to three passive standard calibrators and solving the polarimetric errors [21]. These traditional polarization calibration techniques require high precision of calibrator machining and calibrator posture during measurement (for example, dihedral of a specific angle). The problems of the traditional polarization calibration in the field include difficult calibrator machining and storage, and difficult operation of the calibration process; thus, the real calibration error is large and difficult to apply in practice. Among these techniques, STCT only needs to measure a metal sphere, so it is more practical. However, the machining precision and damage during storage limit the calibration accuracy. To avoid the above problems of passive polarization calibration, this paper proposes a new polarimetric calibration technique based on antenna rotation (PCTAR), which calibrates the systematic polarization error by rotating the antenna and measuring the echoes of the calibrator in different polarization directions. The technique has no requirements for the posture and PSM form of the calibrator, and any target can be used as the calibrator. It is easy to implement in actual field calibration. PCTAR uses the polarization error model of STCT. By rotating the antenna of the radar, a series of PSMs of the calibrator with systematic polarization errors in different directions are obtained. Next, based on the polarization error model and the measured PSMs, the system polarization errors are calculated. Finally, the system polarization errors are calibrated based on the polarization error model. In this paper, the theoretical analysis and simulation of the PCTAR technique are carried out to verify the feasibility. The performance for practical application is also tested based on our self-developed, fully polarimetric radar system.
The subsequent arrangement of this paper is as follows: Section 2 introduces the polarization measurement error model of the fully polarimetric radar; Section 3 gives the details of the PCTAR polarization calibration technique; In Section 4, the detailed performance simulations of PCTAR under different conditions are carried out, and the performances of PCTAR and STCT are compared; In Section 5, field experiments were conducted based on a fully polarimetric entomological radar to verify the effectiveness of the proposed method; and Section 6 and Section 7 are the discussion and conclusion of the paper, respectively.
Mathematical Notations: In this paper, the scalars are indicated with lower case, non-bold letters. Vectors are indicated with bold, upper-case letters. Matrices are indicated with bold, capital letters. ( ) T and ( ) 1 denote the transpose, and the inverse transpose of a matrix. The symbol | | denotes the modulus of a number.

2. Polarization Error Model of the Fully Polarimetric Entomological Radar System

A fully polarimetric entomological radar developed by our team will be taken as an example to introduce the polarimetric error model. The fully polarimetric entomological radar system is a two-dimensional scanning radar. It is capable of transmitting and receiving two polarization orthogonal signals (horizontal polarization, H, and vertical polarization, V) at the same time. The PSM of the target in the coordinate system can be obtained by constructing the coordinate system with mutually perpendicular H- and V-polarization directions. The system block diagram is shown in Figure 1, including a dual-channel RF transceiver module with a bandwidth of 800 MHz and the center frequency of 16.2 GHz, a dual-polarized antenna, and a signal processor. A high resolution range profile (HRRP) for target detection is generated by transmitting, receiving, and processing the step-frequency chirp signals. In a previous paper [15], the internal design block diagram and signal processing of the radar was described in detail.
The signal processor generates two sets of orthogonal polarization baseband signals, each of which consists of a number of stepped frequency chirp pulses. The orthogonal isolation is realized by reverse frequency modulation and different phase encoding among pulses in the stepped frequency cycle [15]. The baseband signals are converted to RF signals by up-conversion, amplifying, filtering, and so on. Next, the signals pass through the circulators. Finally, the signals of two channels are transformed into H-polarization and V-polarization signals through the dual-polarization antenna, and then radiated into space. Due to reverse frequency modulation and phase encoding, the two signals are isolated from each other and can be restored separately by demodulation. Due to the depolarization effect of the target, the H-polarization signal is partly reflected as an H-polarization signal (HH) and partly as a V-polarization signal (VH). Similarly, V-polarization signals are also partly reflected as H-polarization signals (HV) and partly as V-polarization signals (VV). After the processes of antenna receiving, transceiver switching, down-conversion, amplification, and filtering, each signal is demodulated by two kinds of coding and matched filtering, so as to get four kinds of polarization signals (HH, HV, VH, and VV), which are, namely, the elements of PSM.
For an ideal, fully polarimetric radar, the target PSM can be obtained based on the above process. However, there are many unsatisfied factors due to the complex multiple polarization channels of the radar system. First, it is impossible to make a dual-polarization antenna with the ideal polarization isolation. The coupling exists between two polarization ports of the antenna. Second, the two polarization signals pass through two independent frequency conversion amplification circuits, which contain some active devices, such as amplification circuits, filters, mixers, etc. These devices exhibit different amplitude and phase characteristics at different temperatures, which cannot be completely consistent. The above undesirable factors will seriously distort the measured PSM, making it impossible to use the effective polarization information of the target. In order to obtain the correct PSM of a target, it is necessary to establish the systematic polarization error model of the measured PSM, and then solve and compensate for the systematic error based on the model.
Sarabandi and colleagues described the fully polarimetric antenna as a four-port network and establishes a complete systematic polarization error model for a fully polarimetric radar system with a single antenna and dual-polarization channels [18]. This paper directly follows this model:
M = RCS C T T + N
In Equation (1), M refers to the measured (uncalibrated) PSM of the target with noise; R and T refer to the error matrices of the receiving and transmitting channels, respectively, which can be written as:
R = g r h 0 0 g r v
T = g t h 0 0 g t v
where g r h and g r v represent the amplitude and phase errors (gains) of the H-polarization and V-polarization receiving channels, respectively; g t h and g t v represent that of the transmitting channels, respectively. S represents the true value of target PSM, which is written as:
S = s h h s h v s v h s v v
where the subscript i, j of s i j ( i , j h , v ) represent the polarization forms of receiving and transmitting, respectively. According to the reciprocity principle of radar polarization with single station and single antenna, s h v = s v h . C refers to antenna crosstalk matrix, which is written as:
C = 1 C 1 C 2 1
where C 1 and C 2 represent antenna leakage between polarization channels, respectively, and C 1 C 2 1 . C can be assumed to be a symmetric matrix [18]. N = N 11 N 12 N 21 N 22 , representing a complex, white Gaussian noise matrix. During polarization calibration, the effects of additive white noise are usually ignored, so the model is simplified as:
M = R C S C T = g r h 0 0 g r v 1 C 1 C 1 1 s h h s h v s h v s v v 1 C 1 C 1 1 g t h 0 0 g t v = g r h g t h 1 0 0 G r 1 C 1 C 1 1 s h h s h v s h v s v v 1 C 1 C 1 1 1 0 0 G t
where G r = g r v g r h and G t = g t v g t h .
Since the absolute amplitude of the target (RCS) is obtained through RCS calibration, polarization calibration only needs to realize the relative value calibration between different polarization channels of the target. That is, obtaining the amplitude and phase information of each polarization channel relative to the HH-polarization channel is enough. Therefore, the absolute amplitude item g r h g t h can be discarded, and Equation (6) can be further simplified as:
M = 1 0 0 G r 1 C 1 C 1 1 s h h s h v s h v s v v 1 C 1 C 1 1 1 0 0 G t
Polarization calibration is the process of solving the three systematic polarization error matrices of R = 1 0 0 G r , T = 1 0 0 G t , and C. Next, the real target PSM can be obtained by matrix inversion and multiplication:
S = C 1 R 1 M T 1 C 1
So far, the polarization measurement model of a fully polarimetric radar has been established. The following describes how to solve the three unknown matrices of R , T , and C by a polarimetric calibration technique based on antenna rotation designed in this paper.

3. Polarimetric Calibration Technique Based on Antenna Rotation (PCTAR)

For the proposed method, the target used for calibration is not required to be a regular calibrator with a known PSM. It can be an arbitrary, single scatter point target, such as a metal sphere, metal wire, irregular metal object, etc. During measurement, it is necessary to place the target in the center of the radar beam, and rotate the radar beam around the axis of the beam center for at least one rotation period (0°–360°) (Figure 2a). If the antenna cannot rotate, the measurement can also be conducted by rotating the target around the axis of the radar beam center for at least one rotation period (Figure 2b). Thereby, a series of target PSMs at different angles (at least one rotation period) are obtained.
The steps and calculation methods of PCTAR polarization calibration technique are described in detail below. When the rotation angle θ is 0°, taking the H-polarization direction as the reference direction, the target PSM can be expressed as:
S 0 = s h h 0 s h v 0 s h v 0 s v v 0
The direction of H-polarization (V-polarization) rotates with the rotation of the radar beam, which is equivalent to the rotation of the polarization coordinate system. The change of the target PSM with the rotation angle θ can be expressed by the rotation matrix [18]:
S θ = s h h θ s h v θ s v h θ s v v θ = cos θ sin θ sin θ cos θ s h h 0 s h v 0 s h v 0 s v v 0 cos θ sin θ sin θ cos θ
The cos θ sin θ sin θ cos θ is the rotation matrix, and
s h h θ = s h h 0 cos 2 ( θ ) + 2 s h v 0 cos ( θ ) sin ( θ ) + s v v 0 sin 2 ( θ ) s h v θ = s h h 0 cos ( θ ) sin ( θ ) + s h v 0 cos 2 ( θ ) s h v 0 sin 2 ( θ ) + s v v 0 cos ( θ ) sin ( θ ) s v h θ = s h h 0 cos ( θ ) sin ( θ ) s h v 0 sin 2 ( θ ) + s h v 0 cos 2 ( θ ) + s v v 0 cos ( θ ) sin ( θ ) s v v θ = s h h 0 sin 2 ( θ ) 2 s h v 0 cos ( θ ) sin ( θ ) + s v v 0 cos 2 ( θ )
Equation (10) represents the real target PSM. The measured PSM can be obtained by substituting Equation (10) into Equation (7), marked as M t θ = M h h t θ M h v t θ M v h t θ M v v t θ , where
M h h t = s h h t + 2 C 1 s h v t + C 1 2 s v v t M h v t = G t ( 1 + C 1 2 ) s h v t + C 1 ( s h h t + s v v t ) M v h t = G r ( 1 + C 1 2 ) s h v t + C 1 ( s h h t + s v v t ) M v v t = G t G r ( s v v t + 2 C 1 s h v t + C 1 2 s h h t )
The average of each element of M t θ over different angles (0°~360°) can be represented as:
{ M h h t | 0 360 = s h h t | 0 360 + 2 C 1 s h v t | 0 360 + C 1 2 s v v t | 0 360 M h v t | 0 360 = G t [ ( 1 + C 1 2 ) s h v t | 0 360 + C 1 ( s h h t | 0 360 + s v v t | 0 360 ) ] M v h t | 0 360 = G r [ ( 1 + C 1 2 ) s h v t | 0 360 + C 1 ( s h h t | 0 360 + s v v t | 0 360 ) ] M v v t | 0 360 = G t G r ( s v v t | 0 360 + 2 C 1 s h v t | 0 360 + C 1 2 s h h t | 0 360 )
where, refers to the operation to take the average of the values in angle brackets. In the process of beam rotation, H-polarization and V-polarization can be regarded as two linear polarizations with an angle difference of 90°; thus, the sum of the measured signals of H-polarization and V-polarization in one period is equal, namely s h h | 0 360 = s v v | 0 360 .
By taking [Equation (13)-④]/[Equation (13)-①] and [Equation (13)-③]/[Equation (13)-②], the following equation can be obtained:
M v v t | 0 360 M h h t | 0 360 = G t G r M v h t | 0 360 M h v t | 0 360 = G r G t
According to Equation(14), two sets of solutions for G t and G r with opposite signs, can be calculated:
G t = M v v t | 0 360 M h v t | 0 360 M h h t | 0 360 M v h t | 0 360 G r = M v v t | 0 360 M v h t | 0 360 M h h t | 0 360 M h v t | 0 360
or
G t = M v v t | 0 360 M h v t | 0 360 M h h t | 0 360 M v h t | 0 360 G r = M v v t | 0 360 M v h t | 0 360 M h h t | 0 360 M h v t | 0 360
Next, we will introduce how to solve the cross-polarization coupling matrix parameter C 1 . Through G t G r [Equation (12)-①] + [Equation (12)-④], and G r [Equation (12)-②], it can be found that:
M h h t G t G r + M v v t = G t G r ( ( s h h t + s v v t ) ( 1 + C 1 2 ) + 4 C 1 s h v t ) M h v t G r = G t G r [ ( 1 + C 1 2 ) s h v t + C 1 ( s h h t + s v v t ) ]
Substitute Equation (11) into Equation (17), and we can obtain:
M h h t G t G r + M v v t = G t G r ( ( s h h t + s v v t ) ( 1 + C 1 2 ) + 4 C 1 [ 1 2 ( s v v t s h h t ) sin 2 θ + s h v t cos 2 θ ] ) M h v t G r = G t G r [ ( 1 + C 1 2 ) [ 1 2 ( s v v t s h h t ) sin 2 θ + s h v t cos 2 θ ] + C 1 ( s h h t + s v v t ) ]
As at least one rotation period was measured, the terms with cos 2 θ in Equation (18) can be eliminated by taking the average of the term in one period or its integral in multiple periods. The average of M h h t G t G r + M v v t and M h v t G r at each angle is defined as two constant items, u, and v, where
u = 1 n M h h t G t G r + M v v t v = 1 n M h v t G r
where n is the number of angles in one rotation period. By taking the average of both sides of equation in Equation (18), we can get:
u = G t ( s h h t + s v v t ) ( 1 + C 1 2 ) v = G t C 1 ( s h h t + s v v t )
Next, it can be found that
u v = 1 + C 1 2 C 1
We get:
C 1 = u 2 v + u 2 v 2 1
or:
C 1 = u 2 v u 2 v 2 1
Two solutions of C 1 can be obtained. Since C 1 represents the polarization isolation that 1 , the smaller one is usually the real solution.
The polarization calibration can be realized by taking the values of G t , G r , and C 1 into Equation (7):
S u = 1 C 1 C 1 1 1 1 0 0 G r 1 M u 1 0 0 G t 1 1 C 1 C 1 1 1
In Equation (24), M u is the measured value of the target PSM and S u is the truth value of relative PSM. The truth value of the PSM can be obtained by an additional traditional RCS calibration, which is not described here.
As two sets of solutions of G t and G r with opposite signs are obtained, two calibrated matrices can be obtained. The difference between these two matrices is that the signs of the cross-polarization elements are opposite. Similar to the traditional STCT polarization calibration technique, the true value can be determined by measuring the metal wire with an orientation of 45°. Two solutions are taken to calibrate the PSM of metal wire. If the orientation (the calculation method refers to [12]) calculated by the calibrated PSM is around 45°, this solution is the correct one. Otherwise, if the orientation is around 45°, the selected solution is wrong.

4. Simulation Analysis

In the derivation of the proposed method, it is assumed that the radar measurement is free of noise, and the calibrator is always located on the normal of the antenna beam during the rotation of the antenna. However, in the actual calibration process, noise interference is inevitable. In addition, the calibrator may deviate from the normal of the beam due to the interference of environmental factors and rotation control accuracy. In this section, firstly, the effect of noise on the proposed method is simulated. Second, the calibration error caused by beam pointing error is analyzed. Finally, the performances of STCT and PCTAR under difference noise levels are compared.
In this section, assumptions of radar system parameters in all simulations are as follows: the matrix for the amplitude and phase imbalances between the receiving channels is R = [ 1 0 0 1.2 e 1.57 j ] ; the matrix for the amplitude and phase imbalances between the transmitting channels is T = 1 0 0 1.1 e 1.05 j ; the isolation of the fully polarimetric antenna is −25 dB, that is, C = 1 0.0562 0.0562 1 ; and the antenna half-power beam width ( θ 3 dB ) is 1.5 . The assumptions fit the polarization errors of the self-developed radar system.

4.1. Effects of Noise on the PCTAR

The performance of the proposed method under different signal-to-noise ratios (SNRs) is analyzed in this section. The expectation (E) and standard deviation ( δ ) of the amplitude and phase errors of each element in the PSM are used to evaluate the performance of the PCTAR technique. The amplitude and phase errors are defined as the absolute value, Δ S u , of the amplitude and phase difference between the calibrated PSM, S e u , and the real PSM, S u , of the target, respectively.
The simulation takes an ideal metal sphere as the calibrator, and its theoretical PSM is S = 1 0 0 1 . The SNR changes from 10 dB to 40 dB. Since the power of the noise is constant during the rotation of the radar antenna, and the power of the target echo signal usually changes with the angle, the SNR usually changes with the rotation angle θ . For convenience, the SNR in the simulation is defined as the SNR of the HH channel at θ = 0 o . Depending on the different SNR, there are several groups of simulation conditions. The Monte Carlo simulation is used to generate 10 5 PSMs S u of insects randomly in each condition. By using Equation (13), a series of PSMs measurements at 0°~360° are obtained. After adding the complex, white Gaussian noise matrix N, the target PSMs M e | 0 360 with errors are obtained. Equations (15), (16), (18), (22), and (23) are used to calculate R e , T e , and C e , and the calibrated targets PSMs S e u are calculated based on Equation (24). S u and S e u are normalized by the element of HH polarization, and can then be used to calculate the amplitude and phase errors. The expectation, E, and standard deviation, δ , of the amplitude and phase errors of the three polarization channels (HV, VH, and VV) for 10 5 targets are calculated. Simulation results are shown in Figure 3.
Figure 3 shows the influence of SNR on the performance of PCTAR. It can be seen that both expectation, E, and standard deviation, δ , decrease with SNR, which means the calibration accuracy gradually increases with the increasing SNR. Since the HH channel signal is set as the reference, the VV channel error represents the co-polarization channel error in the simulation results. Cross-channel meets the consistency of reciprocity error accuracy, and, therefore, the amplitude and phase error of the HV or VH channel represents cross-polarization channel amplitude and phase error. It can be seen that the calibration precision of the co-polarization channel (VV) is higher than that of the cross-polarization channels (VH) under the same noise level. When the SNRs are larger than 15 dB, for all polarization channels, E a m p 0.007 dB ( δ a m p 0.6 dB), and E p h a 0.04 ° ( δ p h a 3.8 ° ). This means that the calibration is of high accuracy.

4.2. Effects of Beam Alignment Error on PCTAR

The PCTAR technique has no requirements for the prior knowledge of the calibrator PSM. During the calibration, we need to place the target in the center of the antenna beam, and rotate the antenna at least one rotation period. However, errors will be introduced if the target is not on the center of the antenna beam. There are two main factors leading to the offset to the center of antenna beam. One is that the rotation axis is not in alignment with the radar beam center due to the errors of mechanical process or an installation error. The other is that the actual position of the target is not on the rotation axis caused by the experimental error. The effects of these two factors on the measured PSM are analyzed in Appendix A. In the following part of this section, the effect of beam alignment error on the proposed method is simulated.
The ideal metal sphere is used as the calibrator in the simulation. The simulation parameters are set as follows: the misalignment angle between the rotation center axis and the antenna beam center axis is φ 2 [ 0 ° , 1.5 ° ] with the interval step of 0.1°; OA [ 0 m , 5 m ] with the interval step of 1 m. There are several groups of error conditions. In each condition, the Monte Carlo simulation method is used to randomly generate PSMs of 10 5 insects. The theoretical PSM of each insect is represented as S u . A series of PSMs measured at different angles (0°~360°) with a beam alignment error (represented as M e r r o r 0 ~ 360 ) are obtained by using Equations (6), (13), and (A6). The calibration matrices with errors R e , T e , and C e are calculated by using Equations (15), (16), (18), (22) and (23). The calibrated target PSM S e u is calculated by Equation (24). Next, the amplitude and phase errors, as well as their expectations (E) and standard deviations ( δ ), are calculated, respectively. The co-polarization (VV) and cross-polarization (HV) channel errors are shown in Figure 4 and Figure 5, respectively.
Figure 4 shows the simulation results of a co-polarization channel under the different distance errors, OA, and the angle errors, φ 2 . If φ 2 is not zero, as the OA increases, the E a m p V V , E p h a V V , δ a m p V V , and δ p h a V V increase accordingly. If OA is not zero, as φ 2 increases, the E a m p V V , E p h a V V , δ a m p V V , and δ p h a V V increase accordingly. When φ 2 1.5 and OA 5 m, the calibration error of the co-polarization channel can be obtained as: E a m p V V 1   dB , δ a m p V V 0.062   dB , E p h a V V 4.8 , and δ p h a V V 0.42 .
Figure 5 shows the simulation results of cross-polarization channels under the different distance errors, OA, and the angle errors, φ 2 . The relationships between OA and φ 2 and the calibration precision of the cross-polarization channel are similar to that of the co-polarization channel. When φ 2 1.5 and OA 5 m, the amplitude and phase errors of the cross-polarization channel can be obtained as: E a m p H V 0.5   dB ( δ a m p H V 0.038   dB ), and E p h a H V 2.5 ( δ p h a H V 0.25 ).
The simulation indicates that the beam alignment error deteriorates the calibration precision of PCTAR. Therefore, in the polarization calibration measurement, the angle between the antenna beam and the rotation axis should be reduced as much as possible, and the calibrator should be placed on the center of the beam as much as possible to improve the polarization calibration precision.

4.3. Comparison of PCTAR and STCT

The performances of PCTAR and STCT in difference noise levels are compared in this section. The simulation and evaluation methods are similar to that in Section 4.1. First, the influence of SNR on the STCT is analyzed, as shown in Figure 6. Next, the performances of PCTAR and STCT are compared in Figure 7.
It can be seen from Figure 6 that the relationships between the amplitude and phase errors and SNR of STCT are similar to that of PCTAR. When the SNR reaches 15 dB, the amplitude calibration errors of all polarization channels are E a m p 0.35   dB and δ a m p 5   dB ; the phase calibration errors of all polarization channels are E p h a 28 and δ p h a 60 . To compare the performance of PCTAR and STCT, the differences of the expectation, E, and standard deviation, δ , of the amplitude and phase errors at different SNRs and different polarization channels of STCT (Figure 6) and PCTAR (Figure 3) are calculated. The differences of the expectation and the standard deviation of the amplitude and phase errors are denoted as Δ E a m p , Δ δ a m p , Δ E p h a , and Δ δ p h a , respectively. If the differences are greater than 0, it means that the precision of PCTAR is better than that of STCT. Otherwise, it means that the precision of STCT is better than that of PCTAR. The comparison results are shown in Figure 7:
The accuracies of PCTAR are much higher than STCT at SNRs lower than 26 dB. When SNR > 26   dB , Δ E a m p and Δ E p h a are close to zero, which indicates that the error expectations of the two methods are close. However, the standard deviations Δ δ a m p and Δ δ p h a of PCTAR are smaller than that of STCT. This indicates that the PCTAR is more robust to the system noise than STCT. Therefore, PCTAR can achieve higher accuracy than STCT under different SNRs. The reason is that STCT uses only one PSM of the calibrator, while PCTAR uses a series of PSMs (0°~360°), which has an accumulation effect on the echo signal, and the SNR is improved.

5. Experimental Verification

In order to verify the feasibility of PCTAR, a practical experiment is carried out based on a fully polarimetric radar. The radar is a Ku-band, high-resolution, fully polarimetric entomological radar set up in Jiangcheng, Yunnan, China. The details of the experiments and results are shown as follows.

5.1. Experiment

The experiment was carried out in calm weather. Before the measurement, the radar was set on a flat place to ensure the base level. An inclination sensor was used to measure the inclination angle of the radar base at different directions, and the radar was adjusted until the inclination angles in all directions are smaller than 0.5 ° . Next, the antenna beam was adjusted to point vertically to the sky. Before the measurement, the radar worked for about thirty minutes so that it could warm up to bring the amplitude and phase characteristics of the radar to a thermal equilibrium state. After that, the calibrator was placed in the center of the radar beam. Two DJI M600 unmanned aerial vehicles (UAVs) were used to suspend the target at an altitude of 300 m above the radar. The two UAVs were 120 m away, connected by a polyethylene (PE) line with low reflectivity. The calibrator was suspended in the middle of the line. The position of the calibrator was adjusted to be on the axis of radar beam. The radar antenna rotated along the azimuth angle for one cycle at a constant speed of 1 ° / s , and the PSMs at different angles were recorded. A schematic diagram of the experiment is shown in Figure 8. The experiment scenario is shown in Figure 9.
The calibration experiment flow is shown as follows:
1. Start;
2. Set up the radar;
3. Warm up the radar;
4. Place the calibrator on the axis of the radar beam left;
5. Activate radar transmission and record data;
6. Rotate the radar antenna for one cycle;
7. Calculate the error matrices C, R, and T;
8. Shut down the radar;
9. End.
Four targets were measured in the experiment (Figure 10):
  • Target I: A metal sphere with a diameter of 122 mm used as a calibrator for STCT;
  • Target II: A 120 mm wrench used as a calibrator for PCTAR;
  • Target III: A metal sphere with a diameter of 60 mm used as a target;
  • Target IV: An insect (bollworm) with 15 mm length used as a target.

5.2. Experimental Results

Target I, the sphere with a diameter of 122 mm, was used as the calibrator for STCT. The error matrices obtained by STCT are as follows: G r | = 0.8276 e 0.3566 i ; G t | = 1.1564 e 2.3007 i ; and C 1 | = 0.0065 e 0.1365 i . Target II, the 120 mm wrench, was used as a calibrator for PCTAR. The error matrices obtained by PCTAR are as follows: G r = 0.7644 e 0.4103 i ; G t = 1.1835 e 2.2619 i ; and C 1 = 0.0064 e 0.0814 i .
Target III, the sphere with a diameter of 60 mm, was used to validate the calibration performance of these two methods. For this purpose, the polarization pattern is introduced here, which represents the echo signal power at different polarization directions σ θ , and can be calculated from the PSM as [13]:
σ θ = h T θ S h θ
where h θ represents the effective antenna length, and can be written as:
h θ = cos θ sin θ
The ideal polarization pattern of a sphere is a standard circle. Thus, the shape of the polarization pattern of the sphere can be used to qualitatively evaluate the calibration accuracy.
The polarization patterns of Target III without calibration (blue), calibrated with STCT (red), and calibrated with PCTAR (black) are shown in Figure 11, respectively. It can be seen that the polarization patterns after calibration are close to a standard circle, which means the calibrations of the two methods are effective. The polarization pattern calibrated with PCTAR is closer to a standard circle that that with STCT, which means the accuracy of PCTAR is higher than that of STCT in the experiment.
The polarization patterns of Target IV, an insect, without calibration (blue), calibrated with STCT (red), and calibrated with PCTAR (black) are shown in Figure 12, respectively. As the true polarization pattern of the insect is unknown, the polarization pattern calibrated with STCT can be treated as the true value because the STCT has been fully verified. In fact, this is a typical shape of an insect polarization pattern [13]. It can be seen that the polarization pattern calibrated with PCTAR is very close to that with STCT. This means that the calibration of PCTAR highly accurate, and can be used for an entomological radar.

6. Discussion

The traditional polarization calibration method requires a standard calibrator with a known theoretical PSM. The machining accuracy and posture errors of the calibrator will affect the calibration precision. This paper proposes a polarization calibration technique based on antenna rotation, which has no assumption for the target PSM. Therefore, this method avoids the errors introduced by calibrator. However, during the measurement, if the beam center is not aligned with the calibrator or the beam center is inconsistent with the rotation axis, PSM measurement error will be introduced and the calibration accuracy will be affected.
The effectiveness of the method is verified by simulation and an experiment. It should be pointed out that, although the calibrator used in the simulation is metal sphere, this method has no requirements for the calibrator form, and any form of calibrator can be used. For example, a wrench was used as the calibrator in the experiment, which realized a polarization calibration of high quality. Nevertheless, in order to reduce the influence of noise, it is recommended to use a target with a strong radar echo as the calibrator. Furthermore, it should be noted that, before the polarization calibration of the radar, it is necessary to calculate whether the radar works in the linear gain region for the calibrator and avoid errors caused by the saturation of the radar receiving channel.
The proposed method requires the radar to rotate along the beam center axis for at least one cycle. For an entomological radar with rotation function, the equipment cost and device cost of polarization calibration can be reduced. For the fully polarimetric radar without rotation function, it can be realized by rotating the calibrator, but this is usually difficult to realize.
The proposed method provides a novel idea for the polarization calibration of a fully polarimetric radar. As an extended application, the polarization calibration and correction of a radar system may be realized by using nonstandard calibration in nature, which will be a good future research topic.

7. Summary

The fully polarimetric entomological radar used to monitor the biological characteristics of airborne migratory insects has a requirement for the accuracy of the polarization measurement. A polarimetric calibration technique based on antenna rotation is designed in this paper, aiming to meet the requirement of high polarization calibration precision and solve the problem of inconvenient experimental operations for measuring the specific posture of a calibrator in the field environment. Based on this new method, the highly accurate polarization measurement of the fully polarimetric entomological radar is realized. Compared with other passive calibration techniques, there is no assumption of the calibrator PSM, and, therefore, the errors of the mechanical process introduced by the hypothetical calibrator PSM are avoided. In this paper, the fully polarimetric system model is first described, and the PCTAR method is proposed based on the model. The performance of the proposed method at different noise levels is analyzed by simulation. The influence on the calibration accuracy due to the offset of calibrator to the beam axis is analyzed based on theoretical derivation and simulation. In addition, the performances of the proposed PCTAR and the traditional STCT is compared under different SNRs. The simulation results show that the accuracy of the PCTAR technique is much higher than STCT. Finally, the effectiveness of PCTAR is verified based on a field experimental conducted with a fully polarimetric entomological radar. The experiment indicates that this method can be used for the real fully polarimetric radar.

Author Contributions

T.Y., M.L. and W.L. developed the method and wrote the manuscript. T.Y. and H.M. designed and carried out the experiments. T.Y., R.W., C.H. and T.L. reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Special Fund for Research on National Major Research Instruments under Grant 31727901; Major Scientific and Technological Innovation Project (2020CXGC010802); The National Natural Science Foundation of China (Grant No. 62001021).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Effects of Beam Alignment Error on the Measured PSM

This section will theoretically analyze the PSM measurement errors of the calibrator in rotation caused by the misalignment between the rotation axis and the radar beam center, and the position offset of the target relative to the rotation axis.
The radar beam model is shown in Figure A1. Figure A1a is the side view of the radar beam model, where the radar position is set as the origin point O . The calibrator position is A. The height is H. The rotation center axis is O O . O is the intersection point between the horizontal plane of the calibrator and the rotation center axis. O B is the beam center axis. B is the intersection point between the beam center axis and the horizontal plane of calibrator. The angle between the rotation center axis O O and the antenna beam center axis O B is φ 2 . Let O O A = φ 1 and A O B = φ 3 . During the measurement, point B rotates around O. AOB is θ . φ 2 is the misalignment angle between the rotation axis and the radar beam center, and OA is the position offset of the calibrator relative to the rotation axis. Figure A1b is the top view of the radar beam model. The vector OA is defined as the x-axis. The direction perpendicular to the x-axis is defined as the y-axis. O is the origin. In this coordinate system, the beam center rotates around point O in the measurement.
Figure A1. Diagram of the radar beam model: (a) side view; (b) top view.
Figure A1. Diagram of the radar beam model: (a) side view; (b) top view.
Remotesensing 14 01551 g0a1
It can be easily found that OA = H tan φ 1 and OB = H tan φ 2 . The vector AB can be expressed as:
AB = ( OB cos θ OA , OB sin θ ) ; ( θ 0 ~ 360 ° )
Take OA and OB into Equation (A1), and we can get:
AB = H tan φ 2 cos θ H tan φ 1 , H tan φ 2 sin θ ; ( φ 1 0 ° φ 2 0 ° θ 0 ~ 360 ° )
Then,
A B = H tan φ 2 cos θ H tan φ 1 2 + H tan φ 2 sin θ 2
Since H is generally larger than the blind range of the radar, it is usually much larger than AB . Thus, the approximate expression of the angle between the calibrator position and the radar beam center axis, φ 3 , can be expressed as:
φ 3 atan AB H = atan H tan φ 2 cos θ H tan φ 1 2 + H tan φ 2 sin θ 2 H = atan tan φ 2 cos θ tan φ 1 2 + tan φ 2 sin θ 2 ; ( θ 0 ~ 360 ° )
Assuming that the antenna pattern is of the Gaussian function, the antenna gain can be written with φ 3 as:
f G φ 3 = exp 4 l n 2 φ 3 2 θ 3   dB 2
f G φ 3 represents the gain loss caused by φ 3 . According to the radar equation, the echo power loss cause by φ 3 is f G 2 φ 3 , the square of f G φ 3 . In the rotation measurement of the calibrator, the measured PSM of the calibrator can be written as:
M e r r o r θ = f G 2 φ 3 M θ
where M e r r o r θ represents the measured PSM with an error at rotation angle θ ; M θ represent the true value of the PSM.
It can be seen from Equation (A6) that the test error of PCTAR is mainly caused by the nonzero angle φ 3 , and that the error angle φ 3 is determined by the beam pointing angle error φ 2 , the calibrator height H, and the error of target position from beam center OA. With an increase of OA and a decrease of H, the precision of PCTAR gradually decreases.

References

  1. Wang, R.; Mao, H.; Hu, C.; Zeng, T.; Long, T. Joint association and registration in a multiradar system for migratory insect track observation. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 4028–4043. [Google Scholar] [CrossRef]
  2. Chapman, J.W.; Drake, V.A.; Reynolds, D.R. Recent insights from radar studies of insect flight. Ann. Rev. Entomol. 2011, 56, 337–356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Smith, A.D.; Riley, J.R.; Gregory, R.D. A method for routine monitoring of the aerial migration of insects by using a vertical-looking radar. Philos. Trans. Biol. Sci. 1993, 340, 393–404. [Google Scholar]
  4. Harman, I.T.; Drake, V.A. Insect monitoring radar: Analytical time-domain algorithm for retrieving trajectory and target parameters. Comput. Electron. Agric. 2004, 43, 23–41. [Google Scholar] [CrossRef]
  5. Wang, H.K.; Drake, V.A. Insect monitoring radar: Retrieval of wingbeat information from conical-scan observation data. Comput. Electron. Agric. 2004, 43, 209–222. [Google Scholar] [CrossRef]
  6. Cui, K.; Hu, C.; Wang, R.; Sui, Y.; Mao, H.; Li, H. Deep-learning-based extraction of the animal migration patterns from weather radar images. Sci. China Inf. Sci. 2020, 63, 10. [Google Scholar] [CrossRef] [Green Version]
  7. Chao, Z.; Wang, R.; Hu, C. Equivalent point estimation for small target groups tracking based on MLE. Sci. China Inf. Sci. 2019, 63, 8. [Google Scholar]
  8. Hu, C.; Wang, Y.; Wang, R.; Zhang, T.; Cai, J.; Liu, M. An improved radar detection and tracking method for small UAV under clutter environment. Sci. China Inf. Sci. 2019, 62, 29306. [Google Scholar] [CrossRef] [Green Version]
  9. Mao, H.; Wang, R.; Hu, C.; Yang, J. Fully polarimetric radar observing insects flight. In Proceedings of the IEEE International Conference on Signal, Information and Data Processing (ICSIDP), Chongqing, China, 11–13 December 2019; pp. 1–4. [Google Scholar]
  10. Hu, C.; Li, W.; Wang, R.; Long, T.; Drake, V.A. Discrimination of parallel and perpendicular insects based on relative phase of scattering matrix eigenvalues. IEEE Trans. Geosci. Remote Sens. 2020, 58, 3927–3940. [Google Scholar] [CrossRef]
  11. Long, T.; Hu, C.; Wang, R.; Zhang, T.; Zeng, T. Entomological radar overview: System and signal processing. IEEE Aerosp. Electron. Syst. Mag. 2020, 35, 20–32. [Google Scholar] [CrossRef]
  12. Hu, C.; Li, W.; Wang, R.; Liu, C.; Zhang, T.; Li, W. Accurate insect orientation extraction based on polarization scattering matrix estimation. IEEE Geosci. Remote Sens. Lett. 2017, 14, 1–5. [Google Scholar] [CrossRef]
  13. Hu, C.; Li, W.; Wang, R.; Long, T.; Liu, C.; Drake, V.A. Insect biological parameter estimation based on the invariant target parameters of the scattering matrix. IEEE Trans. Geosci. Remote Sens. 2019, 57, 6212–6225. [Google Scholar] [CrossRef]
  14. Drake, V.A.; Chapman, J.W.; Lim, K.S.; Reynolds, D.R.; Riley, J.R.; Smith, A.D. Ventral-aspect radar cross sections and polarization patterns of insects at X band and their relation to size and form. Int. J. Remote Sens. 2017, 38, 5022–5044. [Google Scholar] [CrossRef]
  15. Yu, T.; Wang, R.; Li, M.; Hu, C. Research on the design and calibration of wideband fully polarized vertical insect radar. J. Signal Process. 2021, 37, 222–233. [Google Scholar]
  16. Giuli, D. Polarization diversity in radars. Proc. IEEE 1986, 74, 245–269. [Google Scholar] [CrossRef]
  17. Fujita, M.; Masuda, T.; Fujino, Y.; Satake, M. Polarimetric calibration of the SIR-C C-band channel using active radar calibrators and polarization selective dihedrals. IEEE Trans. Geosci. Remote Sens. 1998, 36, 1872–1878. [Google Scholar] [CrossRef]
  18. Sarabandi, K.; Ulaby, F.T. A convenient technique for polarimetric calibration of single-antenna radar systems. IEEE Trans. Geosci. Remote Sens. 1990, 28, 1022–1033. [Google Scholar] [CrossRef]
  19. Boerner, W.M. Recent advances in polarimetry and polarimetric interferometry. In Proceedings of the IEEE International Symposium on Geoscience and Remote Sensing, Denver, CO, USA, 31 July–4 August 2006; pp. 49–51. [Google Scholar]
  20. Wiesbeck, W.; Riegger, S. A complete error model for free space polarimetric measurements. IEEE Trans. Antennas Propag. 1991, 39, 1105–1111. [Google Scholar] [CrossRef]
  21. Whitt, M.W.; Ulaby, F.T.; Polatin, P.; Liepa, V.V. A general polarimetric radar calibration technique. IEEE Trans. Antennas Propag. 1991, 39, 62–67. [Google Scholar] [CrossRef]
Figure 1. Block diagram of a fully polarimetric entomological radar system. The abbreviations in the figure are as follows: FM, frequency modulation; RF, radio frequency; ADC, analog-to-digital converter; LO, radio frequency; DC, down-conversion.
Figure 1. Block diagram of a fully polarimetric entomological radar system. The abbreviations in the figure are as follows: FM, frequency modulation; RF, radio frequency; ADC, analog-to-digital converter; LO, radio frequency; DC, down-conversion.
Remotesensing 14 01551 g001
Figure 2. Schematic diagram of PCTAR. (a) Rotate the antenna; (b) Rotate the calibrator.
Figure 2. Schematic diagram of PCTAR. (a) Rotate the antenna; (b) Rotate the calibrator.
Remotesensing 14 01551 g002
Figure 3. Relationships between SNRs and expectation, E, and standard deviation, δ , of the amplitude and phase errors of different polarization elements for PCTAR. (a) Expectations, E a m p , of amplitude errors; (b) standard deviations, δ a m p , of amplitude errors; (c) expectations, E p h a , of phase errors; (d) standard deviations, δ p h a , of phase errors.
Figure 3. Relationships between SNRs and expectation, E, and standard deviation, δ , of the amplitude and phase errors of different polarization elements for PCTAR. (a) Expectations, E a m p , of amplitude errors; (b) standard deviations, δ a m p , of amplitude errors; (c) expectations, E p h a , of phase errors; (d) standard deviations, δ p h a , of phase errors.
Remotesensing 14 01551 g003
Figure 4. The effects of beam alignment errors on the calibration precision of co-polarization channel (VV). (a) Expectation of amplitude errors; (b) standard deviation of amplitude errors; (c) expectation of phase errors; (d) standard deviation of phase errors.
Figure 4. The effects of beam alignment errors on the calibration precision of co-polarization channel (VV). (a) Expectation of amplitude errors; (b) standard deviation of amplitude errors; (c) expectation of phase errors; (d) standard deviation of phase errors.
Remotesensing 14 01551 g004
Figure 5. The effects of beam alignment errors on the calibration precision of the cross-polarization channel (HV). (a) Expectation of amplitude errors; (b) standard deviation of amplitude errors; (c) expectation of phase errors; (d) standard deviation of phase errors.
Figure 5. The effects of beam alignment errors on the calibration precision of the cross-polarization channel (HV). (a) Expectation of amplitude errors; (b) standard deviation of amplitude errors; (c) expectation of phase errors; (d) standard deviation of phase errors.
Remotesensing 14 01551 g005
Figure 6. Relationships between SNRs and expectation, E, and standard deviation, δ , of the amplitude and phase errors of different polarization elements for STCT. (a) Expectations, E a m p , of amplitude errors; (b) standard deviations, δ a m p , of amplitude errors; (c) expectations, E p h a , of phase errors; (d) standard deviations, δ p h a , of phase errors.
Figure 6. Relationships between SNRs and expectation, E, and standard deviation, δ , of the amplitude and phase errors of different polarization elements for STCT. (a) Expectations, E a m p , of amplitude errors; (b) standard deviations, δ a m p , of amplitude errors; (c) expectations, E p h a , of phase errors; (d) standard deviations, δ p h a , of phase errors.
Remotesensing 14 01551 g006
Figure 7. The differences of the expectation and standard deviation of the amplitude and phase errors of STCT and PCTAR at different SNRs and different polarization channels. (a) Differences of amplitude error expectation Δ E a m p ; (b) differences of standard deviation of amplitude error Δ δ a m p ; (c) differences of phase error expectation Δ E p h a ; (d) differences of standard deviation of phase error Δ δ p h a .
Figure 7. The differences of the expectation and standard deviation of the amplitude and phase errors of STCT and PCTAR at different SNRs and different polarization channels. (a) Differences of amplitude error expectation Δ E a m p ; (b) differences of standard deviation of amplitude error Δ δ a m p ; (c) differences of phase error expectation Δ E p h a ; (d) differences of standard deviation of phase error Δ δ p h a .
Remotesensing 14 01551 g007
Figure 8. Schematic diagram of the experiment.
Figure 8. Schematic diagram of the experiment.
Remotesensing 14 01551 g008
Figure 9. Field experiment for polarimetric calibration (the use of UAVs had been reported to the local air traffic control department).
Figure 9. Field experiment for polarimetric calibration (the use of UAVs had been reported to the local air traffic control department).
Remotesensing 14 01551 g009
Figure 10. Calibrators and targets for validation. (a) Target Ⅰ: sphere with 122mm diameters; (b) Target Ⅱ: wrench with 120mm length; (c) Target Ⅲ: sphere with 60mm diameters; (d) Target Ⅳ: bollworm with 15mm length.
Figure 10. Calibrators and targets for validation. (a) Target Ⅰ: sphere with 122mm diameters; (b) Target Ⅱ: wrench with 120mm length; (c) Target Ⅲ: sphere with 60mm diameters; (d) Target Ⅳ: bollworm with 15mm length.
Remotesensing 14 01551 g010
Figure 11. Polarization pattern of target Ⅲ.
Figure 11. Polarization pattern of target Ⅲ.
Remotesensing 14 01551 g011
Figure 12. Polarization pattern of target IV.
Figure 12. Polarization pattern of target IV.
Remotesensing 14 01551 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yu, T.; Li, M.; Li, W.; Mao, H.; Wang, R.; Hu, C.; Long, T. Polarimetric Calibration Technique for a Fully Polarimetric Entomological Radar Based on Antenna Rotation. Remote Sens. 2022, 14, 1551. https://doi.org/10.3390/rs14071551

AMA Style

Yu T, Li M, Li W, Mao H, Wang R, Hu C, Long T. Polarimetric Calibration Technique for a Fully Polarimetric Entomological Radar Based on Antenna Rotation. Remote Sensing. 2022; 14(7):1551. https://doi.org/10.3390/rs14071551

Chicago/Turabian Style

Yu, Teng, Muyang Li, Weidong Li, Huafeng Mao, Rui Wang, Cheng Hu, and Teng Long. 2022. "Polarimetric Calibration Technique for a Fully Polarimetric Entomological Radar Based on Antenna Rotation" Remote Sensing 14, no. 7: 1551. https://doi.org/10.3390/rs14071551

APA Style

Yu, T., Li, M., Li, W., Mao, H., Wang, R., Hu, C., & Long, T. (2022). Polarimetric Calibration Technique for a Fully Polarimetric Entomological Radar Based on Antenna Rotation. Remote Sensing, 14(7), 1551. https://doi.org/10.3390/rs14071551

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